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Gasoline  is  a  significant  source  of  groundwater 
contamination  in  Florida.   This  results  from  the  large 
numbers  of  gasoline  storage  tanks,  high  rainfall,  reliance 
on  groundwater-based  potable  water  supplies  and  the 
hydrogeology  of  Florida.   Sorption,  biodegradation  and 
hydrolysis  of  dissolved  aromatic  hydrocarbons  (all  isomers 
of  CgHg-CgH,.)  were  determined  in  multicomponent  experiments 
with  natural  aquifer  materials  under  saturated  conditions. 
Hydrogen  peroxide,  air,  oxygen  gas  and  ammonium  chloride 
treatments  were  evaluated  as  methods  to  enhance  microbial 
degradation  of  aromatic  hydrocarbons.   The  solutes  and 
sorbents  were  from  a  gasoline  contaminated  aquifer  in 
central  Florida.   This  site  was  typical  of  sandy  surficial 

xiii 


aquifers  in  Florida  with  a  low  organic  carbon  content 
(0.015%).   The  aquifer  was  composed  primarily  of  fine  to 
medium  grained  sands. 

Hydrolysis  was  not  a  significant  removal  mechanism  for 
the  selected  aromatic  solutes.   Equilibrium  batch  isotherms 
and  column  studies  determined  sorption  coefficients  for 
aromatic  solutes  ranging  between  0.045  and  0.1  with 
retardation  values  between  1.36  and  2.40.   Column 
breakthrough  curves  exhibited  minimal  effects  of  adsorption 
non-equilibrium.   Sorption  isotherms  were  linear  through  the 
concentration  range  tested  and  no  significant  hysteresis  was 
noted.   Partitioning  and  surface  dependent  adsorption  were 

evaluated  by  regression  of  column  K    data  with  literature 

1  oc 

values  of  K  ,    water  solubility  and  first  order  connectivity 

ow  J 

indices.   No  single  model  fully  described  the  sorption 
process,  and  the  sorption  mechanism  appeared  to  be  a 
combination  of  several  processes.   Competitive  solute 
interactions  were  not  shown  to  be  significant. 

Column  biodegradat ion  experiments  with  acclimated 
microorganisms  were  performed  at  flow  velocities  close  to 
those  from  the  contaminated  aquifer.   Half  lives  ranged  from 
0.940  hr  for  benzene  to  0.086  hr  for  n-propy lbenzene  at 
0.680  cm/min.   Branched  aromatic  solutes  were  more  easily 
degraded  in  column  studies. 

Batch  studies  demonstrated  the  ability  of  field 
microbes  to  degrade  aromatic  hydrocarbons  to  less  than  0.5 
ug/L  given  sufficient  oxygen.   Microbes  were  not  phosphorus 


xiv 


or  nitrogen  limited.   Hydrogen  peroxide  increased  dissolved 
oxygen/  but  did  not  lead  to  increased  hydrocarbon  removal  in 
lab  studies.   Ammonium  chloride  produced  nitrifying 
conditions.   Oxygen  augmentation  with  air  and  oxygen  gas  was 
shown  to  enhance  biological  removal  of  aromatic 
hydrocarbons  . 
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CHAPTER  I 
INTRODUCTION 


Groundwater  contamination  is  a  topic  of  great 
scientific  interest  and  public  concern.   Groundwater 
provides  approximately  100  million  people  with  potable  water 
in  the  United  States  ( Hoag  and  Marley,  1986)  and  nearly 
every  state  contains  some  number  of  contaminated  wells 
(Barbash  and  Roberts/  1986).   The  sources  of  groundwater 
contamination  are  numerous.   These  include  seepage  from 
lagoons  and  impoundments/  landfills,  agricultural  and 
silvicultural  practices,  accidental  spills  and  leaking 
storage  tanks  and  transfer  equipment. 

Gasoline  and  petroleum  products  are  some  of  the  most 
common  groundwater  pollutants.   The  potential  magnitude  of 
this  problem  is  evident  from  the  volume  of  petroleum  used  in 
the  United  States.   Approximately  110  billion  gallons  of 
motor  fuels  are  stored  underground  each  year,  in  an 
estimated  1.4  million  underground  storage  tanks,  85%  of 
which  are  unprotected  steel  tanks  with  finite  lifetimes 
(Hoag  and  Marley,  1986).   It  is  expected  that  10  to  30%  of 
these  tanks  may  leak  ( Dowd ,  1984). 

Gasoline  contamination  of  groundwater  in  Florida  is  a 
particularly  serious  problem.   This  results  from  the 
confluence  of  three  factors:  the  large  number  of  petroleum 
storage  tanks  in  the  state,  the  reliance  on  groundwater 


based  potable  water  supplies,  and  the  hydrogeology  of 
Florida . 

The  major  sources  of  petroleum  contamination  in  Florida 
are  leaking  storage  tanks  and  pipes.   The  high  water  table 
in  the" state  leads  to  conditions  favorable  for  corrosion. 
As  of  February  1986  there  were  455  known  storage  tank 
incidents  resulting  in  368  cases  of  groundwater 
contamination.   The  total  volume  of  spilled  gasoline  exceeds 
4.2  million  gallons  (FLDER,  1986).   The  remaining  60,000 
petroleum  storage  tanks  in  the  state  provide  potential 
sources  for  future  groundwater  pollution. 

These  sources  of  contamination  are  particularly 
significant  owing  to  the  importance  of  groundwater  in 
Florida.  Groundwater  withdrawal  for  potable  water  use  is 
approximately  1400  Mgal/d,  comprising  87%  of  public  water 
and  94%  of  rural  water  supplies.   It  is  noteworthy  that 
nearly  2  million  residents  drink  untreated  water  from 
shallow  private  wells  which  are  particularly  prone  to 
contamination  from  underground  storage  tanks  (Fernald  and 
Patton,  1984). 

Hydrogeology  is  the  third  factor  which  contributes  to 
the  sensitivity  of  Florida's  water  supplies  to  gasoline 
contamination.   Most  of  the  potable  water  aquifers  are 
surficial  or  intermediate  in  depth,  and  are  susceptible  to 
contamination.   In  addition,  the  generally  porous  nature  of 
top  soil  in  the  state  enhances  pollutant  transport  to  the 
underlying  aquifers.   Most  soils  in  Florida  are  sandy  loam, 


sandy  clay  and  sandy  clay  loams,  all  of  which  are  noted  for 
their  relatively  high  permeabilities  (Fernald  and  Patton, 
1984).   The  sandy  deposits  of  the  Pliocene  and  Pleistocene 
ages  common  to  Florida  are  also  marked  by  low  organic  carbon 
and  clay  content  (Fetter,  1980),  resulting  in  high 
permeability  and  low  sorptive  capacity. 

Given  the  magnitude  of  this  problem,  the  transport  and 
environmental  reactions  of  dissolved  gasoline  components  in 
shallow  sandy  aquifers  is  an  important  area  of  study.   This 
is  particularly  true  of  the  aromatic  constituents  of 
gasoline,  owing  to  their  toxicity,  concentrations  in 
gasoline  and  their  high  aqueous  solubility.   Gasoline 
products  which  are  released  into  the  vadose  zone  travel 
downward  under  the  influence  of  capillary  and  gravitational 
forces.   When  the  sorptive  capacity  of  the  soil  is  exceeded, 
gasoline  moves  onto  the  groundwater  table,  where  it  spreads 

laterally  across  the  top  of  the  saturated  zone.   This  is  a 

3 
result  of  the  density  of  gasoline  (0.7-0.75  g/cm  ). 

Gasoline  components  partition  into  the  water  to  the  extent 

of  their  water  solubility,  and  move  in  the  direction  of  the 

water  table  gradient. 

Physical,  chemical  and  biological  factors  must  all  be 

considered  in  the  determination  of  the  fate  and  transport  of 

dissolved  gasoline  hydrocarbons  in  groundwater.   The 

interaction  of  these  factors  may  be  conveniently  examined  in 

the  context  of  a  generalized  mass  transport  equation.   A 

one-dimensional  form  of  this  equation  is 
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3C/3t 


Dh  (   3C/  3x  )  -  v  (  3C/  3x  )  -  p/  9(  3S/  3t )  -  Q . 


[1.1] 


where  C 
S 
6 

t 
x 

P 
Dh 

V 

Qi 


solution  phase  concentration  of  solute  (ug/L) 

adsorbed  phase  concentration  of  solute  (ng/g) 

volumetric  water  content  (mL/cm  ) 

time  (rain) 

horizontal  distance  (cm) 

bulk  density  (g/cm  ) 

hydrodynamic  dispersion  coefficient  (cm  /rain) 

average  pore  water  velocity  (cm/min) 

degradation  rate  (rain   ) 


The  components  of  this  equation  are  convection/ 
dispersion/  sorption  and  degradation  terms.   Convection 
describes  the  movement  of  a  dissolved  contaminant  with  the 
groundwater.   Dispersion  describes  the  spreading  of  the 
solutes  during  flow  through  the  aquifer  material.   Sorption 
terms  account  for  the  retardation  of  the  dissolved  solutes 
by  interaction  with  the  aquifer  matrix,  and  degradation 
terms  evaluate  the  removal  or  transformation  of  the 
contaminants.   Research  has  shown  that  sorption  and 
biological  degradation  are  the  major  attenuation  mechanisms 
for  organic  solutes  in  soils  and  groundwater  (Woodburn, 
1985). 

Mathematical  models  based  on  such  equations  are 
important  tools  for  the  prediction  of  contaminant  movement 
(Pinder,  1984).   However,  the  adequacy  of  these  predictions 
is  directly  related  to  a  knowledgeable  and  accurate 
quantification  of  the  processes  involved  (MacKay  et  al . , 
1985).   The  remediation  of  groundwater  contamination  also 


requires  detailed  and  usually  site  specific  data  for  these 
processes . 

This  dissertation  presents  a  detailed  investigation  of 
the  attenuation  of  selected  dissolved  aromatic  gasoline 
hydrocarbons  in  a  typical  sandy  surficial  aquifer  in 
Florida,  using  a  variety  of  batch  and  column  techniques. 
Sorption  coefficients  and  biological  and  abiotic  degradation 
rates  from  laboratory  studies  are  presented.   These  studies 
simulated  conditions  at  a  contaminated  field  site  and 
included  experiments  to  assess  various  treatment 
alternatives.   Actual  field  data  are  summarized,  and 
compared  to  the  laboratory  data. 

The  improved  understanding  obtained  from  the  collection 
and  analysis  of  these  data  should  aid  in  the  formation  and 
improved  use  of  predictive  models  describing  the  movement 
and  reaction  rates  of  water  soluble  components  of  gasoline 
in  shallow  groundwater  systems  and  aid  in  the  selection  of 
appropriate  groundwater  reclamation  technologies. 


CHAPTER  II 
OBJECTIVES 


The  main  objectives  of  this  study  were: 


(1)  To  evaluate  the  sorption  coefficients  for  12 
selected  aromatic  hydrocarbons  found  in  water  at  the  Lake 
Alfred  research  site/  employing  batch  isotherms  and  soil 
columns; 

(2)  To  determine  the  rates  of  hydrolysis  of  the  12 
selected  aromatic  hydrocarbons; 

(3)  To  determine  the  rates  of  biodegradation  of  the  12 
selected  aromatic  hydrocarbons  under  simulated  field 
conditions,  and  after  treatment  with  hydrogen  peroxide, 
oxygen  gas  and  ammonium  chloride; 

(4)  To  determine  the  most  appropriate  predictive  model 
for  sorption  of  the  12  selected  aromatic  hydrocarbons  in  a 
sandy  surficial  aquifer  in  Florida; 

(5)  To  correlate  molecular  properties  of  the  selected 
aromatic  hydrocarbons  with  sorptive  and  biological 
parameters; 

(6)  To  evaluate  field  data  based  on  the  laboratory 
measurements  of  sorption  and  biodegradation  and 

(7)  To  extrapolate  the  laboratory  data  for  application 
of  aquifer  remediation  practices. 


CHAPTER  III 
LITERATURE  REVIEW 


3.1   Introduction 


This  chapter  presents  a  review  of  the  pertinent 
literature  for  the  reactions  of  gasoline  derived  aromatic 
hydrocarbons  in  groundwater.   The  major  areas  of  discussion 
are  the  environmental  effects  of  gasoline  contamination,  the 
use  of  advection-dispersion  transport  models,  the  sorption 
of  aromatic  compounds  to  aquifer  materials,  and  the 
biodegradation  of  aromatic  compounds  in  groundwater  systems. 

3. 2   Environmental  Effects  of  Gasoline  Contamination 

Gasoline  is  a  complex  mixture  of  many  hydrocarbon 

compounds.   A  typical  gasoline  contains  between  150-250 

identifiable  hydrocarbon  components  (Sanders  and  Maynard, 

1968)  consisting  of  alkane,  alkene,  aromatic  and  napthene 

hydrocarbons.   Automobile  gasolines  are  comprised  of  Ce-C,_ 

5   12 

hydrocarbons  with  boiling  points  in  the  range  32-210  C.   Un- 
leaded gasolines  contain  greater, concentrations  of  aromatic 
hydrocarbons  to  provide  for  anti-knock  protection  and 
branched  hydrocarbons  to  increase  octane  ratings  (Moore  and 
Moore,  1976). 
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Despite  the  large  number  of  hydrocarbons  comprising 

gasoline,  much  of  the  environmental  concern  focuses  on  the 

water  soluble  components  of  gasoline/  particularly  the 

single  ring  aromatic  compounds.   These  compounds  are  of 

concern  based  on  their  toxicity,  aqueous  solubility  and 

concentration  in  gasoline  (Barker  and  Patrick,  1985).   Acute 

toxicity  is  associated  with  the  water  soluble  fraction  of 

oils  (Blumer  et  al . ,  1973)  and  the  major  components  of  the 

water  soluble  fraction  are  aromatic  (Coleman  et  al . ,  1984). 

Data  from  the  work  of  Coleman  et  al.  (1984)  showed  that 

although  aromatic  components  made  up  only  50%  of  the 

unleaded  gasoline  product  in  their  study,  87-95%  of  the 

components  in  the  water  soluble  fraction  were  aromatic. 

Thus  in  a  spill  situation  a  significant  amount  of  the 

contaminants  in  the  water  phase  will  be  aromatic.   Selected 

physical  properties  of  the  compounds  used  in  this  study  are 
listed  in  Table  3-1. 

The  health  effects  from  the  use  of  gasoline 

contaminated  water  may  be  significant.   Benzene  is  a 

carcinogen  in  rats  and  mice  and  exposure  is  linked  with 

leukemia  (USPHS,  1981).   The  maximum  contaminant  level  for 

benzene  in  community  drinking  water  supplies  is  1  ppb  in 

Florida.   Toluene,  ethylbenzene  and  m-xylene  affect  the 

central  nervous  system  (Windholtz,  1976).  Unleaded  gasoline 

induces  renal  and  hepatocellular  carcinomas  in  rats  and  the 

use  of  petroleum  contaminated  water  can  produce  elevated 

levels  of  indoor  air  pollutants  allowing  chronic  exposure  to 
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hydrocarbons  (Shehata,  1985).  Dermal  absorption  of  volatile 
organic  contaminants  from  gasoline  may  also  be  a  significant 
exposure  (Brown  et  al . ,  1984). 

Fire  and  explosion  hazards  are  also  a  risk  factor  in 
the  release  of  gasoline  to  the  environment.   Volatilization 
and  subsequent  gas  phase  transport  of  hydrocarbons  in  the 
unsaturated  zone  have  destroyed  buildings  (Hoag  and  Marley, 
1986). 

3.3  Convecti ve-Dispersi ve  Models 


The  cogent  evaluation  of  contaminant  plumes,  remedial 
action  alternatives,  and  risk  assessment  for  organic 
compounds  in  groundwater  requires  a  thorough  understanding 
of  the  behavior  of  these  contaminants  in  groundwater 
systems.   This  includes  an  assessment  and  quantification  of 
the  relevant  processes  which  influence  their  fate  and 
transport  (Miller  and  Weber,  1984). 

The  interaction  of  these  processes  may  be  examined  in 
the  context  of  convecti ve-dispersive  models.   These  models 
have  been  reviewed   (Anderson,  1979,  Freeze  and  Cherry, 
1979)  and  are  marked  by  their  computational  simplicity, 
reasonable  data  requirements  and  sufficiently  accurate 
output  (Roberts  et  al . ,  1985).  Although  the  adequacy  of 
convective-dispersi ve  models  for  describing  solute  transport 
has  been  questioned  (Anderson,  1979,  Smith  and  Schwartz, 
1980),  particularly  with  regard  to  dispersivity 
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approximations,  these  models  provide  a  convenient  framework 
for  understanding  the  transport  of  dissolved  solutes  in 
groundwater . 

The  general  form  of  the  solute  transport  equation  under 
saturated  flow  conditions  is  given  by  Bear  (1979).   A  one - 
dimensional  form  of  this  equation  for  conservative 
contaminants  under  steady  flow  conditions  is 


3C/3t  -  Dh(  32C/9x2)-  v(  3C/3x) 


[3.1] 


where  C  = 
S  = 

K 

X  = 
V  = 


solution  phase  concentration  of  solute  (ug/L) 

sorbed  phase  concentration  of  solute  (ng/g) 

hydrodymanic  dispersion  coefficient  (cm  /min) 

time  (min) 

horizontal  distance  (cm) 

average  pore  water  velocity  (cm/rain) 


The  major  components  of  this  equation  are  convection  (bulk 
flow)  and  dispersion  (deviation  from  bulk  flow).   A  brief 
discussion  of  dispersion  follows,  with  reference  to 
extrapolation  of  laboratory  data  to  field  scale 
applications  . 

The  hydrodynamic  dispersion  coefficient  describes  the 
spreading  of  a  solute  as  it  moves  through  porous  media. 
Hydrodynamic  dispersion  (D  )  is  the  sum  of  mechanical 
dispersion,  caused  by  differences  in  water  velocity  through 
sinuous  and  tortuous  pores,  and  molecular  diffusion  (Biggar 
and  Nielsen,  1962).   Dispersion  values  reflect  the 
heterogeneity  of  the  aquifer  material.   Dispersion  is 
usually  determined  by  measuring  the  breakthrough  of  a 
conservative  tracer  such  as  chloride  or  tritiated  water. 
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The  physical  and  mathematical  relationships  of  water  and 
solute  transport  were  reviewed  by  Davidson  et  al .  (1983). 
Solute  dispersion  was  noted  to  occur  because  of  macroscale 
spatial  changes  in  the  direction  and  magnitude  of  water 
flow.   The  continuum  approach  to  mathematically  describe 
water  and  solute  transport  in  laboratory  soil  columns  was 
shown  to  be  reasonably  successful. 

In  practice,  laboratory  measurements  and  theory  may  be 
of  little  value  in  predicting  dispersion  in  natural 
aquifers.   Laboratory  columns  give  dispersivity  estimates  on 
the  order  of  centimeters,  whereas  field  scale  dispersion  is 
usually  in  meters  (Bedient  et  al . ,  1985).   This  is  a  result 
of  the  greater  heterogeneity  of  a  field  site  versus  a  small 
homogeneous  laboratory  column.   A  solution  for  equation 
[3.1]  for  a  finite  column  using  dimensionless  variables  was 
presented  by  Brenner  (1962).   The  dimensionless  Peclet 
number  (P  )  was  used  as  a  measure  of  dispersion: 

Pe  =  vL/4Dh  [3.2] 

where  v  is  pore  water  velocity  (cm/min),  L  is  the  length 

(cm)  of  the  soil  column  and  D   is  the  hydrodynamic 

2 
dispersion  coefficient  (cm  /rain).   For  values  of  P   >  100 

e 

dispersion  is  assumed  negligible.   Values  of  P   <  10 

e 

generally  indicate  complete  mixing.   Boundary  conditions  for 
displacement  experiments  through  short  laboratory  columns 
were  reviewed  by  van  Genuchten  and  Parker  (1984).   The 
solution  of  Brenner  (1962)  was  shown  to  correctly  conserve 
mass  in  finite  laboratory  soil  columns,  based  on  mass 
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balance  considerations.   For  a  flux  type  inlet  boundary 
condition  (flowing  concentrations),  Brenner's  solution  was 
applicable  provided  the  column  Peclet  number  was  not  much 
less  than  five.   The  solution  of  Lapidus  and  Amundson  (1952) 
was  recommended  to  evaluate  flux  averaged  concentrations  in 
finite  laboratory  columns  or  semi  infinite  field  profiles. 

3.4  Sorption  of  Aromatic  Compounds 


Sorption  is  a  major  mechanism  in  the  attenuation  of 

organic  solutes  in  the  saturated  zone.   Solutes 

differentially  sorb  onto  aquifer  materials  and  thus  are 

retarded  in  their  movement  through  the  subsurface,  resulting 

in  a  chromatographic  like  separation  of  the  soluble 

constituents  of  a  plume,  with  groundwater  as  the  mobile 

phase . 

Sorption  describes  the  transfer  of  solutes  from  a 

liquid  phase  to  a  solid  phase  (Miller  and  Weber,  1984).   In 
this  literature  review  the  liquid  phase  is  assumed  to  be 

water,  containing  solubilized  organic  solutes  and  the  solid 

phase  is  the  aquifer  material  under  saturated,  steady  flow 

conditions.   Sorption  is  influenced  by  the  physical  and 

chemical  characteristics  of  the  aquifer  (ie.,  soil  type, 

fraction  of  organic  carbon),  and  the  solute  ( ie . , 

solubility,  volatility,  density). 

Although  sorption  is  a  major  component  in  the 

attenuation  of  solutes  in  the  subsurface,  the  fundamental 


14 

processes  of  solute-soil  interaction  and  the  thermodynamics 
of  this  process  are  not  completely  characterized.  Therefore, 
sorption  is  used  in  this  study  as  a  generic  term  to  describe 
solute  retention  ( ie .  uptake  of  solute),  regardless  of 
whether  the  process  is  one  of  adsorption,  absorption  or 
partitioning  (Woodburn,  1985).   Desorption  is  used  here  to 
describe  solute  removal  from  the  solid  phase. 
3.4.1   Sorption  Processes 

The  attractive  forces  acting  to  effect  sorption  of 
hydrophobic  compounds  onto  natural  sorbents  were  reviewed  by 
Voice  and  Weber  (1983).   The  major  theory  is  discussed 
below. 

Bonding  forces  in  sorption  may  be  both  physical  and 
chemical,  though  both  are  basically  electrostatic  in  nature. 
Physical  sorption  results  from  Van  der  Waals  forces.   The 
strength  of  these  interactions  is  generally  on  the  order  of 
1-2  Kcal/mole.   These  energies  may  be  augmented  by  a 
thermodynamic  gradient  driving  hydrophobic  molecules  out  of 
solution.   This  is  based  on  entropic  considerations 
(solvophobic  theory). 

Chemical  sorption  is  the  interaction  between  specific 
sites  of  the  sorbent  and  individual  solute  molecules.   This 
approximates  a  true  chemical  bond  with  heats  of  adsorption 
between  15-50  Kcal/mole.   Voice  and  Weber  (1983)  point  out 
that  it  is  difficult  to  assess  the  importance  of  each  type 
of  bonding.   The  heterogeneous  nature  of  natural  sorbent 
materials  is  largely  unknown,  and  sorption  processes 
probably  involve  all  types  of  interactions. 
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3.4.2   Sorption  Equilibria 

Two  experimental  techniques  are  widely  used  to  evaluate 
the   gS/  9t   term  in  equation  [1.1],   These  are  batch 
equilibrium  and  soil  column  methods.   Batch  studies  allow 
the  evaluation  of  the  linearity  of  the  sorption  isotherm  and 
their  use  is  well  documented  ( Schwarzenbach  and  Westall, 
1981,  Chiou  et  al . ,  1979).   The  most  widely  used  models  to 
describe  sorption  equilibria  in  groundwater  systems  are  the 
linear  [3.3]  and  Freundlich  models  [3.4]  (Miller  and  Weber, 
1984)  : 

[3.3] 
(n  <  1)  [3.4] 

where  S  (ug/g)  and  C  (ug/L)  are  the  adsorbed  phase  and 
solution  phase  concentrations  respectively  at  equilibrium, 
K   (L/g)  is  the  linear  sorption  coefficient,  K_  (L/g)  is  the 
Freundlich  sorption  coefficient  (both  K,  and  Kf  indicating 
sorption  capacity)  and  n  is  an  empirical  constant 
(indicating  sorption  intensity).   The  linear  model  is  in 
effect,  a  special  case  of  the  Freundlich  model  where  n=l. 
The  Freundlich  equation  is  often  linearized  (log 
transformed)  to  facilitate  calculation  of  variables  Kf  and  n 
in  batch  studies: 


S  =  K,  *  C 
d 

S  =  Kf  *  C1 


log  S  =  n  *  log  C  +  Log  K 


[3.5] 


In  column  studies  K,  is  evaluated  through  the  retardation 
factor  (R).   The  mass  transport  equation  for  reactive 
solutes  under  steady  flow  is  described  by  equation  [3.6]: 
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ac/at  +  p/e  a  s/at    = 


2     2 

D,   9C  /3x   -  v  3  C/8x 


[3.6] 


where  p  is  the  bulk  density,  0  is  the  volumetric  water 
content  and  S  is  the  sorbed  phase  concentration.   Note  that 
equation  [3.6]  is  equivalent  to  equation  [3.1]  with  the 
addition  of  the  sorption  term   33/  St.   Assuming  linear , 
reversible  sorption,  the  sorbed  concentration  of  a  solute  is 
related  to  the  aqueous  concentration  of  the  solute  by  the 
relationship : 

as/at    =    Kd  ac/at  [3.7] 


Substitution  for  9 S/  at  in  equation  [3.6]  with  equation 
[3.7]  yields  the  relationship: 

ac/at  +  k    ac/at   (p/e)   =     Dh    a2c/ax2  -  v  ac/ax  [3.8] 

After  separation  of  variables  equation  [3.8]  becomes 

ac/at   [  l  +  p  K./.e]     =     Dh     s2c/ax2     -     v  ac/ax  [3.9] 


r  =  i  +  p  Kd/e 


and  by  defining  the  retardation  factor  (R)  as 

[3.10] 

substitution  of  equation  [3.10]  into  [3.9]  results  in  the 
incorporation  of  the  retardation  factor  (R)  into  the  mass 
transport  equation  for  solute  transport  under  saturated 
steady  flow  conditions: 

[3.11] 

Analysis  of  equation  [3.10]  indicates  that  the  value  of 
R  is  largely  dependent  on  K   for  a  homogeneous  aquifer 
system  or  laboratory  column.   Determination  of  R  from  soil 


r   ac/at     =     Dh    £c/ ax2  -v  ac/ax 
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column  studies  leads  to  the  evaluation  of  K   from  equation 
[3.10].   Nkedi-Kizza  et  al .  (1987)  compared  techniques  for 
the  calculation  of  R  from  soil  column  leaching  experiments 
and  from  batch  isotherm  experiments.   Values  of  R  calculated 
by  determining  the  area  above  the  breakthrough  curve  were 
shown  to  be  equivalent  to  R  values  calculated  by  using 
equation  [3.10]. 
3.4.3   Sorption  Estimators 

Recently,  approximation  methods  based  on  the  assumption 
of  partitioning  as  the  dominant  method  of  solute  interaction 
have  become  common  (Karickhoff  et  al . ,  1979,  Chiou  et  al . , 
1979,  Kenaga  and  Goring,  1980,  Chiou  et  al . ,  1983).   Their 
use  is  largely  a  result  of  the  time  and  difficulty  in  the 
accurate  measurement  of  sorption  coefficients  (K  ),  and  the 
general  lack  of  data  on  hydrocarbon  sorption  to 
environmental  sorbents  .   These  authors  note  a  correlation 
between  the  fractional  organic  carbon  content  of  the  sorbent 

material  (f   )  and  K , .   The  K,  normalized  to  f    of  the 
oc        a         a  oc 


sorbent  is  described  as  K  where: 

oc 

K  =  K,/f 

oc     d   oc 


[3.12] 


Values  of  Kqc  is  well  correlated  with  aqueous  solubility 
(WS)  (Chiou  et  al . ,  1979)  and  the  octanol-water  partition 
coefficient  (KQW)  (Karickhoff  et  al . ,  1979).   These  authors 
suggest  that  the  solute-sorbent  interaction  is  a 
partitioning  process  rather  than  an  interaction  between 
solute  and  the  mineral  surface.   Evidence  for  partitioning 
is  partially  supported  by  the  hydrophobic  character  of  soil 
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Log  K    =  a  *  Log  K    +  Log  f    +  b 

oc  3  ow      3   oc 

Log  K    =  c  *  Log  WS   +  Loq  f    +  d 

oc  J          3   oc 


organic  matter,  and  by  solvophobic  theory  (Rao  et  al .  , 

1985).  The  general  relationship  between  K    and  K    and  WS 

oc       ow 

takes  the  form  (Curtis  et  al . ,  1986): 

[3.13] 
[3.14] 

where  a,b,c,  and  d  result  from  regression  analysis  of 
laboratory  isotherm  data  and  depend  on  the  solute-sorbent 
system. 

However,  there  are  significant  limitations  on  the  use 
of  these  estimators,  and  the  basis  of  partitioning  as  a 
sorption  mechanism  is  questionable  (Milgelgrin  and  Gerstl, 
1983).   In  a  strict  sense,  these  relationships  hold  only  for 
those  compounds  and  sorbents  used  in  the  original  studies 
(i.e.,  these  are  empirical  relationships).   This  is 
reflected  in  orders  of  magnitude  variation  in  estimates  of 
sorption  using  these  relationships.   Application  of  the 
partitioning  models  may  not  be  appropriate  in  experimental 
systems  with  solutes  and  sorbents  which  are  different  from 
those  used  to  develop  these  models.   In  addition,  these 
equations  may  not  apply  at  organic  carbon  fractions  less 
than  0.1%  (Curtis  et  al . ,  1986).   Rao  and  Jessup  (1983) 
noted  that  the  use  of  Kqc  to  estimate  sorption  can  lead  to 
significant  errors  with  soils  with  very  low  (less  than  0.1%) 
organic  carbon  contents. 

Milgelgrin  and  Gerstl  (1983)  reviewed  the  evidence  for 
partitioning  and   noted  that  a  correlation  between  the 
organic  carbon  content  of  the  soil  and  sorption  was  not 
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universally  significant.   These  authors  cited  several 

studies  where  removal  of  organic  carbon  from  a  soil  actually 

increased  the  amount  of  sorption/  or  had  no  negative  effect 

on  the  sorption  values.   These  authors  suggested  that 

molecular  structure  of  the  solute  may  be  a  better  predictor 

of  sorption  to  sediments  than  water  solubility  or 

octanol/water  partition  coefficients.   This  results  from  the 

observation  that  with  a  relatively  rigid  adsorbing  surface, 

the  conformation  of  the  solute  molecule  will  greatly  affect 

its  adsorption  (i.e./  steric  effects),  but  not  its 

partitioning  between  an  organic  phase  and  water. 

Recently/  first  order  molecular  connectivity  indexes 

(  X)  were  shown  to  be  well  correlated  with  K    values 

oc 

(SabljiC/  1984,  Sabljic,  1987).   Molecular  connectivity  is 

described  as  a  quantitative  measure  of  the  area  occupied  by 

the  projection  of  the  non-hydrogen  skeleton  of  a  molecule. 

The  correlation  between  K    and  the  first  order  molecular 

oc 

connectivity  index  supports  the  contention  that  the  process 
of  soil  sorption  may  be  viewed  as  an  attractive  interaction 
between  two  planes,  with  the  magnitude  of  the  interaction 
directly  proportional  to  the  surface  area  of  the  molecule. 
This  suggests  that  the  soil  sorption  and  partitioning 
process  reflect  different  mechanisms.   An  accurate  model  of 
sorption  may  include  both  partitioning  and  surface  area 
dependent  affects. 


The  relationship  between  K    and   X  is  (Sabljic, 

oc  J 


1987) 
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Log  K    =  0.53  *   X  +  0.54 
oc 


[3.15] 


This  relationship  is  based  on  literature  values  of  K    from 

oc 

laboratory  experiments  with  72  compounds  covering  a  broad 
range  of  polarities  and  classes,  and  a  variety  of  sorbent 
systems.   The  correlation  coefficient  is  0.976  which 
explains  95.2%  of  the  variance. 
3.4.4   Desorption 

In  most  cases  sorption  is  considered  to  be  completely 
reversible,  that  is,  the  adsorption-desorption  isotherms  are 
reversible  and  single  valued.   However,  several 
investigators  report  desorptions  which  display  hysteresis  in 
batch  studies  (Bailey  and  White,  1970,  Boucher  and  Lee, 
1972,  Carringer  et  al . ,  1975,  DiToro  and  Horzempa,  1982). 
Van  Genuchten  et  al .  (1974)  found  that  the  exponent  for 
desorption  is  concentration  dependent,  and  described  the" 
hysteretic  behavior  by  using  separate  isotherm  equations  for 
sorption  and  desorption: 

[3.16] 
[3.17] 

where  subscripts  s  and  d  indicate  sorption  and  desorption 
respectively.   Hysteresis  in  column  studies  was  noted  by 
Schwarzenbach  and  Westall  (1981),  although  the  reaction  was 
termed  reversible,  since  all  the  solute  was  eventually 
eluted  from  the  column. 


S   =  K,   C 
s     ds   s 


Sd  "  Kdd  Cd 


ns 


nd 
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3.4.5  Sorption  Kinetics 

Hysteretical  behavior  may  actually  be  a  manifestation 
of  sorption-desorption  kinetics.   Rao  and  Jessup  (1983) 
noted  that  the  influence  of  non-singular  isotherms  (ie./ 
isotherms  which  display  hysteresis)  on  solute  movement  may 
be  less  significant  than  the  effects  of  sorption 
nonequilibria .   In  a  study  of  the  transport  of  pesticides  at 
high  concentrations/  Rao  and  Davidson  (1979)  noted  that  the 
position  of  an  adsorbed  solute  in  a  breakthrough  curve  was 
governed  by  the  nature  of  the  equilibrium  adsorption 
isotherm  equation,  whereas  the  shape  of  the  curve  was 
defined  by  the  kinetics  of  the  sorption-desorption  process. 

Sorption  reactions  between  hydrophobic  pollutants  and 
sediments  are  generally  rapid  and  not  rate  limited  (Weber  et 
al  .  ,  1983).   Rao  and  Davidson  (1980)  concluded  that  many 
sorption  reactions  are  complete  within  one  minute  in  batch 
slurry  experiments,  although  longer  times  to  equilibrium 
were  noted  in  several  studies  (Karickhoff  et  al .  ,  1979, 
Miller  and  Weber,  1984).   Schwarzenbach  and  Westall  (1981), 
in  a  comparison  of  solute  breakthrough  at  various  flow 

velocities,  concluded  that  K   values  from  column  experiments 

_3 
where  velocity  was  less  than  10    cm/second  were  similar  to 

K   values  from  18  hour  equilibrium  batch  studies. 

3.4.6  Aromatic  Sorption  Values  From  The  Literature 
There  are  few  data  in  the  literature  addressing 

sorption  of  dissolved  gasoline  components  in  the  subsurface. 
Much  of  the  research  involved  aromatic  compounds  in  single 
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solute  experiments,  simple  mixtures,  or  data  from  crude  oil 
studies . 

Houzim  (1978)  observed  a  decrease  in  sorption  in  the 
order  alkenes  >  aromatics  >  cycloalkanes  >  alkanes. 
Nathawani  and  Phillips  (1977)  in  a  study  of  hexadecane,  o- 
xylene,  toluene  and  benzene  in  crude  oil  on  soils  of  varying 
organic  matter  presented  sorption  coefficients  based  on 
Freundlich  isotherms.   Rodgers  et  al .  (1980)  reported  the 
adsorption  and  desorption  of  benzene  on  several  soils  and 
clays  at  25  C.   The  aqueous  phase  concentration  range  was  10 
to  1000  ug/L.   Sorption  of  benzene  was  minimal,  except  on 
aluminum  saturated  clay.   These  data  are  summarized  in  Table 
3-2. 

Wilson  et  al .  (1981)  evaluated  the  sorption  of  toluene 
on  a  fine  sand  in  a  column  study.   A  retardation  factor  less 
than  2  for  the  concentration  range  of  200-900  ug/L  was 
reported.   This  indicates  the  relatively  low  retardation 
potential  of  sandy  aquifers.   The  retardation  factor 
describes  the  extent  of  solute  transport  relative  to  water. 
The  retardation  factor  for  water  is  defined  as  unity. 
Solutes  with  large  retardation  factors  are  less  mobile  and 
and  their  movement  is  retarded,  relative  to  that  of  water. 

Schwarzenbach  and  Westall  (1981)  presented  data  for  the 
sorption  of  several  chlorinated  and  alkyl  benzenes  on  twelve 
natural  aquifer  materials  with  varying  amounts  of  organic 
carbon.   The  initial  concentrations  of  the  alkylbenzene 
components  were  20  ug/L.   Sorption  coefficients  from  batch 
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studies  of  a  soil  with  low  organic  carbon  (0.0015  g   /g 
soil)  are  shown  in  Table  3-3. 

As  may  be  noted  from  this  short  review/  most  of  the 
above  studies  involve  data  from  individual  components  or 
from  oil  based  products.   Given  the  differences  in 
composition  among  these  petroleum  products  and  gasoline/ 
extrapolation  may  be  insufficient  to  provide  accurate  data 
(Brookman  et  al . ,  1985). 

3. 5  Biodegradation  of  Aromatic  Hydrocarbons  in  Groundwater 


Biological  activity  is  an  important  process  in  the 

attenuation  of  gasoline  hydrocarbons  in  the  subsurface 

environment.   This  realization  is  only  recent.   Early 

techniques  for  the  enumeration  of  microbes  in  the  subsurface 

(Waksman,  1916)  underrepresented  the  numbers  of  microbes  in 

the  subsurface,  showing  a  decline  in  population  with  depth. 
These  data  resulted  from  the  use  of  nutrient  rich  growth 

media,  inappropriate  for  the  enumeration  of  groundwater 

bacteria  (Wilson  and  McNabb ,  1983). 

Recent  work  shows  that  more  substantial  populations  of 

heterotrophic  organisms  exist  in  shallow  water  table 

aquifers  than  were  previously  thought.   Wilson  et  al . 

(1983a)  demonstrated  that  the  numbers  of  organisms  were 

relatively  constant  to  a  depth  of  six  meters  in  a  shallow 

water  table  aquifer.   The  populations  of  heterotrophic 

bacteria  were  estimated  to  be  approximately  10 

organisms/gram  dry  weight  soil  (Ghiorse  and  Balkwill,  1985). 
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Table  3-3.   Sorption  coefficients  of  selected  aromatic 
hydrocarbons  on  low  organic  carbon  soil . 


Kd 


Compound  average    standard  deviation 


Toluene                        0.3  7  0.12 

p-Xylene                       0.50  0.10 

1,3,5-Trimethylbenzene         1.00  0.16 

1,2,3-Trimethylbenzene         0.95  0.11 

Source:   Schwarzenbach  and  Westall,  1981. 
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A  review  of  the  techniques  for  the  enumeration  and 

estimation  of  microbial  biomass  were  presented  by  Atlas 
(1982)  and  Webster  et  al .  (1985).   Bouwer  and  McCarty  (1884) 
noted  that  the  majority  of  bacterial  activity  was  associated 
with  bacteria  attached  to  surfaces.   This  results  in  the 
formation  of  biofilms,  which  are  favored  in  low  substrate- 
high  surface  area  conditions.   The  biofilm   may  also  present 
an  active  surface  with  solutes  sorbing  to  the  surfaces  of 
microbial  cells.   In  terms  of  the  advection-dispersion 
models,  the  rates  of  biological  degradation  are  incorporated 

into  the  model  through  the  sink  term  Q. /  which  describes  the 

l 

microbial  degradation  of  solutes  from  the  aqueous  phase. 
Q.  is  defined  as: 

Q.  -   -k  9  C  [3.18] 

where  k  is  the  rate  of  biological  degradation  ( T~  ) ,  9  is 

3 
the  volumetric  water  content  (ml/cm  )and  C  is  the  solution 

phase  concentration  of  a  solute  (ug/L). 

3.5.1   Environmental  Factors  Affecting  Biodegradation 

Many  factors  can  affect  the  transformation  of  organic 

contaminants  in  the  subsurface.   McCarty  (1984)  included  low 

substrate  concentrations,  toxic  conditions,  molecular 

structure  of  the  substrate,  inaccessibility  of  the 

substrate,  and  absence  of  essential  arowth  facfTs 

Biological  activity  is  often  limited  bv  certain  mpfaHnH" 

requirements  of  the  cell,  supplied  from  the  environment. 

Important  geochemical  properties  include  dH  •  redox 

potential,  nitrogen  and  phosphorus  concentrations  and  the 
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availability  of  an  appropriate  electron  acceptor-   Oxygen  is 

used  as  the  ultimate  electron  acceptor  for  aerobic 

degradation  processes  and  is  often  a  limiting  factor  in  the 

degradation  of  hydrocarbons.   Molecular  oxygen  is  also 

essential  to  the  aerobic  metabolism  of  aromatic  compounds, 

because  it  is  incorporated  into  the  structure  of  the 

metabolic  products  (Evans,  1977).   The  biochemistry  of  the 

aerobic  metabolism  of  aromatic  compounds  is  well  established 

(Dagley,  1975).    The  first  step  in  this  metabolic  pathway 

is  the  removal  of  side  chains,  followed  by  the  enzyme 

(oxygenases)  mediated  hydroxylation  of  the  aromatic  ring. 

Assuming  50%  conversion  of  carbon  to  biomass  and  incomplete 

oxidation  of  the  hydrocarbon  molecules,  two  parts  of  oxygen 

are  required  for  the  degradation  of  each  part  hydrocarbon 

(Wilson  et  al . ,  1986).   The  complete  oxidation  of 

hydrocarbon  molecules  to  CO_  and  H„0  may  require  three  to 
four  parts  of  oxygen  per  part  hydrocarbon. 

There  is  some  evidence  for  the  anaerobic  biodegradation 

of  aromatic  compounds  in  the  environment.   In  the  absence  of 

oxygen,  nitrates,  sulfates  and  CO   become  electron 

acceptors.   Bouwer  and  McCarty  (1984)  presented  a  review  of 

these  processes.   Nitrate  respiration  ( Psuedomonas  and 

Moraxella  sp . )  and  methanogenic  fermentation  processes  can 

reduce  the  benzene  nucleus  followed  by  hydrolysis  to  yield 

aliphatic  acids  (Evans,  1977).   Wilson  and  Rees  (1985) 

showed  the  anaerobic  degradation  of  benzene,  toluene, 

xylenes  and  alkylbenzenes  under  methanogenic  conditions. 
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Over  a  six  week  period  only  toluene  showed  substantial 
degradation,  but  after  40  weeks,  benzene  was  reduced  by  72%, 
toluene  by  99%,  ethylbenzene  by  74%  and  o-xylene  by  78%. 
Nutrient  addition  decreased  the  rate  of  hydrocarbon  removal. 
The  metabolic  products  from  the  anaerobic  degradation  of  the 
aromatic  molecules  were  not  investigated.   Nitrate 
respiration  of  xylene  in  a  river  alluvium  was  demonstrated 
by  Kuhn  et  al .  (1985).   However,  anaerobic 

biotransformations  occur  extremely  slowly  (months),  relative 
to  aerobic  processes  which  may  be  completed  in  a  matter  of 
hours  (Wilson,  1985). 

Physical  properties  of  the  aquifer  also  play  an 
important  role  in  determining  the  extent  of  microbial 
degradation.   Porosity  and  hydraulic  conductivity  are 
significant  parameters  since  the  resupply  of  oxygen, 
substrate  and  nutrients  to  the  microbial  cells  must  come  via 
the  groundwater. 

The  concentration  of  the  contaminant  substrates  is  an 
important  factor  in  the  extent  of  biodegradation .   High 
concentrations  may  result  in  incomplete  degradation 
resulting  from  rapid  depletion  of  oxygen  and  high  substrate 
concentrations  may  also  lead  to  increased  acclimation  times. 
Jensen  et  al.  (1985)  demonstrated  an  increase  in  the  time 
required  for  acclimation  (lag  time)  of  bacterial  cultures 
with  increasing  concentrations  of  napthelene.   Lag  times 
prior  to  substantial  microbial  degradation  of  a  solute  or 
nutrient  reflects  the  time  required  by  the  indigenous 
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microflora  to  adapt  to  the  added  substance.   Adaption  is  a 
phenomenon  rather  than  a  mechanism  or  process,  and  the  term 
refers  to  an  increase  in  the  rate  of  biotransformation  of  a 
substance  resulting  from  exposure  to  that  substance  (Wilson 
et  al .  ,  1983b).   Low  solute  concentrations  may  result  in  the 
occurrence  of  a  threshold  limit,  below  which  the  microflora 
are  unable  to  utilize  the  solute  with  out  a  cosolute  (Wilson 
and  McNabb,  1983).   Jensen  et  al .  (1985)  demonstrated  the 
degradation  of  aromatic  molecules  to  less  than  1  ug/L , 
implying  that  there  was  a  very  low  threshold  limit  for  the 
aromatic  hydrocarbons.   The  relationship  between 
concentration  and  biodegradation  was  reviewed  by  Alexander 
(1985).   He  stressed  the  importance  of  studying  contaminant 
levels  that  exist  in  the  environment. 
3.5.2   Aromatic  Biodegradation  Values  From  the  Literature 

McKee  et  al.  (1972)  reported  the  oxidation  of  gasoline 
by  Pseudomonas  and  Arthrobacter  under  aerobic  but  not 
anaerobic  conditions.   Degradation  of  gasoline   by 

Pseudomonas  was  reported  by  Williams  and  Wilder  (1971),  and 
Litchfield  and  Clark  (1973)  showed  significant  numbers  (10 
cells/mL)  of  hydrocarbon  degrading  bacteria  in  groundwater 
contaminated  with  petroleum  hydrocarbons  from  twelve  sites. 
Bacterial  populations  appeared  to  be  related  to  the 
concentrations  of  hydrocarbons.   These  data  indicate  the 
adaptation  of  microbial  communities  to  the  changing  nutrient 
source  (i.e.,  gasoline).   The  two  major  mechanisms  of 
adaptation  are  induction  of  metabolic  pathways,  or  the 
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activation  or  transfer  of  plasmids  (Litchfield,  1986).   This 
ability  of  microorganisms  to  adapt  to  the  presence  of 
contaminants  forms  the  basis  of  in-situ  biodegradation . 

Several  researchers  have  reported  the  biodegradation  of 
aromatic  compounds  in  groundwater.   One  shortcoming  of  most 
of  this  research  is  the  lack  of  degradation  rate  coefficient 
data,  required  for  use  in  groundwater  transport  models,  and 
insufficient  data  on  solute  concentrations. 

Jamison  et  al .  (1976)  reported  the  use  of  benzene  as  a 

sole  carbon  source.   No  rate  coefficient  data  were  given. 

HcKenna  and  Heath  (1976)  noted  the  slow  oxidation  of  benzene 

ky  P»  putida.   Delfino  and  Miles  (1985)  showed  the 

degradation  of  benzene  in  16  days  under  aerobic  conditions 

in  Floridan  groundwater  with  an  eight  day  lag  phase. 

Ethylbenzene  was  degraded  as  a  sole  carbon  source  (Gibson 

and  Yeh  ,  1973),  but  no  rate  data  were  given.   Schwarzenbach 
et  al .  (1983)  found  toluene  rapidly  degrades  within  several 

meters  in  a  study  of  river  water  infiltration  to  groundwater 

but  rate  and  initial  concentration  data  were  not  specified. 

Kappeler  and  Wuhrmann  (1978b)  in  a  study  of  gas  oil 

degradation  reported  that  nitrogen  and  oxygen  were  the  . 

limiting  factors  in  hydrocarbon  degradation.   Addition  of 

NH  CI  resulted  in  further  microbial  degradation,  and  cell 

densities  were  on  the  order  of  10  /mL.   Lag  times  of  5-6 

days  were  noted  in  the  batch  experiments.   Kappeler  and 

Wuhrmann  (1978a)  showed  that  microbes  from  uncontaminated 

groundwater  can  attack  gas  oil  components.   Lag  times  of  1 
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day  at  25  C  and  5  days  at  10  C  were  reported  in  column 
studies  with  mixed  autochthonous  flora  from  clean 
groundwater.   The  meta  and  para  isomers  of  xylene  and  1,2,4- 
trimethylbenzene  were  degraded  more  rapidly  than  o-xylene , 
1,  2,3-t.rimethylbenzene  or  1 ,  3  ,  5-trimethylbenzene  .   These 
studies  by  Kappeler  and  Wuhrmann  (1978a, b)  make  up  the  bulk 
of  the  work  on  degradation  of  alkyl  substituted  benzenes. 
3.5.3   In-situ  Biodegradation 

Application  of  in-situ  bioremediation  technology  for 
the  renovation  of  hydrocarbon  contaminated  aquifers  is  based 
primarily  on  the  work  of  Raymond  et  al .  (1975a, b)  and 
Raymond  et  al .  (1977)  at  Suntech.   Nutrients  and  oxygen  were 
introduced  with  injection  wells,  and  circulated  through  the 
aquifer  with  pumping  wells.   This  technique  and  other  bio- 
remediation methods  were  reviewed  by  Wilson  et  al .  (1986). 
These  authors  noted  that  studies  are  needed  to  investigate 
the  effectiveness  of  natural  biorestoration  and  to  evaluate 
whether  enhancement  of  natural  processes  is  possible  or 
desirable . 

Transport  of  sufficient  oxygen  to  subsurface  microbes 
is  a  major  technical  problem.   Oxygen  is  only  slightly 
soluble  in  water  and  is  quickly  depleted  during  aerobic 
biodegradation.   Oxygen  addition  by  air  sparging,  oxygen 
sparging  and  the  use  of  hydrogen  peroxide  are  documented  in 
the  literature  (Lee  and  Ward,  1984).   The  use  of  hydrogen 
peroxide  appears  particularly  advantageous  (TRI,  1982). 
Hydrogen  peroxide  is  relatively  inexpensive,  nonpersistent 


2  2 

[2°2 


OH*  +   HO    >   HO   +   H+   +   0  *        [3.21] 


Britton  (1985)  reported  that  hydrogen  peroxide  was 
relatively  stable  in  combination  with  phosphates,  even  in 
the  presence  of  moderate  iron  concentrations/  and  that 
bacterial  populations  can  tolerate  EUO,  concentrations  up  to 

500  mg/L„   Hydrogen  peroxide  was  shown  (Britton,  1985)  to 

2 
increase  microbial  counts  by  10  ,  but  there  was  no  reported 

increase  in  hydrocarbon  removal. 

3.5.4   Measurement  of  Microbial  Activity 

The  reduction  of  INT  ( 2-p-iodophenol-3-p-nitrophenyl-5- 

phenyl  tetrazolium  chloride)  to  INT-formazan  by  the  electron 

transport  system  is  a  function  of  cell  respiration,  and  is 

widely  used  as  a  general  measure  of  microbial  activity. 
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and  is  more  soluble  in  water  than  air  or  molecular  oxygen. 
However,  it  is  also  cytotoxic  and  may  be  chemically  reduced, 
especially  in  the  presence  of  iron  salts.   The  biological 
decomposition  of  hydrogen  peroxide  is  enzymatic: 

catalase 

2H2°2   >   2H2°   +   °2  [3.18] 

peroxidase 
H202   +   XH2   >   2H20  +   X         [3.19] 

where  X  is  a  biological  reducing  agent.   Non-enzymatic 
decomposition  occurs  most  frequently  in  the  presence  of  iron 
salts : 

Fe  '   +  HO    >   Fe+++  +  OH-  +  OH*         [3.20] 
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This  technique  is  recommended  as  an  index  of  general 
microbial  activity  of  soil  microorganisms  (Klein  et  al., 
1971). 

Reduction  of  INT  to  INT-formazan  is  a  sensitive  assay 
for  dehydrogenase  activity.   The  INT-formazan  is  easily 
extracted  from  sediments  and  soils  by  methanol,  and  the  INT- 
formazan  complex  is  stable.   Trevors  et  al.  (1982)  found  a 
high  correlation  between  electron  transport  system  activity 
and  oxygen  consumption.   Klein  et  al .  (1971)  presented  a 
rapid  and  simple   procedure  for  the  determination  of 
dehydrogenase  activity  using  INT  in  soils  with  low  organic 
carbon . 

3.6  Summary 


This  literature  review  has  presented  some  of  the  basic 
principles  required  as  a  basis  for  the  discussion  of  the 
experimental  work  reported  in  this  dissertation,  and  has 
highlighted  some  of  the  important  findings  relative  to  the 
dispersion,  sorption  and  biodegradation  of  aromatic 
compounds  in  groundwater  systems. 


CHAPTER  IV 
MATERIALS  AND  METHODS 


4.1   Introduction 


This  chapter  discusses  the  materials  and  experimental 
methods  employed  during  this  study.   The  field  site/  and  the 
solutes  and  sorbents  are  described  followed  by  a  description 
of  the  chromatographic  systems.   Laboratory  experiments  for 
the  determination  of  hydrolysis/  sorption  and  biodegradation 
parameters  are  discussed.  Finally/  the  field  procedures  and 
experiments  are  discussed. 

4. 2  Site  Description 

The  field  research  site  used  for  a  portion  of  this 
study  was  located  at  the  Citrus  Research  and  Education 
Center  (CREC)  at  Lake  Alfred,  Fl  .   The  site  was  located  in 
the  Trail-Ridge  Lake  Wales  Ridge  system  of  hills  containing 
deep  internally  drained  lake  basins.   Unconsolidated 
deposits  in  the  area  consisted  of  sand  and  sandy  clays  up  to 
150  ft  thick  above  the  limestone  bedrock.   The  geology  was 
marked  by  many  sinkholes  formed  through  subsidence  of  the 
unconsolidated  deposits  into  solution  cavities  in  the 
limestone  (Spangler,  1984). 
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The  research  site  was  located  on  the  rim  of  an  ancient 
sinkhole.   The  surficial  aquifer  was  composed  of  sand  and 
clayey  sands.   An  continuous  clayey  confining  layer  of 
uneven  depth  was  present  between  7  to  12  ft  below  land 
surface.   This  layer  supported  a  saturated  zone  between  3  to 
6  ft  in  thickness.   Local  relief  was  from  156  ft  (above  mean 
sea  level)  at  the  top  of  the  hill  at  the  eastern  boundary  of 
the  site,  to  131  ft  in  the  wetland  area  at  the  west  edge  of 
the  site.   A  site  map  is  shown  in  Figure  4-1.   The  surficial 
aquifer  was  comprised  of  medium  angular  grained  sands  and 
fill  material.   The  hydrology  of  the  site  was  discussed  by 
Killan  (1987). 

The  surficial  aquifer  was  contaminated  during  the 
spring  of  1983  by  the  loss  of  7500-8000  gallons  of  leaded 
gasoline  from  a  storage  tank.   Free  floating  gasoline  was 
removed  by  surface  skimming  as  of  May  1985.   The  outline  of 
the  contaminated  area  as  of  October,  1986  is  shown  in  Figure 
4-2.   The  plume  was  defined  by  determination  of  explosive 
gas  concentrations  in  bore  holes  throughout  the  site.  These 
data  were  confirmed  by  GC  analysis  of  soil  cores  and  the  use 
of  ground  penetrating  radar.  These  techniques  were  described 
in  detail  by  Killan  (1987). 

4. 3  Aquifer  Material 


Aquifer  materials  used  in  this  research  were  obtained 
from  the  field  research  site  at  the  IFAS-Citrus  Research  and 
Education  Center  at  Lake  Alfred  Florida.   A  site  map  is 
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Figure  4-2.   Extent  of  the  hydrocarbon  plume  at  the 
field  research  site  as  of  October,  1986 


38 

shown  in  Figure  4-1.   All  experiments  were  carried  out  with 
subsamples  of  the  same  aquifer  material.   The  sample  was 
collected  with  a  stainless  steel  auger  just  below  the  water 
table  at  a  depth  of  about  4  feet,  approximately  10  feet  east 
of  Well  RAP-1.   Care  was  taken  to  avoid  contamination  with 
surface  materials  by  removing  one  foot  of  top  soil,  and 
through  careful  handling  of  the  auger.   The  aquifer  material 
was  oven  dried  at  105  C  for  24  hours,  sieved  through  2  mm 
standard  sieve  and  stored,  covered,  at  room  temperature. 
Prior  to  use,  the  aquifer  material  was  autoclaved  for  90 
minutes  on  each  of  three  consecutive  days  to  sterilize  the 
materials . 

Prior  to  sterilization  and  drying,  selected  physical 
and  chemical  properties  of  the  aquifer  materials  used  in  the 
laboratory  studies  were  characterized  (pH,  particle  density, 
particle  size  analysis,  percent  organic  carbon,  bulk 
density,  hydraulic  conductivity,  and  water  content)  using 
standard  methods  of  soil  analysis  (Black,  1965). 

4.4  Choice  of  Solutes 


Gasoline  contaminated  well  water  from  the  Lake  Alfred 
research  site  was  used  as  the  source  of  dissolved  solutes 
for  the  majority  of  experiments  in  this  study.   With  the 
exception  of  a  single  solute  column  sorption  experiment  with 
benzene  spiked  in  to  RAP-2  water,  all  experiments  were 
performed  with  mixtures  of  dissolved  hydrocarbons  at 
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concentrations  occurring  in  the  field.   These  concentrations 
are  the  result  of  the  solubilization  and  subsequent 
weathering  of  gasoline  hydrocarbons  into  groundwater.   Well 
OHM-4  was  used  as  the  source  of  water  for  these  experiments. 
This  well  was  chosen  based  on  consistently  high  levels  of 
dissolved  aromatic  hydrocarbons.   Hydrocarbon  free 
groundwater  was  obtained  from  a  non-contaminated  portion  of 
the  aquifer  (Well  RAP-2).   Hydrocarbon  concentrations  in 
these  waters  were  monitored  monthly.   Well  RAP-2  remained 
free  of  aromatic  hydrocarbons  throughout  the  course  of  these 
experiments.   Concentrations  of  aromatic  hydrocarbons  varied 
in  Well  OHM-4  but  remained  high  enough  to  provide  samples 
for  laboratory  experiments. 

Water  samples  were  collected  with  a  5.1  cm  (2")  id  poly 
vinyl  chloride  (PVC)  bailer7  following  removal  of  five  well 
volumes  to  allow  collection  of  a  representative  sample. 
Well  volumes  were  calculated  based  on  the  diameter  of  the 
well/  and  the  depth  of  water  in  the  well.   These  water 
samples  were  collected  in  four  liter  brown  glass  bottles 
transported  on  ice/  and  stored  at  4  C  upon  arrival  at  the 
laboratory.   The  pH  of  these  well  waters  ranged  from  6  to  7. 
The  conductivity  was  approximately  300  umhos .   Total 
phosphate  was  0.4  mg/L  for  RAP-2  and  was  0.65  mg/L  for  well 
OHM-4.   Nitrate  was  0.29  mg/L  in  well  RAP-2  and  0.20  in  well 
OHM-4. 
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The  single  solute  column  experiment  with  benzene  used 
RAP-2  water  spiked  with  benzene  (Aldrich,  gold  label  99.9%) 
to  yield  a  solution  of  4000  ug/L  benzene. 

4. 5   Hydrocarbon  Analyses 

Gas  chromatographic  analyses  of  hydrocarbons  for  this 
study  were  performed  on  two  systems.   These  are  described 
below. 
4.5.1   GC/MS  Analyses 

Field  samples  collected  before  September,  1986,  and  the 
initial  hydrolysis  vials  were  analyzed  for  volatile  aromatic 
constituents  using  a  Hewlett  Packard  model  5985B  GC/MS/COMP 
system  equipped  with  a  10  port  Tekmar  Automatic  Liquid 
Sampler  (ALS)  and  Liquid  Sample  Concentrator  (LSC)   purge 
and  trap  system.   EPA  method  624  was  used.   Separation  of 
analytes  was  achieved  with  a  0.32  mm  i.d.,  30  meter  long, 
DB-5  (1  um  film  thickness)  fused  silica  capillary  column 
(J  &   W  Scientific),  with  manual  liquid  nitrogen 
cryofocusing.   The  1,4  isomer  of  dichlorobenzene  was  used  as 
an  internal  standard.   Response  factors  for  benzene, 
toluene,  ethylbenzene  and  o-xylene  were  determined  relative 
to  the  internal  standard  and  used  for  quantitation. 
Response  factors  for  meta  and  para  xylene  were  assumed  to  be 
the  same  as  the  ortho  isomer.   A  response  factor  of  1  was 
assumed  for  the  CgH12  hydrocarbons.   Chromatographic 
conditions  were  as  follows: 


41 


mass  range  45-450  amu 

Temp  1  30  C 

Temp  2  280  C 

Rate  5  C/min 

Hold  time  10  minutes 
Cryofocus  time    5  minutes 

Pre-cool  2  minutes 


4.5.2   GC  Anal 


yses 


Hydrocarbon  analyses  were  performed  on  a  Perkin  Elmer 
model  8410  gas  chromatograph  with  a  flame  ionization 
detector  and  microprocessor  data  system.   Samples  were 
concentrated  by  purge  and  trap  with  a  Tekmar  LSC/ALS  system 
employing  a  modified  version  of  EPA  method  602.   Analytical 
separation  was  achieved  with  a  0.53  mm  i.d.,  30  meter  long, 
fused  silica  Megabore  DB-1  (100%  methylpolysiloxane )  column 
(J  &  W  Scientific)  with  a  3  urn  film  thickness. 

Benzene,  toluene,  ethylbenzene ,  o-xylene  and  m,p-xylene 
were  quantified  using  the  internal  standard  method  (1,4 
dichlorobenzene)  during  August  and  September  1986  for 
monthly  analysis  of  field  samples  and  for  day  =  0  hydrolysis 
ampules.   After  this  date,  eight  isomers  of  C  H    were 
identified  and  confirmed  by  analysis  of  individual  standards 
and  were  quantified  in  all  chromatograms  along  with  BTEX 
(benzene  +  toluene  +  ethylbenzene  +  m, p, o-xylene )  compounds. 
The  internal  standard  was  changed  from  1 , 4-dichlorobenzene 
to  chlorobenzene  to  avoid  co-elution  problems.   A  complete 
description  of  this  analytical  method  and  a  summary  of  the 
quality  control  parameters  for  this  method  are  in  Appendix 
A. 
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The  raeta  and  para  isomers  of  xylene  were  not  resolved 
on  either  chromatographic  system  employed  in  this  research. 
The  combination  of  these  analytes  was  reported  as  m,p- 
xylene.   Likewise,  3-ethyl toluene  and  4-ethyltoluene  were 
not  resolved  with  the  analytical  system  employed  in  this 
study,  and  the  combined  concentrations  of  these  analytes 
were  reported  in  this  study  with  the  abbreviation  3,4- 
ethyltoluene . 

4.6   Hydrolysis  Studies 


Hydrolysis  studies  were  performed  in  5  mL  glass  ampules 
(Fisher  Scientific).   Ampules  were  rinsed  with  methanol  and 
oven  dried  at  105  C  .   Ten  microliters  of  hydrocarbon 
contaminated  groundwater  were  spiked  into  ampules  containing 
5  mL  of  buffer  solution.   Buffer  solutions  were  prepared 
with  non  contaminated  well  water,  and  the  pH  was  adjusted  to 
pH  ■  2.0,  7.0,  9.2,  and  12.0  with  0.01  M  phosphate  buffers. 
The  ampules  were  sealed  with  an  ampule  sealer  ( Oceanographic 
International,  College  Station,  TX) ,  and  autoclaved  (1  hour 
at  120  C).   One  set  of  ampules  was  analyzed  at  time  zero. 
Another  set  of  ampules  was  stored  in  the  dark  at  20,  40  and 
60  C  and  analyzed  by  gas  chromatography  (GC)  after  60  days. 
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4.7   Batch  Sorption  Studie: 


Sorption  batch  studies  were  performed  in  40  mL  VOA 
vials  with  Teflon  coated  septa  (Fisher  Scientific).   Vials 
were  first  filled  with  60  g  of  aquifer  material,  and  then 
autoclaved  at  120  C  for  1  hour  on  each  of  three  consecutive 
days  . 

Water  from  well  OHM-4,  containing   a  mixture  of 
dissolved  aromatic  hydrocarbons,  was  used  in  the  batch 
sorption  experiments.   As  a  result,  all  these  experiments 
are  multisolute,  at  concentrations  representative  of  those 
found  across  the  aquifer  at  Lake  Alfred. 
4.7.1   Sorption  Experiments 

Water  used  in  the  sorption  experiments  was  filter 
sterilized  through  0.2  um  membrane  filters  (Gelman  Metricel) 
and  then  added  to  each  vial.   The  range  of  solute 
concentrations  was  achieved  by  dilution  of  Well  OHM-4  water 
with  Well  RAP-2  water  at  ratios  between  1:1  to  1:1000.   Each 
dilution  was  performed  in  triplicate.   Non-soil  controls 
(solute  water  with  no  soil)  were  also   set  up  in  triplicate. 
To  minimize  headspace,  the  vials  were  premixed  on  a  rotary 
tumbler  for  approximately  1  hour  to  remove  interstitial  air 
and  to  disperse  the  foam  that  formed  during  mixing.   The 
vials  were  then  opened,  completely  filled  with  sample  and 
recapped.   A  high  solids  to  solution  ratio  (2.9  g/  g  )  was 
used  to  maximize  the  fractional  decrease  in  solution 
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concentration  owing  to  sorption,  and  to  more  closely 
simulate  natural  aquifer  conditions. 

Vials  used  in  sorption  experiments  were  equilibrated  at 
room  temperature  (20  +  2  C)  on  a  rotary  tumbler  at 
approximately  20  rpm  for  18  hours,  and  then  centrifuged  at 
800  G  for  30  minutes.   Samples  were  analyzed  by  purge  and 
trap/gas  chromatography.   Vials  used  in  the  batch  sorption 
kinetic  rate  study  were  sampled  at  1,  2,  4,  8,  16,  24,  36 
and  48  hours. 
4.7.2   Desorption  experiments 

Desorption  experiments  were  conducted  subsequent  to  a 
sorption  experiment.   Following  centrifugation  and  sampling 
for  sorption  losses,  approximately  10  mL  of  supernatant  were 
removed  and  replaced  with  hydrocarbon  free  water  (Well  RAP- 
2).   The  vials  were  re-equilibrated  for  24  hours  on  the 
rotary  tumbler,  centrifuged  and  sampled.   Each  vial  was  only 
desorbed  one  time.   These  experiments  were  not  designed  to 
calculate  desorption  isotherms  or  test  isotherm 
nonsingularity. 
4-7.3   Calculation  of  Sorption  Coefficients 

The  amount  of  solute  sorbed  to  the  aquifer  material  ( ng 
solute/gram  soil)  was  calculated  by  determining  the 
difference  between  the  solution  concentration  of  the  non- 
soil  blanks  and  the  soil  containing  vials.   The  amount  of 
solute  lost  was  divided  by  the  solution  to  soil  ratio  to 
normalize  the  data  to  a  ng/gram  basis.   Sorption 
coefficients  were  calculated  by  fitting  isotherm  data  to 
three  models;  linear,  linear  with  intercept  forced  through 
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zero/  and  the  log  normalized  (Freundlich)  models  (Miller  and 
Weber,  1984). 

4. 8   Column  Sorption  Studies 

4.8.1   Experimental  Procedures 

Leaching  column  experiments  were  performed  with  a  25  x 
250  mm  glass  preparative  chromatography  column  (Altex   cat. 
no.  252-18)  with  a  Teflon  coated  adjustable  plunger  (Nkedi- 
Kizza  et  al . ,  1987).   Aquifer  material  was  dry  packed  into 
the  column  which  was  then  autoclaved  at  120   C  for  1  hour. 
The  solutes  were  pumped  from  2.6  L  Teflon  gas  sampling  bags 
(Alltech  Associates,  Deerfield,  IL)  with  a  Gilson  model  302 
HPLC  pump  fitted  with  a  model  5s  pump  head  (Gilson  Medical 
Electronics,  Middleton-,  WI).   The  flow  range  of  this  system 
was  0.005  -  5.00  mL  per  minute.   All  transfer  lines  and 
connections  were  Teflon  or  stainless  steel  to  minimize 
interaction  of  solutes  with  reactive  surfaces.   Column 
length  was  adjusted  to  5.0  cm.   Flow  rates  through  the 
column  were  set  at  1  ml/min  (0.204  cm/min)  for  sorption 
studies.   Column  effluent  breakthrough  curves  (BTCs)  were 
measured  under  steady  saturated  water  flow  conditions  with 
continuous  application  of  solute  containing  water. 

Effluents  from  the  sorption  columns  were  collected 
manually  in  1  ml  crimp  seal  vials.   These  column  effluents 
were  either  analyzed  immediately  or  stored  at  4  C  in  1  mL 
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crimp  seal  vials  with  Teflon  coated  septa  for  later 
analysis.   All  samples  were  analyzed  within  48  hours. 

The  breakthrough  of  an  unretained  solute  was  determined 
for  each  column  using  calcium  chloride  (1  ml/min  columns). 
Breakthrough  curves  were  determined  by  spiking  hydrocarbon 
free  groundwater  from  Lake  Alfred  (RAP-2)  with  chloride  (600 
mg/L  CaCL2).   Chloride  analyses  were  performed  with  a 
chloridometer  automatic  titrator  (Buchler-Cotlove ) . 
Chloride  ion  was  not  expected  to  be  adsorbed  owing  to  the 
low  cation  exchange  capacity  of  the  Lake  Alfred  soil. 

Well  water  used  in  the  sorption  experiments  was 
filtered  through  0.2  urn  membrane  filters  (Gelman  Metricel ) 
directly  into  the  Teflon  bags.   The  bags  were  autoclaved 
prior  to  each  use.   Columns  were  saturated  with  filter 
sterilized  water  from  well  RAP-2  prior  to  the  input  of 
solute  containing  water. 
4-8.2   Estimation  of  Retardation  Factor  (R)  in  Columns 

Three  methods  were  used  to  estimate  the  value  of  R  from 
the  column  data.   In  method  1,  retardation  factors  (R  )  were 
calculated  by  fitting  the  solution  of  Brenner  (1962)  to  the 
column  effluent  curves.   Peclet  numbers  used  for  these 
calculations  were  determined  from  the  breakthrough  of  the 
non-retained  solutes  having  retardation  factors  equal  to 
unity.   Method  2  was  based  on  the  conservation  of  mass 
principle.  This  method  calculated  retardation  factors  (R  ) 


by  evaluating  the  area  above  the  breakthrough  curve  usi 


ng 


Simpson  s  Rule  (Swokowski,  1975).   The  R   value  was  assumed 
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equal  to  the  area  above  the  BTC  when  the  effluent 

concentration  (C)  divided  by  the  influent  concentration  ( CQ ) 

was  plotted  vs  pore  volume  as  described  by  equation  [4.1] 

pv 

^  max 

r   =      f   n-C/C_]  dpv  [4.1] 


/   [1-C/C  ]  dpv 

0 

where  pv     is  the  total  number  of  pore  volumes  displaced 
r  max 

through  the  column,  and  pv  is  pore  volumes  (Nkedi-Kizza  et 
al./  1987).   This  method  assumed  a  mass  balance  existed  in 
the  soil  columns.  The  third  method  set  the  retardation 
factor  (R   )  to  equal  the  number  of  pore  volumes  required 
for  the  effluent  concentration  of  each  analyte  to  reach  0.5 
of  the  influent  concentration.   The  use  of  this  method 
assumes  that  the  breakthrough  curve  is  symmetrical  and 
sigmoidal,  and  that  equilibrium  conditions  exist  between  the 
solution  and  sorbed  concentrations  during  leaching  through 
the  column  (Nkedi-Kizza  et  al . ,  1987).   The  value  of  Kd  was 
calculated  from  the  various  R  values  with  equation  [3.10]. 

4, 9    Hydrogen  Peroxide  Evaluation 


The  reaction  rate  of  hydrogen  peroxide  in  the  aquifer 
environment  was  simulated  by  monitoring  the  dissolved  oxygen 
(DO)  (YSI  model  5739  probe  and  YSI  model  54A  DO  meter), 
redox  potential  (platinum  redox  electrode,  Fisher 
Scientific)  and  pH  (gel  membrane  electrode,  Fisher 
Scientific)  of  well  water  and  aquifer  material  in  a  3  arm 
500  mL  reaction  flask.   Contaminated  well  water  was 
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equilibrated  at  room  temperature  (20  +  2  C)  in  the  sealed 
flask.   Hydrogen  peroxide  (50%)  was  added  undiluted  in 
microliter  quantities  and  at  various  dilutions.   Aquifer 
material  was  then  added  to  assess  the  ability  of  the 
material  to  catalyze  the  reaction.   The  50%  hydrogen 
peroxide  stock  was  titrated  with  0.01N  potassium 
permanganate  (Dupont,  1984)  to  check  the  strength  of  the 
stock  solution.   The  standardized  stock  was  then  used  to 
make  the  appropriate  dilutions  without  further  calibration. 

4. 10   Batch  Biodegradation  Studies 


4.10.1   Experimental  Procedure 

Batch  biodegradation  experiments  were  performed  in  40 
mL  VOA  vials  as  described  for  the  batch  sorption 
experiments . 

Well  water  from  OHM-4  was  used  as  the  source  of  both 
dissolved  aromatic  hydrocarbons  and  bacteria  in  these 
studies.   The  water  was  not  filtered  prior  to  use.   This 
experiment  was  designed  to  evaluate  the  ability  of  adapted 
groundwater  bacteria  to  degrade  mixtures  of  dissolved 
aromatic  solutes  at  field  scale  concentrations.   The 
experimental  design  for  batch  biodegradation  experiment 
number  1  is  shown  in  Table  4-1.   Seven  treatments  were  set 
up,  with  15  replicate  vials  per  treatment.   Water  from  Well 
OHM-4  was  added  (350  mL)  to  a  500  mL  erlenmeyer  flask,  and 
then  amended  with  hydrogen  peroxide  (50%),  ammonium  chloride 


Table  4-1.   Experimental  design  for  batch  biodegradation 
experiment  #1. 
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Treatment 

Hydrogen 

Peroxide 

(mg/L) 

NH  CI 
(mg/L) 

1A 

none 

none 

IB 

17 

none 

1C 

68 

none 

ID 

none 

18 

IE 

17 

18 

IF 

68 

18 

1G 

none 

none 

Sodium 

Azide 

(mg/L) 


none 
none 
none 
none 
none 
none 
1.25 
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(Reagent  grade,  Fisher  Scientific)  and  10%  (w/v)  aqueous 
solution  of  sodium  azide  (Fisher  Scientific)  as  outlined  in 
Table  4-1.   Triplicate  samples  were  analyzed  for  each 
treatment  at  0/  3/  7,  15,  and  31  days.   Treatment  number  1G 
was  a  sterile  control.   Hydrogen  peroxide  was  added  based  on 
data  from  the  hydrogen  peroxide  evaluation  experiment  and  on 
data  from  Britton  (1985)/  who  demonstrated  that  cytotoxicity 
was  minimal  at  hydrogen  peroxide  concentrations  less  than 
100  mg/L .   Ammonium  chloride  was  added  based  on  data  from 
Mitchell  (1974)  who  found  that  ammonia  nitrogen  is 
assimilated  quickly  during  microbial  growth.   Ammonia  (as 
NH  CL)  was  added  to  achieve  quantities  calculated  to  meet 
nitrogen  requirements  of  the  bacteria. 

Biodegradation  experiment  number  2  was  designed  to 
evaluate  the  efficacy  of  oxygen  gas  in  addition  to  hydrogen 
peroxide  (Table  4-2).   Sterile  controls  were  maintained  in 
treatments  2C/  2F  and  21.   Ammonium  chloride  and  hydrogen 
peroxide  were  added  as  in  biodegradation  experiment  #1. 
Oxygen  was  added  by  bubbling  0_  gas  into  a  closed  3  liter 
flask  filled  with  1200  mL  of  contaminated  well  water.   A 
valve  allowed  for  pressure  relief.   Water  was  released 
through  a  glass  tube  at  the  bottom  of  the  flask/  fitted  with 
a  teflon  stopcock.   Vials  were  filled  as  described  in  the 
sorption  experiments.   The  vials  used  in  both  batch 
biodegradation  experiments  were  placed  in  an  incubator  (20  ± 
1  C)  and  inverted  once  every  two  days  to  provide  mixing. 
Samples  were  taken  at  0/  2,  7,  14,  21  and  35  days. 


Table  4-2.  Experimental  design  for  batch  biodegradation 

experiment  #2. 
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Treatment  # 

Oxyg 
addi 

en 
tion 

NH4C1 
(mg/L) 

Sodium 
Azide 
(mg/L) 

2A 

air 

none 

none 

2B 

air 

18 

none 

2C 

air 

none 

1.25 

2D 

60 

mg/L 

H2°2 

none 

none 

2E 

60 

mg/L 

H2°2 

18 

none 

2F 

60 

mg/L 

H2°2 

none 

1.25 

2G 

°2 

saturation 

none 

none 

2H 

°2 

saturation 

18 

none 

21 

°2 

saturation 

none 

1.25 
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Following  sample  removal  for  GC  analysis,  dissolved  oxygen 
was  measured  in  each  batch  biodegradation  vial  (batch 
experiments  1  and  2)  with  a   YSI  model  5739  DO  probe  and  YSI 
model  54A  DO  meter  (Yellow  Springs  Instruments  Co.). 

Microbial  activity  was  assessed  through  the  measurement 
of  INT  reduction  to  INT-formazan  (Klein  et  al . ,     1971).   Ten 
grams  of  soil  from  each  vial  were  placed  in  sterile  50  mL 
Erlenmeyer  flasks.   Each  flask  was  amended  with  1  mL 
distilled  water  and  1.5  mL  of  0.4%  (w/v)  aqueous  solution  of 
filter  sterilized  (0.2  urn  Gelman  Metricel  membrane  filters) 
INT  (Eastman-Kodak  Co.).   The  soil  was  mixed  with  a  sterile 
glass  rod,  capped  with  aluminum  foil  and  incubated  at  20  C 
for  72  hours.    Sterile  controls  were  prepared  by 
autoclaving  several  flasks  for  3  consecutive  days  for  90 
minutes.   Approximately  3  grams  (dry  weight)  of  soil  were 
removed  from  each  flask  following  incubation  and  placed  in  a 
test  tube.   Ten  mL  of  methanol  were  added  to  each  tube  and 
the  contents  were  mixed  on  a  vortex  mixer  for  1  minute,  then 
centrifuged  at  800G  for  20  minutes.   The  INT-  formazan  in 
the  methanolic  extract  was  measured  spectrophotometrically 
at  480  nm  against  a  methanol  extract  of  soil  containing  no 
INT.   The  INT-formazan  concentration  was  derived  from  a 
standard  curve  of  INT-formazan  in  methanol. 


4-10.2   Calculation  of  Biological  Rate  Constants 

Aqueous  concentration  data  from  the  batch 
biodegradation  vials  were  used  with  the  regression  equations 
from  the  Freundlich  fit  of  the  batch  desorption  data  to 
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calculate  the  amount  of  solute  lost  to  sorption  in  each 
batch  vial.   The  extent  of  the  sorption  loss  correction 
varied  with  the  concentration  of  the  analyte  and  the 
sorption  parameters  Kfd  and  n.   This  correction  factor  added 
as  much  as  20%  to  the  measured  concentration  values.   These 
predicted  losses  resulting  from  sorption  were  added  to  the 
aqueous  concentration  for  each  analyte  in  each  vial  to 
calculate  the  total  concentration  of  each  solute  in  the  vial 
(Ct).   These  C   values  were  employed  to  obviate  the  need  for 
simultaneous  calculation  of  biodegradation  on  both  sorbed 
and  aqueous  concentrations  and  any  calculation  of  rates  of 
sorption-desorption  during  biodegradation.   The  C   values 
were  used  to  model  the  biological  rate  coefficients.   The 
rate  data  were  fitted  to  zero  order,  first  order,  second 
order  and  mixed  order  rate  equations  (Levenspiel,  1972),  and 
to  the  Thomas  slope  method  (Thomas,  1950). 

4. 11   Column  Biodegradation  Studies 


4.11.1   Experimental  Procedure 

Column  biological  degradation  experiments  were 
performed  with  the  same  column  system  described  in  section 
4.8.  Two  flow  rates  were  used  in  these  experiments.  These 
were  1  ml/min  (0.680  cm/min)  and  0.90  ml/hr  (0.010  cm/min). 
Columns  operated  at  0.90  ml/hr  were  fitted  with  a  low  dead 
volume  in-line  septa.  Effluent  was  withdrawn  with  a  50  ul 
syringe  (Hamilton  Co.)  and  analyzed  immediately.   Effluents 
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from  the  1  mL/min  columns  were  collected  and  analyzed  as 
described  in  section  4.8.   Breakthrough  curves  for  a  non 

retained  solute  were  obtained  with  tritiated  water  (0.01  uci 

3 
H  for  0.90  ml/hr  columns)  or  CaCl2  (1  ml/min  columns). 

3 
Analyses  of   H20  effluents  were  performed  on  a  Delta  300 

model  6890  Liquid  Scintillation  Counter  (Cearle  Analytical) 

qsing  Scintiverse  II  scintillation  cocktail  (Fisher 

Scientific).   Columns  operated  for  biodegradation 

experiments  were  set  up  in  the  same  manner  as  the  sorption 

columns.   Solute  containing  water  was  filtered  through  0.45 

urn  membrane  filters  (Gelman  Metricel)  to  remove 

particulates.   A  standing  microbial  population  was  developed 

in  the  columns  operated  at  1  ml/min  by  inoculation  with 

water  from  well  OHM-4. 

4.11.2   Calculation  of  Rate  Constants 

Rate  constants  for  the  biological  degradation  of 

aromatic  solutes  from  column  breakthrough  curves  were 

determined  by  application  of  the  first  order  rate  equation 

to  the  breakthrough  curve  data  at  steady  state.   Microbial 

degradation  processes  are  often  assumed  to  be  first  order 

(Bossert  and  Bartha,  1984).   Substitution  of  equation  [3.18] 

into  equation  [3.11]  incorporates  the  degradation  term  into 

the  one  dimensional  mass  transport  equation: 


R  30/  3t   -  D  32C/8x2  -  v  3C/  3x  -  k  C      [4.2] 
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Dividing  all  terms  of  equation  [4.2]  through  by  R  and 
defining  D^  =  D/R,  v*  =  v/R  and  k*  -  k/R,  equation  [4.2] 
becomes 


3  C/  3  t   =  Dh*  32C/  9X2  -  V*  gc/9x  -  k*  C      [4.3] 


at  steady  state  conditions  where  3C/3t  =  0,  equation  4.3 
reduces  to: 


Dh    32C/  3  x2  -  v  3C/  3  x  +  k  C   =  0 


[4.4] 


For  steady  state  conditions  there  is  no  sorption  effect 
since  the  R  term  cancels  out  and  the  rate  of  biodegradation 
(k)  may  then  be  calculated  by  application  of  the  first  order 
rate  equation  to  the  portion  of  the  BTC  which  is  at  steady 
state.  The  first  order  rate  equation  for  this  system  is: 


C/Co   =   e   kt  [4.5] 

For  a  column  of  length  x,  and  with  a  pore  water  velocity  of 
v,    time  may  be  expressed  as 

t  =  x/v  [4.6] 

and  substitution  of  [4.6]  into  [4.5]  and  rearrangement 
allows  calculation  of  the  rate  constant  for  biodegradation: 


-  (  In  C/C   )  v/x 

o 


[4.7] 
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Average  C/C   values  were  calculated  from  the  regions  of  th 


o 


solute  breakthrough  curves  where  3C/ St  =  0.   This  derivation 
assumes  that  microbial  degradation  occurs  only  from  the 
aqueous  phase,  and  that  dispersion  is  negligible. 

4.12   Field  Studies 

4.12.1   Aquifer  Characterization 

A  tracer  experiment  was  conducted  to  measure  seepage 
velocities  and  obtain  better  estimates  of  aquifer  hydraulic 
conductivity  and  field  scale  dispersion.   RAP-9  was  used  as 
the  dosing  well.   The  following  steps  outline  the 
experimental  procedure: 

1.  A  tracer  solution  was  prepared  by  dissolving 
50  lb  (23  kg)  of  technical  grade  ammonium  chloride  in  55  gal 
(208  L)  of  tap  water.    The  resulting  concentration  was 
109,000  mg/L  ammonium  chloride. 

2.  The  ammonium  chloride  solution  was  injected 
into  the  dosing  well  and  simultaneously  diluted  with  tap 
water  at  a  metered  rate  of  1  gallon  per  minute  (gpm). 
Dosing  continued  for  15.8  hours,  resulting  in  a  total  dose 
volume  of  1,035  gallons  of  tracer  solution. 

3.  Detection  of  the  tracer  was  monitored  in 
wells  RAP-10  and  RAP-11  using  a  conductivity  meter  with  a 
field  probe.    Measurements  were  obtained  at  one  half-  to 
one-hour  intervals  for  the  first  24-hour  period.    Wells  P- 
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6,  P-7,  UF-1E,  RAP-4,  OHM-4,  UF-2M  and  UF-3W  were  also 
monitored  periodically  for  the  following  two  weeks. 

The  breakthrough  of  the  tracer  was  calculated  in  pore 
volumes  using  the  equation 

pv  =  vt/L  [4.8] 

where  pv  is  pore  volumes,  t  is  time  (hours),  L  is  the 
distance  between  RAP-9  and  RAP-10  (5  ft)  and  v  is  the 
seepage  velocity  ( ft/hour ) from  work  by  Killan  (1987). 
4.12.2   Water  Quality  Monitoring 

Samples  for  hydrocarbon  analysis  were  taken  from 
selected  monitoring  wells  monthly  from  February  1,  1986  to 
June  1987.   Sampling  procedures  are  detailed  in  Appendix  B. 
Hydrocarbon  analyses  were  as  described  previously  using  the 
three  chromatographic  methods  as  they  became  available.   The 
pH ,  temperature,  dissolved  oxygen,  and  conductivity  were 
measured  periodically  in  selected  wells.   Sampling 
procedures  are  detailed  in  Appendix  B. 

Total  phosphorus  concentrations  were  determined  in  all 
monitoring  wells  (EPA  method  365.1).   All  phosphorus  forms 
were  converted  to  orthophosphate  by  autoclaving  with 
potassium  permanganate  in  an  acidic  medium.   Orthophosphate 
was  determined  spectrophotometrically  at  880  nm  with  a 
Perkin  Elmer  model  552  spectrophotometer. 

Chloride  was  determined  colorometr ically  in  each 
monitoring  well  over  several  months  to  determine  background 
levels  of  chloride  via  EPA  method  325.1.   Average  background 
concentrations  were  20+8  mg/L .   Nitrate  was  measured  with 
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an  Orion  nitrate  electrode  via  standard  method  418b  (APHA, 

1980). 

4.12.3   Microbial  analyses 

Viable  microbial  cells  were  enumerated  by  plate  count 
technique  using  dilute  soil  extract  agar  (DSEA)  media.   This 
technique  was  developed  based  on  work  by  Ghiorse  and 
Balkwill  (1983)  and  Wilson  et  al .   (1983). 

DSEA  was  prepared  by  autoclaving  100  g  of  surface  soil 
in  100  mL  of  distilled  water  for  one  hour  at  120  C.   The 
supernatant  was  filtered  (Whatman  glass  fiber  filters)  to 
remove  particulates  and  diluted  ten  fold  with  distilled 
water  and  amended  1.5%  (w/v)  with  agar  (Fisher  Scientific). 

Ten  grams  of  subsurface  material  were  suspended 
aseptically  in  100  mL  0.1%  sodium  pyrophosphate  (Fisher 
Scientific)  then  appropriate  dilutions  were  plated  in 
triplicate  on  DSEA  media.   All  plates  were  incubated 
aerobically  at  27-30  C  for  ten  days. 


CHAPTER  V 
RESULTS  AND  DISCUSSION 


5.1   Introduction 


This  chapter  will  review  the  results  of  all  experiments 
performed  as  part  of  this  dissertation  research.   The 
presentation  of  results  and  the  interpretation  of  the  data 
for  each  major  experimental  section  are  grouped  together  to 
avoid  loss  of  continuity.   Hydrolysis  of  aromatics  in 
groundwater  is  discussed  first,  followed  by  the  results  and 
discussion  of  the  batch  and  column  sorption  experiments. 
Batch  and  column  biodegradation  experiments  are  addressed 
next>  followed  by  the  presentation  of  field  data,  and  the 
correlation  of  field  data  with  laboratory  experiments. 

5. 2   Hydrolysis  of  Aromatic  Hydrocarbons 

The  results  of  initial  measurements  (time  =  0)  of  the 
hydrolysis  ampules  were  inconclusive,  resulting  from  over- 
dilution  of  the  samples.   Several  analytes  were  below  the 
limit  of  detection  of  the  GCMS  analytical  system,  therefore 
no  conclusive  statements  may  be  made  relative  to  the  rates 
of  aromatic  hydrolysis. 

Statistical  analysis  of  data  from  ampules  after  60  days 
of  temperature  controlled  storage  indicated  that  for  a  given 
temperature,  there  was  no  significant  (student's  t-test, 
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0.05  level)  change  in  concentration  of  analytes  over  the  pH 
range  tested  (pH  2/7/9,12). 

Only  one  compound,  1 , 2 , 3-trimethylbenzene ,  demonstrated 
a  significant  (student's  t-test ,  0.01  level)  temperature 
effect  at  pH  values  of  7,9  and  12.   This  change  in 
concentration  occurs  only  at  40  C,  and  the  concentration 
values  at  20  and  60  C  are  equivalent  for  all  pH  values. 
These  data  are  not  consistent  with  data  from  other  aromatic 
compounds  in  this  study,  which  showed  no  change  in 
concentration  with  varying  temperatures.   The  apparent  loss 
of  solute  seen  for  1 , 2 , 3-trimethylbenzene  was  likely  the 
result  of  working  near  the  detection  limit  of  the  analytical 
system,  or  the  result  of  experimental  error. 

The  absence  of  concentration  differences  across  a  wide 
range  of  pH  and  temperature  for  the  60  day  ampules  implies 

that  hydrolysis  was  not  a  significant  mechanism  for  removal 
of  aromatic  hydrocarbons.   This  was  expected,  owing  to  the 
resistance  of  aromatic  structures  to  nucleophilic  attack  by 
water.   This  results  from  the  electronegativity  associated 
with  the  delocalization  of  electrons  in  the  pi  bonds  of  the 
aromatic  nucleous.   McCarty  (1984)  has  noted  that  chemical 
hydrolysis  may  occur,  but  that  for  most  compounds  this 
process  was  slow  relative  to  biological  removal  rates.   In 
addition,  hydrolysis  results  in  simple  changes  in  the 
molecular  structure,  whereas  biological  transformations 
often  result  in  the  mineralization  of  organic  compounds  to 
carbon  dioxide  and  water. 
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5.3   Characterization  of  Aquifer  Materials 

An  analysis  of  a  sub-sample  of  the  aquifer  materials 
used  in  the  laboratory  experiments  is  presented  in  Table  5- 
1.   All  experiments  with  aquifer  material  were  performed 
with  subsamples  of  well-mixed  aquifer  material.   Single  size 
fractions  of  aquifer  material  were  not  used  since 
extrapolation  from  one  size  fraction  to  another  has  been 
shown  to  lead  to  errors  in  the  estimation  of  sorption  values 
(Abdul  et  al.,  1986).   Unwashed,  natural  sorbent  material 
from  the  Lake  Alfred  site  was  used  in  this  study  to  more 
closely  approximate  field  conditions.   The  organic  carbon 
content  of  this  material  was  low,  and  particle  size  analysis 
indicated  the  dominance  of  fine  to  medium  grained  sands. 
The  pH  of  the  aquifer  material  was  in  the  range  suitable  for 
biological  degradation,  and  was  consistent  with  the  pH 
values  in  the  well  water. 

5.4   Batch  Sorption  Studies 


5.4.1   Sorption  Rate  Studies 

Rate  studies  were  conducted  to  determine  the  sorption 
kinetics  of  the  selected  aromatic  solutes  with  Lake  Alfred 
aquifer  material.   These  experiments  established  the 
equilibration  time  for  the  sorption  isotherms.   The  approach 
to  equilibrium  is  shown  in  Figure  5-1.   These  curves  are 
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Table  5-1.   Selected  physical  and  chemical  Properties 
of  the  Lake  Alfred  aquifer  material. 


Parameter  Value 

PH  7.4  (0.01M  CaCl  ) 

Particle  Density  2.6  g/mL 

Water  Content  (by  weight)       24% 
Organic  Carbon  0.015% 

Bulk  Density  1.4  g/mL 

Particle  Size  Analysis 

clay  1.8% 

silt  1.7% 

very  fine  sand  3.0% 

fine  sand  38.2% 

medium  sand  47.0% 

coarse  sand  8.2% 

very  coarse  sand  0.2% 
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marked  by  an  initial  rapid  sorption,  and  equilibrium 
conditions  are  established  within  several  (4  to  8)  hours. 
These  data  are  in  agreement  with  Weber  et  al.  (1983)  who 
stated  that  sorption  reactions  with  natural  sorbents  were 
generally  rapid  and  not  rate  limited.   Based  on  these  data 
an  equilibration  time  of  18  hours  was  chosen.   Eighteen 
hours  was  chosen  to  maximize  the  time  for  sorption  yet 
minimize  the  time  for  losses  from  the  system  (ie./  via 
diffusion  of  solutes  through  the  Teflon  septum).   This  was 
equivalent  to  time  scales  used  in  previous  studies  (Chiou  et 
al.,  1979,  1983;  Schwarzenbach  and  Westall,  1981).   Longer 
equilibration  times  were  not  possible  using  this 
experimental  technique  since  losses  in  non-soil  blanks  after 
3  days  made  it  difficult  to  differentiate  between  sorption 
and  loss  from  the  system.   The  equilibration  time  used  in 
this  study  did  not  guarantee  that  the  sorption  process  was 
complete,  but  that  it  was  complete  to  the  extent  that  it 
could  be  accurately  measured. 

5.4.2   Batch  Sorption  Isotherm  Data 

Data  for  the  equilibrium  batch  sorption  isotherms  are 
presented  in  Appendix  C.   Solution  concentrations  are  in 
ug/L ,  and  sorbed  concentrations  are  in  ng/g.   Three  models 
were  fitted  to  these  data  using  the  method  of  least  squares 
regression  analysis.   These  models  were  the  linear,  linear 
with  suppressed  fit  (forced  through  the  origin),  and  the 
Freundlich  (log-log  transformed).   The  results  of  these 
analyses  for  the  linear  models  are  presented  in  Table  5-2 
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Table  5-2. 


Regression  parameters  for  the  analysis  of 
average  values  of  equilibrium  batch  isotherm 
sorption  data  with  the  linear  model. 


Compound 

Na 

(ug?£) 

Benzene 

14 

950 

Toluene 

14 

4200 

m,p-Xylene 

16 

4300 

o-Xylene 

16 

2500 

3  or  4  ETe 

11 

935 

l/3,5-TMBf 

11 

460 

2-ET9 

11 

373 

l/2,4-TMBh 

10 

1600 

1/2,3-TMB1 

10 

558 

K 


d' 


std 


y-mt  , 


std 


0.066,  0.005 
0.049,  0.010 
0.095,  0.005 
0.097,  0.004 
0.087,  0.012 
0.142,  0.008 
0.106,  0.011 
0.131,  0.006 
0.124,     0.010 


1.5,  4.28  0.914 
10.9,  65.6  0.694 
3.4,  28.7  0.961 
3.1,  10.8  0.979 
3.8,     12.6  0.861 

1.6,  3.92  0.973 
0.81,    4.66  0.910 

4.6,  10.6  0.981 

1.7,  6.10  0.951 


number  of  data  points 


maximum  concentration 


'standard  deviation 


y-intercept 
'3  or  4  Ethyltoluene 

1,3, 5-Trimethylbenzene 
'2-Ethyltoluene 

1 ,  2,4-Trimethylbenzene 


1,2,  3-Trimethylbenzene 
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and  Table  5-3.   The  Freundlich  regression  parameters  are 
presented  in  Table  5-4. 

There  is  no  significant  difference  (student's  t-test , 
0.05  level)  between  sorption  coefficients  predicted  with  the 
linear  model  and  those  predicted  with  the  linear  model  with 
suppressed  intercept.   The  sorption  coefficients  are 
determined  form  the  slope  of  the  linear  isotherms.   In 
addition,  statistical  determination  of  the  confidence 
intervals  of  the  y-intercepts  indicate  that  there  is  no 
significant  difference  between  the  predicted  value  of  the  y- 
intercept  in  the  linear  model  and  zero  at  the  0.05 
significance  level,  confirming  that  these  two  models  are 
analogous.   Equivalence  between  these  two  models  is  expected 
since  the  sorbed  concentration  should  equal  zero  when  no 
solute  is  added  to  the  system.   A  non  zero  intercept  is  an 
indication  of  nonlinearity  in  the  isotherm.   In  most  cases 
both  linear  models  fit  the  data  well  as  evidenced  by  the 
relatively  high  coefficients  of  determination  (r2).   Based 
on  these  analyses,  the  isotherms  for  the  sorption  of 
aromatic  solutes  from  Lake  Alfred  water  onto  Lake  Alfred 
aquifer  material  were  concluded  to  be  linear.   The 
coefficients  of  determination  for  toluene  in  both  linear 
models  were  substantially  lower  than  for  other  compounds  in 
this  study.   The  linear  model  with  suppressed  intercept 
accounted  for  only  63.6%  of  the  total  sum  of  squares 
deviations  about  the  means  for  the  14  values  in  the  toluene 
isotherm.   This  suggested  that  this  model  was  not 
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Table  5-3 


Regression  parameters  for  the  analysis  of 
average  values  of  equilibrium  batch  isotherm 
sorption  data  with  the  linear  model 
(suppressed  intercept). 


Compound 

Na 

c       b 

«it:= 

stdC 

y-mt 

2 

r 

=  =  =  =  =  : 

sssscssssssis::: 

======== 

=========== 

======= 

Benzene 

14 

950 

0.069, 

0.005 

0 

0.904 

Toluene 

14 

4200 

0.051, 

0.009 

0 

0.636 

m,p-Xylene 

16 

4300 

0.096, 

0.004 

0 

0.960 

o-Xylene 

16 

2500 

0.099, 

0.003 

0 

0.978 

3    or    4    ETe 

11 

935 

0.093, 

0.010 

0 

0.850 

l,3,5-TMBf 

11 

460 

0.146, 

0.007 

0 

0.969 

2-ET9 

11 

373 

0.108, 

0.009 

0 

0.908 

l,2,4-TMBh 

10 

1600 

0.135, 

0.005 

0 

0.978 

1/2,3-TMB1 

9 

558 

0.128, 

0.008 

0 

0.948 

number  of  data  points 


maximum  concentration 


'standard  deviation 


y-intercept 
'3  or  4  Ethyltoluene 
'1,3, 5-Trimethylbenzene 
'2-Ethyltoluene 

1, 2,4-Trimethylbenzene 


1,2,3-Trimethylbenzene 
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Table  5-4. 


Regression  parameters  for  the  analysis  of 
average  values  of  equilibrium  batch  isotherm 
data  with  the  Freundlich  model. 


Compound 

IS" 

============ 

=  =  =  =  = 

Benzene 

14 

Toluene 

14 

m,p-Xylene 

16 

o-Xylene 

16 

3  or  4  ETf 

11 

1,3,5-TMB9 

11 

2-ETh 

11 

1,2,4-TMB1 

10 

1,2,3-TMB11 

9 

C  =v      log    log 

<B|7u   k  *   stdc 


std 


950  -0.857,  0.217  0.901,  0.071  0.952 

4200  -0.783,  0.260  0.876,  0.076  0.917 

4300  -0.758,  0.110  0.916,  0.031  0.984 

2500  -0.873,  0.098  0.966,  0.026  0.990 

935  -0.702,  0.158  0.904,  0.054  0.958 

460  -0.605,  0.081  0.921,  0.029  0.991 

373  -0.951,  0.255  0.993,  0.088  0.934 

1600  -0.641,  0.099  0.937,  0.035  0.989 

558  -0.623,  0.078  0.918,  0.028  0.995 


number  of  data  points 


maximum  concentration 


'log  standard  deviation  of  Kf  values 

Freundlich  exponent 

standard  deviation  of  Freundlich  exponent 
"3  or  4  Ethyltoluene 
'1,3, 5-Trimethylbenzene 

2-Ethyltoluene 

1 , 2,4-Trimethylbenzene 


1,2,3-Trimethylbenzene 
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appropriate  for  estimation  of  sorption.   However,  analysis 
of  variance  with  a  global  F-test  indicated  that  the  model 
was  useful  for  predicting  sorption  at  the  0.01  significance 
level.  Therefore,  the  linearity  of  all  isotherms  was 
confirmed.   Linear  isotherms  have  been  noted  by  several 
authors  ( Schwar zenbach  and  Westall ,  1981,  Chiou  et  al . , 
1979,  Karickhoff  et  al . ,  1979)  and  the  data  presented  in 
this  study  are  in  agreement  with  these  studies. 

Curtis  et  al .  (1986)  noted  that  the  use  of  the  linear 
regression  technique  was  not  statistically  rigorous  since 
variance  in  the  dependent  variable  was  not  distributed 
uniformly  across  the  observed  concentration  range.   These 
authors  suggested  that  a  least  squares  fit  on  the  log 
transformed  data  (Freundlich  model)  gives  a  better 
approximation  by  providing  a  more  uniform  distribution  of 
variance.   The  Freundlich  model  provided  a  good  fit  to  the 
data  in  this  study  (Table  5-4)  as  evidenced  by  the  high 
coefficients  of  determination  for  the  Freundlich  model.   The 
Freundlich  isotherm  explained  between  91.7  to  99.5%  of  the 
variance  in  the  data,  and  provided  a  slightly  improved  fit 
to  the  isotherm  data  relative  to  the  linear  models.   The  r2 
value  for  toluene  was  0.917,  which  was  much  improved  over 
the  coefficient  of  determination  for  toluene  in  the  linear 
model.   The  log  Kf  values  for  the  study  compounds  are  also 
presented  in  Table  5-4.   Values  for  several  components 
(ethylbenzene ,  and  the  propylbenzenes )  were  not  included  in 
the  table  owing  to  their  low  concentrations  in  Well  OHM-4 
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water  on  the  day  samples  were  collected  for  this  study.   The 
linearity  of  these  isotherms  was  confirmed  by  the  values  of 
the  regression  coefficients  and  the  Freundlich  exponents 
(n),  both  of  which  were  close  to  unity.   As  n  approaches 
unity  the  models  should  converge  since  the  linear  model  is 
in  effect  a  special  case  of  the  Freundlich  model. 

A  comparison  of  the  Freundlich  and  linear  models  is 
presented  in  Table  5-5.   The  deviation  between  predicted 
amounts  of  sorption  for  the  linear  model  with  suppressed 
intercept  and  Freundlich  models  are  expressed  as  ratios 
between  the  calculated  sorbed  concentrations.   This  method 
was  used  to  evaluate  the  predictive  equivalence  of  both 
models  over  the  concentration  ranges  encountered  in  this 
study.   This  method  was  chosen,  since  direct  comparison  of 
K   and  K   values  may  be  misleading  owing  to  the  log 

transformation  of  the  data  in  the  Freundlich  isotherm  model. 
The  largest  deviations  occurred  for  toluene  at  1  ug/L.   In 
general  the  ratios  approached  unity  as  the  concentrations 
increased,  but  diverged  in  the  range  between  1  to  50  ug/L. 
This  comparison  indicates  that  the  models,  were  essentially 
similar,  as  is  predicted  from  the  values  of  the  Freundlich 
exponent  (n).   As  n  approaches  unity,  the  Freundlich 
isotherm  approaches  the  linear  isotherm.   The  convergence  of 
these  models  is  confirmed  by  an  examination  of  2- 
ethyltoluene  in  Table  5-5.   This  compound  has  the  highest 
Freundlich  constant  (n  =  0.993)  and  the  ratios  of  predicted 
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Table  5-5.   Ratio  of  sorbed  concentrations  calculated 

from  Freundlich  and  linear  equilibrium  models 


concen 

trations 

(ug/L) 

Solute 

1 

50 

100 

500 

1000 

Benzene 

2.01a 

1.73 

1.27 

1.09 

1.01 

Toluene 

3.23 

1.97 

1.81 

1.48 

1.36 

ni/p-Xylene 

1.81 

1.30 

1.23 

1.08 

1.01 

o-Xylene 

1.35 

1.18 

1.16 

1.09 

1.07 

3  or  4  ETb 

2.14 

1.48 

1.38 

1.18 

1.10 

1,3,5-TMBC 

1.70 

1.25 

1.18 

1.04 

0.99 

2-ET 

1.04 

1.01 

1.00 

0.99 

0.98 

1,2,4-TMBS 

1.69 

1.34 

1.28 

1.15 

1.10 

l,2,3-TMBf 

1.88 

1.35 

1.27 

1.12 

1.05 

athe  ratio  of  the  amount  sorbed  as  calculated  from  the 
Freundlich  model  to  the  amount  sorbed  predicted  from  the 
linear  model  with  suppressed  intercept,  |t  the  same 
solution  concentration. 

3  or  4  Ethyltoluene 

cl,3, 5-Trimethylbenzene 


2-Ethyltoluene 
!1 , 2,4-Trimethylbenzene 
1/2, 3-Trimethylbenzene 
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sorption  from  the  two  models  are  consistently  close  to  one 
over  the  entire  concentration  range  tested. 

Freundlich  isotherms  for  benzene  and  toluene  are  shown 
in  Figures  5-2  and  5-3.  Graphs  of  Freundlich  isotherms  for 
the  remaining  solutes  are  presented  in  Appendix  C.   Average 
values  are  plotted  in  Figures  5-2  and  5-3  and  error  bars 
showing  one  standard  deviation  in  the  experimental 
determination  of  the  sorbed  concentrations  are  presented  to 
give  an  indication  of  the  variance  in  these  data.   Standard 
deviations  for  all  compounds  are  shown  in  Appendix  C.   The 
influence  of  dissolved  organic  carbon  was  not  assessed 
during  this  study.  However,  based  on  the  work  of  Curtis  et 
al.,  (1986)  with  a  sandy  aquifer  material  (0.02%  organic 
carbon),  organic  carbon  in  this  study  was  not  expected  to 
decrease  the  values  of  K,  by  more  than  5%.   Water  from  the 
Lake  Alfred  aquifer  was  used  in  these  experiments,  and  the 
organic  carbon  in  solution  was  assumed  to  be  in  equilibrium 
with  the  organic  carbon  on  the  aquifer  material.   Therefore, 
dissolution  of  additional  organic  carbon  into  solution 
should  have  been  minimal,  and  K   should  not  be  greatly 
affected.   This  hypothesis  was  confirmed  by  evaluation  of 
the  partitioning  model  as  a  predictive  technique  for 
sorption  of  aromatic  solutes  to  the  Lake  Alfred  aquifer 
material  in  section  5.6.   The  interaction  between  organic 
carbon  and  the  solutes  was  shown  to  be  low. 
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5.4.3   Batch  Desorption  Experiments 

Desorption  data  are  also  presented  in  Figures  5-2  and 

5-3,  fit  with  the  Freundlich  type  model.   Visual  inspection 

of  the  desorption  data  suggests  some  degree  of 

irreversibility  or  some  difference  in  desorption  kinetics 

based  on  the  upward  displacement  of  the  desorption 

regression  lines.   However,  the  calculated  values  of  the 

partition  coefficient  for  desorption  with  the  linear  type 

model  (Table  5-6)  and  for  the  linear  type  model  with 

suppressed  intercept  (Table  5-7)  were  not  significantly 

different  from  sorption  values  (Kd)  at  the  0.05  probability 

level.   Desorption  coefficients  from  the  Freundlich  type 

model  (Kjrj)  were  also  not  significantly  different  from  Kf 
f  d 

values  at  the  0.05  level  (Table  5-8).   Statistical  analyses 
of  the  models  used  to  evaluate  the  desorption  coefficients 
indicated  that  all  three  models  gave  excellent  fit  to  the 
data,  as  evidenced  by  the  high  coefficients  of 
determination. 

These  data  suggested  the  reversibility  of  the  sorption 
process,  and  demonstrated  that  the  hysteretical  behavior  of 
the  desorption  data  were  not  significant.   This  was 
consistent  with  a  majority  of  the  published  literature  on 
sorption  of  organic  compounds  to  natural  sorbents  (Miller 
and  Weber,  1984) . 

For  purposes  of  discussion,  the  Kd  values  from  the 
linear  model  with  suppressed  intercept  are  used  in  the 
following  sections.   As  discussed  earlier,  these  data  were 


Table  5-6.   Regression  parameters  for  the  analysis  of 
average  values  of  equilibrium  batch 
desorption  data  with  the  linear  model. 
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Compound 


N 


(ug7£)    K 


dd' 


std 


y-int    std 


Benzene 

5 

950 

Toluene 

8 

4200 

m,p-Xylene 

9 

4300 

o-Xylene 

9 

2500 

3  or  4  ETe 

9 

935 

l,3,5-TMBf 

8 

460 

2-ET9 

7 

373 

l,2/4-TMBh 

10 

1600 

1,2,3-TMB1 

9 

558 

0.248/  0.006 

0.303,  0.011 

0.186,  0.006 

0.152,  0.012 

0.250,  0.018 

0.194,  0.004 

0.566,  0.010 

0.199,  0.005 

0.219,  0.032 


2.2,  4.18  0.998 
4.13,  55.8  0.992 
13.3,  30.1  0.993 
11.0,  31.1  0.955 

-2.9,     18.4  0.968 

1.5,    1.62  0.998 

0.78,     2.64  0.999 

3.3,  7.62  0.996 
1.1,     20.3  0.871 


number  of  data  points 

maximum  concentration 
cstandard  deviation 

y-intercept 
e3  or  4  Ethyltoluene 

1,3,5-Trimethylbenzene 
g2-Ethyltoluene 

1, 2,4-Trimethylbenzene 
1 1 , 2 , 3-Trimethylbenzene 
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Table  5-7.    Regression  parameters  for  the  analysis  of 
average  values  of  equilibrium  batch 
desorption  data  with  the  linear  model 
(suppressed  intercept). 


Compound 

Na 

c   b 

Kdd' 

stdC 

y-mt 

2 
r 

Benzene 

5 

950 

0.251, 

0.004 

0 

0.998 

Toluene 

8 

4200 

0.304, 

0.008 

0 

0.992 

m,p-Xylene 

9 

4300 

0.190, 

0.005 

0 

0.992 

o-Xylene 

9 

2500 

0.159, 

0.010 

0 

0.951 

3  or  4  ETe 

9 

935 

0.256, 

0.014 

0 

0.967 

l,3,5-TMBf 

8 

460 

0.1968 

0.003 

0 

0.997 

2-ET9 

7 

373 

0.570, 

0.007 

0 

0.999 

l,2,4-TMBh 

10 

1600 

0.202, 

0.004 

0 

0.996 

1,2,3-TMB1 

9 

558 

0.221, 

0.022 

0 

0.870 

number  of  data  points 

maximum  concentration 
cstandard  deviation 

y-intercept 
e3  or  4  Ethyltoluene 

1,3, 5-Trimethylbenzene 
g2-Ethyltoluene 

1, 2,4-Trimethylbenzene 
"""l ,  2,3-Trimethylbenzene 
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Table  5-8. 


Regression  parameters  for  the  analysis  of 
average  values  of  equilibrium  batch 
desorption  data  with  the  Freundlich  model 


Compound 


<5|7l> 


log  log 
Kf,  std( 


<3    i.  ^e 

n  ,  std 


Benzene 

5 

950 

Toluene 

8 

4200 

m,p-Xylene 

9 

4300 

o-Xylene 

9 

2500 

3  or  4  ETf 

9 

935 

1,3,5-TMB9 

8 

460 

2-ETh 

7 

373 

1,2,4-TMB1 

10 

1600 

1,2,3-TMB3 

9 

558 

-0.635,  0.168  1.02,  0.079  0.982 

-0.478,  0.039  0.992,  0.019  0.998 

-0.561,  0.096  0.963,  0.043  0.986 

-0.159,  0.263  0.774,  0.099  0.891 

-0.692,  0.068  1.029,  0.030  0.994 

-0.569,  0.051  0.956,  0.025  0.996 

-0.657,  0.079  0.983,  0.046  0.989 

-0.658,  0.122  0.997,  0.056  0.981 

-0.624,  0.092  0.988,  0.038  0.990 


number  of  data  points 

b  .  . 

maximum  concentration 

clog  standard  deviation  of  K,  values 

Freundlich  exponent 
standard  deviation  of  Freundlich  exponent 

3  or  4  Ethyltoluene 
^1,3,5-Trimethylbenzene 

2-Ethyltoluene 
"""l^^-Trimethylbenzene 
-^1,2,3-Trimethylbenzene 
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or  linear  models,  and  this  model  is  more  convenient  for  the 

application  of  equation  [3.10]. 

5.5   Breakthrough  Curves  for  Aromatic  Solutes 

5.5.1   Measurement  of  Column  Dispersion 

Breakthrough  curves  (BTCs)  for  a  non-retained  solute 
were  determined  for  each  column  used  in  these  experiments. 
Analysis  of  these  data  allowed  the  determination  of  the 
Peclet  number  used  to  model  the  breakthrough  of  the  aromatic 
solutes.   Evaluation  of  these  data  also  allowed  the 
calculation  of  dispersion  in  the  column.   Chloride  and 
tritiated  water  were  used  in  these  experiments. 

Dispersion  (D  )  was  calculated  from  the  slope  of  a  plot 
of  C/Co  vs  pore  volumes  (pv)  at  pv  =  1  according  to  the 
equation  (Rao,  1985): 


D,  =  [  v  L/  4  pi  B^] 


[5.1] 


where  D   is  the  hydrodynamic  dispersion  coefficient 

2 

(cm  /min),  v  is  the  flow  velocity  (cm/min),  L  is  the  length 

of  the  column  (cm)  and  B  is  the  slope  of  the  BTC  at  C/Co  =1. 
This  assumes  a  sigmoidal  shaped  curve,  and  this  assumption 
was  valid  for  these  breakthrough  curves.   A  typical 
breakthrough  curve  is  shown  in  Figure  5-4.   Some  values  of 
Dh  are  presented  in  Table  5-9.   Columns  with  flow  rates  of  1 
ml/min  exhibited  higher  values  of  D   since  dispersion 
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Table  5-9.    Values  of  dispersion  coefficients  calculated 
from  the  breakthrough  curves  of  unretained 
solutes  in  laboratory  columns. 


Tracer 

Flow 
(mL/min ) 

Velocity 
( cm/min ) 

Dispersion  (D,  ) 
( cm  /min )      avg 

std 

======= 

—  =  =  =  =  =  :=  —  — 

ssaassssss 

========================= 

========= 

V 

0.015 

0.003 

0.00053 

3H2<D 

0.015 

0.003 

0.00066 

«s 

v 

0.015 

0.003 

0.00040 

0.00053 

0.00013 

CaCl2 

1.0 

0.204 

0.051 

CaCl2 

1.0 

0.204 

0.013 

CaCl 

1.0 

0.204 

0.069 

0.044       0.029 


average  values  of  dispersion  measurements 
standard  deviation  of  dispersion  measurements 
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ml/min  exhibited  higher  values  of  D   since  dispersion 

increases  with  increasing  pore  water  velocity  (Roberts  et 

al . ,  1985).   It  may  be  noted  that  the  pore  water  velocity  of 

0.680   cm/min  was  equivalent  to  the  seepage  velocity  in  some 

portions  of  the  aquifer  at  the  Lake  Alfred  field  site. 

These  data  are  compared  to  field  dispersion  data  in  section 

5.10. 

5.5.2   Aromatic  Solute  Breakthrough  Curves 

Breakthrough  curves  for  selected,  dissolved  aromatic 

solutes  in  the  column  effluent  (Well  OHM-4  water)  are  shown 

in  Figure  5-5  (benzene)  Figure  5-6  (toluene)  and  Figure  5-7 

( n-propylbenzene ) .   Breakthrough  curves  for  these  solutes 

are  presented  because  they  show  the  the  breakthrough  of  the 

least  retained  compounds  (benzene  and  toluene)  and  the  most 

retained  (n-propylbenzene).   These  solutes  are  presented 

separately  to  avoid  overlap  on  a  single  plot,  but  are  part 

of  the  multi-component  mixture  resulting  from  the 

solubilization  of  gasoline  into  groundwater  at  the  Lake 

Alfred  site.   Graphical  representations  of  the  remaining 

solutes  in  the  column  effluent  are  shown  in  Appendix  D.   The 

changes  in  effluent  concentration  near  the  end  of  each 

breakthrough  curve  was  consistent  for  each  solute, 

reflecting  the  same  relative  variability.   These  deviations 

may  be  explained  by  heterogeneities  in  flow  paths  in  the 

porous  media,  or  by  analytical  error. 

Calculated  values  of  R,  K,,  and  K    based  on  the 

d        oc 

analyses  of  the  column  data  by  curve  fitting  to  Brenner 
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(1962)  are  shown  in  Table  5-10.   Breakthrough  curve  data  are 

compiled  in  Appendix  D. 

Retardation  factors  for  these  solute  breakthrough 

curves  were  also  evaluated  by  estimating  the  area  above  the 

curve  using  Simpson's  method-   This  method  of  calculation 

yielded  R  values  which  were  slightly  greater  than  the  fitted 

values  (Table  5-11)/  but  which  exhibited  the  same  relative 

elution  order.   The  R  values  calculated  from  the  area  above 

the  BTC  were  within  0.1  of  the  values  calculated  by  curve 

fitting.   Retardation  factors  were  also  calculated  by 

determination  of  the  number  of  pore  volumes  required  to 

reach  C/C   value  of  0.5.   All  three  methods  of  calculation 
o 

are  compared  in  Table  5-11.   The  breakthrough  curves  for  the 
solutes  in  this  study  were  slightly  asymmetrical  whereas  the 
curves  for  the  unretained  solute  were  sigmoidal  and 
symmetrical.   This  asymmetry  was  attributed  to  sorption 
nonequilibrium  during  transport  through  the  column  and  not 
to  dispersion.   The  shape  of  the  measured  BTC  is  determined 
by  the  kinetics  of  the  sorption-desorption  process. 
Symmetrical  BTC  are  obtained  when  the  sorption  process  is 
instantaneous  and  equilibrium  conditions  exist  between 
sorbed  and  aqueous  concentrations.   Under  non-equilibrium 
conditions  during  flow,  asymmetrical  BTC  are  obtained  (Rao 
and  Davidson,  1979).   This  type  of  response  has  been 
reported  by  several  investigators  for  pesticide  breakthrough 
(Nkedi-Kizza  et  al . ,  1987). 
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Table  5-10  . 


Calculated  values  of  R,  K,  and  K    from 
analysis  of  solute  breakthrough  curves. 


Compound 

__inf 

Rb 

„//„„= 

K   d 
_oc_ 

log  K 
____oc 

Benzene 

4700 

1.36 

0.059 

393 

2.60 

Toluene 

2600 

1.55 

0.091 

607 

2.78 

Ethylbenzene 

1800 

1.65 

0.107 

713 

2.85 

iti/p-Xylene 

1700 

1.85 

0.140 

933 

2.97 

o-Xylene 

2200 

1.60 

0.099 

660 

2.82 

Isopropylbenzene 

1000 

2.00 

0.165 

1100 

3.04 

n-Propylbenzene 

560 

2.40 

0.231 

1540 

3.19 

3  or  4  Ethyltoluene 

1600 

2.25 

0.206 

1373 

3.14 

l,3,5-TMBe 

530 

2.15 

0.190 

1267 

3.10 

2-Ethyltoluene 

990 

1.90 

0.148 

987 

2.99 

l,2,4-TMBf 

770 

2.10 

0.181 

1207 

3.08 

1,2,3-TMB9 

1223 

1.84 

0.138 

920 

2.96 

influent  concentration 

calculated  by  curve  fitting  with  Brenner  (1962) 

with  P   =  8. 

e 

'calculated  from  the  relationship  K,  =  (R-l)  g/p    where  0  - 
0.30  and  p  =  1.82 

'k    =  K,/f     where  f    =  0.00015. 
oc     d   oc  oc 

1,3/5-Trimethylbenzene 


1/2/4-Trimethylbenzene 
gl,2,3-Trimethylbenzene 


Table  5-11.   Retardation  factors  calculated  from  leaching 
column  and  equilibrium  batch  isotherm  data 


Compound 

Benzene 

Toluene 

Ethylbenzene 

m,p-Xylene 

o-Xylene 

Isopropylbenzene 

n-Propy lbenzene 

3  or  4  Ethyltoluene 

1/3/ 5-Tri met hy lbenzene 

2-Ethyltoluene 

1/2/4-Trimethylbenzene 

1/2/ 3-Tri me thy lbenzene 


R    a 
a 

R       b 

R.C 
l 

R    d 

1.45 

1.17 

1.60 

1.36 

1.73 

1.33 

1.72 

1.55 

1.73 

1.49 

nde 

1.65 

1.99 

1.66 

1.76 

1.85 

1.68 

1.39 

1.58 

1.60 

2.08 

1.82 

nd 

2.00 

2.49 

2.29 

nd 

2.40 

2.32 

2.05 

1.86 

2.25 

2.26 

2.02 

1.94 

2.15 

2.02 

1.74 

1.49 

1.90 

2.20 

1.92 

1.99 

2.10 

1.93 

1.67 

2.03 

1.84 

R   retardation  factor  calculated  from  the  area  of  the 
breakthrough  curve. 

R    retardation  factor  calculated  from  the  number  of  pore 

volumes  at  C/C   =0.5 
o 

c 
R.  retardation  factor  calculated  from  equilibrium  batch 

isotherm  data. 


R,  retardation  factor  calculated  from  fitting  column  data 
1 

not  determined 


to  the  solution  of  Brenner  (1962). 
e 
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Nkedi-Kizza  et  al.  (1987)   presented  a  method  to  assess 

the  asymmetry  in  a  BTC  by  measuring  the  difference  in  the  R 

values  calculated  by  the  pore  volume  method  (R   )■ ,  from 

those  calculated  from  the  area  above  the  BTC  (R  ).   An 

a 

empirical  index  for  sorption  nonequilibrium  (ISNE)  was 
defined  as : 


ISNE   a   [(R   -  R   )  /  R  ] 
a     pv     aJ 


[5.2] 


where  R   is  the  retardation  factor  calculated  from  the  area 

a. 

above  the  BTC  and  R    is  the  retardation  factor  calculated 

pv 

by  evaluation  of  the  number  of  pore  volumes  required  for  the 

column  effluent  to  equal  0.5  of  the  influent  concentration. 

For  ISNE  equal  to  zero,  R   is  equal  to  R   /  and  for  ISNE 

a  pv 

equal  to  1  ,    R   is  much  less  than  R   .   The  calculated 
a  pv 

values  of  ISNE  for  all  12  compounds  in  this  study  are 
presented  in  Table  5-12.   Based  on  these  data,  the  solutes 
in  this  study  do  not  show  an  appreciable  amount  of  sorption 
non-equilibrium.   These  solutes  are  not  strongly  sorbed,  and 
thus  are  not  expected  to  exhibit  a  large  degree  of  non- 
equilibrium.   Regression  of  R   values  with  ISNE  values 

a 

2 

yielded  a  model  with  a  r   value  of  0.57.   Statistical 

analysis  of  this  regression  model  indicated  that  the  model 

was  useful  for  predicting  ISNE  from  R   values  at  the  0.01 

a 

significance  level. 

The  cause  of  nonequilibrium  may  be  the  result  of 
several  physical  or  chemical  phenomena  which  limit  the  rate 


90 


Table  5-12.   An  empirical  index  of  sorption  nonequilibrium 
(ISNE)  for  12  selected  aromatic  solutes 
leaching  through  Lake  Alfred  aquifer  material. 


Compound 


ISNE 


Benzene 

Toluene 

Ethylbenzene 

m,p-Xylene 

o-Xylene 

Isopropylbenzene 

n-Propylbenzene 

3  or  4  Ethyltoluene 

1/3,5-Trimethylbenzene 

2-Ethyltoluene 

1/2,4-Trimethylbenzene 

1,2,3-Trimethylbenzene 


0.20 
0.23 
0.14 
0.17 
0.17 
0.13 
0.10 
0.12 
0.11 
0.10 
0.13 
0.13 
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of  sorption.   Physical  limitations  to  equilibrium  include 
diffusion  controlled  adsorption-desorption  processes  (Rao 
and  Davidson,  1979)  or  the  presence  of  physical  barriers 
limiting  the  interaction  of  the  sorbent  and  solute  (eg.  soil 
aggregates,  surface  films).   The  kinetics  of  chemical 
reactions  between  the  sorbent  and  solute  may  be  limiting, 
thereby  explaining  the  nonequilibrium  in  the  column 
breakthrough  curves.   Multisite  models  have  been  proposed  to 
account  for  sorption  nonequilibrium  (Rao  et  al . ,  1979). 
However,  the  physical  and  chemical  processes  are 
mathematically  equivalent  when  written  in  non-dimensional 
form,  thus  the  identification  of  the  process  responsible  for 
the  observed  nonequilibria  is  not  possible  from  breakthrough 
curve  data. 

Sorption  nonequilibria  are  also  a  function  of  the  flow 
velocity.   The  observed  effect  of  increased  flow  velocity  is 
displacement  of  the  elution  curve  towards  a  smaller 
breakthrough  volume,  whereas  the  calculated  effect  of 
increased  flow  velocity  is  a  broadening  of  the  elution 
curve,  owing  to  increased  dispersion,  without  a  change  in 
the  position  of  the  position  of  the  BTC.   This  effect  was 
demonstrated  by  Schwarzenbach  and  Westall  (1981)  where 
retardation  factors  decreased  with  increased  flow  rates, 
indicating  slow  sorption  kinetics. 

Rao  and  Davidson  (1979)  illustrated  that  the 
retardation  factor  may  also  be  a  function  of  concentration, 
with  increased  concentration  of  solutes  leading  to  decreased 
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sorption.   This  results  from  nonlinearity  of  the  sorption 
isotherm,  which  governs  the  position  of  the  breakthrough 
curve . 

The  BTC  data  in  this  study  were  obtained  for  a  single 
flow  velocity  and  at  only  one  range  of  concentrations. 
However,  the  importance  of  increased  flow  velocity  and  high 
concentration  on  the  transport  of  contaminants  may  be 
significant  at  the  Lake  Alfred  field  site.   Killan  (1987) 
reported  seepage  velocities  of  1  to  18  feet  per  day  in  the 
Lake  Alfred  aquifer.  This  high  flow  velocity,  combined  with 
the  high  levels  of  hydrocarbon  contamination  and  oxygen 
limitation  in  certain  areas  at  the  Lake  Alfred  site  suggests 
that  movement  of  solutes  may  be  more  rapid  than  these  column 
experiments  would  predict.   The  sorption  isotherms 
calculated  in  this  study  were  based  on  linear  isotherms,  for 
concentrations  which  were  present  in  the  well  water. 
However,  it  is  possible  that  residual  gasoline  in  the 
aquifer  may  provide  higher  concentrations  in  selected  areas, 
leading  to  more  rapid  leaching  of  aromatic  solutes  owing  to 
higher  concentrations  (eg.,  the  sorption  isotherms  may  be 
nonlinear,  and  R  may  be  concentration  dependent). 
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5-5.3   Comparison  of  Column  and  Equilibrium  Isotherm  Data. 


Retardation  factors  (Table  5-11  )  and  sorption 
coefficients  (Tables  5-10  and  5-7)  from  the  column  data  and 
the  equilibrium  data  isotherm  compare  favorably.   On 
average,  K   values  from  regression  analysis  of  isotherm  data 
using  the  linear  model  overestimate  the  value  of  the 
sorption  coefficients  by  40%,  and  the  Freundlich  model 
underestimates  the  value  of  the  sorption  coefficients  by  33% 
relative  to  the  column  derived  sorption  coefficients. 

Based  on  this  comparison,  the  column  data  appear  to 
fall  within  a  range  of  values  bounded  by  the  batch  isotherm 
data.   In  the  following  discussion,  sorption  coefficients 
from  the  column  data  are  used  to  evaluate  various  sorption 
relationships.   Column  data  are  used  here  since  they  are 
equivalent  to  the  isotherm  data,  and  also  since  the  column 
data  provide  values  for  ethylbenzene  and  the  propylbenzenes 
which  were  not  evaluated  in  the  isotherm  studies  owing  to 
low  concentrations  in  the  well  water. 

The  retardation  values  of  aromatic  hydrocarbons  in 
these  experiments  are  relatively  low.   The  retardation 
values  from  the  column  studies  range  between  1.36  and  2.40. 
These  data  indicate  that  the  most  retained  solute  will 
continue  to  move  at  42%  the  rate  of  water  movement.   Thus 
solutes  may  be  expected  to  move  relatively  rapidly  through 
the  site. 
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5.6   Evaluation  of  Sorption  Models 


The  actual  mechanisms  through  which  sorption  retards 

the  movement  of  solutes  are  not  well  known.   Various 

conceptual  models  are  available  to  help  describe  these 

sorption  process.   These  include  the  partitioning  between 

organic  matter  on  the  aquifer  matrix  (Chiou  et  al .  ,  1979, 

Karickhoff  et  al . ,  1979),  interactions  with  the  mineral 

surfaces  (Sabljic,  1987),  and  solvophobic  theory  (Rao  et 

al .  ,  1985).   To  assess  the  significance  of  these  models, 

sorption  data  from  the  column  experiments  were  compared  to 

several  theoretical  models.   The  column  data  were  normalized 

to  the  organic  carbon  content  of  the  Lake  Alfred  aquifer. 

These  values  of  K    are  shown  in  Table  5-10. 

oc 

5.6.1   Relationship  between  K    and  K 

■ &w &e- 

Several  authors  cite  the  linearity  of  the  sorption 
isotherms  as  evidence  of  the  dominance  of  the  partitioning 
mechanism  (Chiou  et  al . ,  1979,  Chiou  et  al.,  1983).   However 
this  evidence  may  be  suspect,  if  the  range  of  concentrations 
is  far  removed  from  the  maximum  solubility  of  the  compounds. 
In  addition,  many  sorption  models  are  indistinguishable  over 
sufficiently  small  concentration  ranges  (Curtis  et  al . , 
1986).   In  this  study,  solute  concentrations  are  far  below 
the  solubility  limits.   The  least  soluble  compounds  in  this 
work  are  3  ethyltoluene  and  4-ethyl tol uene  with  aqueous 
solubilities  of  40  mg/L.   However,  the  maximum  concentration 
employed  in  this  study  is  935  ug/L,  which  is  only  2%  of  the 
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maximum  solubility  level.   Therefore,  the  fact  that  the 

isotherms  are  linear  in  this  study  does  not  confirm  the 

dominance  of  the  partitioning  theory. 

An  improved  method  to  assess  the  importance  of 

partitioning  in  the  sorption  process  at  the  Lake  Alfred  site 

is  to  compare  Kqc  values  from  this  work  with  octanol-water 

partition  coefficients  (Figure  5-8)  and  water  solubility 

(Figure  5-9)  from  the  literature.   The  regression 

coefficients  for  these  correlations  are  shown  in  Table  5-13. 

These  experimentally  derived  relationships  can  now  be 

compared  with  those  from  previous  studies. 

The  experimentally  derived  relationships  between  K 

oc 

an<3  K   /  and  K    and  WS  were  determined  by  reqression  of 

ow         oc  ■*     3 

literature  values  of  K    and  water  solubility   with  K    data 

ow  J         oc 

from  the  sorption  BTC  data.   These  relationships  are 


log  Kqc   a   0.31  *  log  Kow  +  1.91 


log  K 


oc 


-0.272  *  log  WS  (umoles/L)  +  3.1k 


[5.2] 
[5.3] 


Table  5-14  compares  the  predicted  values  of  K    from 

oc 

the  work  of  several  authors  (Karickhoff  et  al  .  ,  1979,  Means 
et  al  .  ,  1982,  Chiou  et  al . ,  1983,  Kenaga  and  Goring,  1980, 
Briggs,  1981).   The  Kqc  data  from  this  study  (Table  5-10) 
consistently  fall   within  the  upper  range  of  the  predicted 
values  shown  in  Table  5-14.   The  relationships  used  to 
calculate  the  values  in  Table  5-14  were  based  on  a  wide 
range  of  organic  solutes  and  natural  sorbent  materials.   It 
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Table  5-13, 


Regression  coefficients  for  plots  of  log 
Koc  VS'  log  Kow  and  log  Koc  Vs'  log  WS- 


Log  It    vs.  Log  K 

,.  oc       3   ow 


Log  K    vs.  Log  WSC 


Slope 

Std.  error 
of  slope 

Y-intercept 

Std.  error  of 
y-intercept 

2 

r 

Number  of 
observations 

Degrees  of 
freedom 


0.  310 

0.052 
1.909 

0.073 
0.857 

8 


Slope 

-0 

.272 

Std.  error 
of  slope 

0 

.039 

Y-intercept 

3 

785 

Std.  error  of 
y-intercept 

0. 

074 

2 

r 

0. 

828 

Number  of 
observations 

12 

Degrees  of 
freedom 

10 

umoles/L 
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is  notable  that  K    from  this  study  falls  in  the  upper  range 

of  the  predicted  K    values,  given  the  paucity  of  organic 

carbon  in  the  Lake  Alfred  aquifer  material.   A  comparison  of 

predicted  K    values  calculated  from  water  solubility 

relationships  of  several  authors  and  that  of  equation  [5.3] 

also  show  that  the  experimental  K    values  from  this  study 

oc  * 

are  generally  in  the  upper  range  of  these  predicted  values 
as  well. 

The  relatively  high  values  of  K    oredicted  from 

oc  *• 

equations  [5.2]  and  [5.3]  may  be  the  result  of  several 

factors:  1)  error  in  the  measurement  of  organic  carbon,  2) 

increased  hydrophobicity  of  the  organic  matter  at  Lake 

Alfred  compared  to  the  referenced  studies,  3)  sorption  to 

the  mineral  surface,  or  4)  any  combination  of  the  above 

(Curtis  et  al . ,  1986).   It  is  unlikely  that  the  organic 

matter  in  this  study  is  more  hydrophobic  than  that  used  by 

other  researchers.   This  is  confirmed  by  analysis  of  the 

slopes  of  the  regression  lines  in  the  relationships  between 

K    and  K   .   The  slope  of  a  plot  of  K    vs.  K    may  be 
oc       ow  c  t  oc       ow    1 

viewed  as  a  measure  of  the  hydrophobicity  of  the  organic 
phase  in  the  partitioning  model  (Leo  et  al . ,  1971).   In  the 
experimental  data  presented  here  the  slope  of  equation  [5.2] 
was  less  than  that  observed  in  previous  studies  where 
partitioning  is  thought  to  predominate.   This  is  presented 
graphically  in  Figure  5-10.   Therefore,  increased 
hydrophobicity  of  the  organic  matter  at  Lake  Alfred  was 
ruled  out  . 
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It  is  more  likely  that  the  increased  sorption  results 

from  some  affinity  of  the  aromatic  solutes  in  this  study  for 

the  mineral  surface  of  the  aquifer  material.   This  is 

consistent  with  the  findings  of  Schwarzenbach  and  Westall 

(1981)  and  Curtis  et  al  .  ,  (19-86).   These  authors  demonstrate 

that  the  mineral  surface  area  and  the  nature  of  the  mineral 

surface,  exert  a  greater  influence  on  sorption  than  organic 

carbon  for  sorbents  with  low  amounts  of  naturally  occurring 

organic  material. 

5.6.2   Relationship  Between   X  and  K 

■* &e- 

To  assess  the  contribution  of  the  mineral  surface  in 

the  sorption  of  aromatic  solutes  from  the  Lake  Alfred 

aquifer,  K    values  from  this  study  were  correlated  with 

first  order  molecular  connectivity  indices  (  X).   This 

relationship  is  shown  in  Figure  5-11  and  the  regression 

parameters  are  shown  in  Table  5-15.   The  use  of  this 

correlation  was  based  on  the  suggestion  of  Milgelgrin  and 

Gerstl  (1983)  that  molecular  structure  or  topology  may  be 

more  effectively  correlated  with  sorption  than  K    or  WS . 

ow 

Sabljic  (1987)  suggested  the  use  of  first  order  molecular 

connectivity  indices  (  X)  as  an   estimator  of  molecular 

topology.   The  regression  equation  developed  by  Sabljic 

(1987)  is  also  shown  in  Figure  5-11.   This  relationship  was 

based  on  the  regression  of  calculated   X  versus  literature 

values  of  K    data, 
oc 

The  slopes  of  these  lines  are  not  significantly 
different  at  the  0.05  significance  level.   The  correlation 
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Table  5-15.   Regression  coefficients  for  the  relationship 

between  log  K    and  1X. 
3   oc 


Slope  0.360 

Std.  error  of  slope     0.050 
Y-intercept  1.509 


Std.  error  of 
y-intercept  0.072 

„2 


0.839 


Number  of 
observations  12 

Degrees  of  freedom     10 
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coefficients  are  also  comparable:  0.916  for  this  study  vs 
0.976  for  Sabljic  (1987).   Sabljic  (1987)  demonstrated  that 
first  order  molecular  connectivity  was  a  quantitative 
measure  of  the  area  occupied  by  the  projection  of  the  non- 
hydrogen  skeleton  of  a  molecule.   The  goodness  of  fit 
between   X  and  K    data  in  this  study  supports  the 
hypothesis  that  sorption  depends,  at  least  in  part,  on  some 
type  of  surface  interaction. 

However,  comparison  of  correlation  coefficients  between 

the  three  models  (K   ,  WS  and   X)  indicates  that  neither  of 

ow 

these  models  completely  describe  the  sorption  process  (see 
Tables  5-13  and  5-15).   This  suggests  that  the  sorption 
mechanism  is  in  reality  a  combination  of  processes, 
interacting  to  yield  an  overall  sorption  effect.   The  lack 
of  any  dominant  mechanism  may  be  more  pronounced  in  this 
study  resulting  from  the  low  organic  carbon  content  of  the 
Lake  Alfred  aquifer  material.   This  serves  to  reduce  the 
partitioning  effect,  by  eliminating  the  sorption  substrate 
(organic  carbon).   In  addition,  Schwarzenbach  and  Westall 
(1981)  demonstrated  that  organic  poor  sorbents  with  high 

specific  surface  areas  may  still  exhibit  small  K .,  values, 

d 

indicating  that  surface  interactions  alone  did  not 
completely  account  for  sorption.   These  data  support  the 
observation  of  Voice  and  Weber  (1983)  that,  given  the 
heterogeneous  nature  of  natural  sorbent  materials,  sorption 
mechanisms  of  organic  solutes  in  the  environment  probably 
involve  many  types  of  interactions.   The  importance  of  a 
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given  reaction  mechanism  depends  on  the  nature  of  the 
sorbent  surface.   Partitioning  is  probably  more  important  in 
soils  with  high  organic  carbon  contents.   The  varying 
degrees  of  sorption  in  soils  with  low  organic  carbon 
contents  reported  by  Milgelgrin  and  Gerstl  (1984)  reflect 
the  variation  in  the  ability  of  mineral  surfaces  to  sorb 
organic  compounds. 


5.7  Comparison  of  Mixed  Solute  and  Single  Solute 

Retardation  . 


Nkedi-Kizza  et  al .  (1987)  demonstrated  the  influence  of 
organic  co-solvents  on  the  movement  of  hydrophobic  organic 
compounds  through  soils.   In  this  dissertation,  the  study 
compounds  were  a  multicomponent  mixture  of  dissolved 
aromatic  solutes,  resulting  from  the  partial  solubilization 
of  a   leaded  gasoline  product  into  groundwater  at  the  Lake 
Alfred  site.   The  presence  of  multiple  solutes  may  affect 
the  sorption  of  a  single  component  of  the  mixture  either  by 
changing  the  solubility  of  the  component  or  through 
competitive  sorption  (Brookman  et  al . ,  1985).   To  evaluate 
this  possibility,  a  single  solute  (benzene,  4  mg/L  dissolved 
in  RAP-2  well  water)  was  passed  through  a  soil  column.   The 
breakthrough  of  this  solute  is  shown  in  Figure  5-12. 
Evaluation  of  the  retardation  factor  for  this  column  yielded 
an  R  value  of  1.4,  which  is  equivalent  to  the  R  value  for 
benzene  from  the  mixed  solute  sample.   Based  on  these  data, 
no  co-solute  effect  on  benzene  was  observed.   If  a 
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competitive  sorption  effect  is  operating  in  the  mixed  solute 
sample,  then  benzene  as  a  single  solute  should  show 
increased  sorption.   Lack  of  sorption  increase  indicates 
that  competitive  sorption  between  solutes  was  not  important 
in  this  experimental  system  although  this  may  not  be  the 
case  in  portions  of  the  aquifer  with  residual  concentrations 
of  gasoline. 

5 .8   Evaluation  of  Hydrogen  Peroxide  Reactivity 


Hydrogen  peroxide  is  known  to  be  a  viable  method  of 
increasing  the  dissolved  oxygen  in  aquifer  systems  (Britton, 
1985,  TRI  Report,  1982).   The  purpose  of  this  experiment  was 
to  evaluate  the  reaction  kinetics  and  the  extent  of 
conversion  of  hydrogen  peroxide  to  O   in  the  Lake  Alfred 
aquifer  system. 

Initial  experiments  with  distilled  deionized  water  and 
known  additions  of  dilute  hydrogen  peroxide  showed  no 
increases  in  dissolved  oxygen.   Even  after  addition  of  5  mL 
of  50%  hydrogen  peroxide  to  the  reaction  flask,  no  increase 
in  dissolved  oxygen  was  noted  over  a  45  minute  interval. 
The  titer  of  the  50%  stock  solution  (49.7%)  was  confirmed  by 
titration  with  potassium  permanganate.   These  data  indicated 
the  stability  of  hydrogen  peroxide  in  the  absence  of  a 
catalyst . 

The  reaction  of  hydrogen  peroxide  in  filter  sterilized 
Lake  Alfred  water  was  investigated  to  assess  the 
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availability  of  non-biological  catalysts  in  the  aqueous 
phase.   Water  from  wells  OHM-4  (with  aromatic  solutes)  and 
RAP-2  (no  aromatic  solutes)  were  titrated  with  hydrogen 
peroxide  solutions  of  240  and  2400  mg/L  with  no  apparent 
increase  in  dissolved  oxygen,  indicating  an  the  absence  of  a 
H„0„  active  catalyst  in  both  these  water  samples. 

The  response  of  non-filtered  water  from  well  OHM-4  to 
the  addition  of  50%  hydrogen  peroxide  is  shown  in  Figure  5- 
13.   The  pH  of  the  water  in  the  reaction  flask  was  6.5  and 
this  value  remained  constant  throughout  the  course  of  the 
experiment.   The  redox  potential  increased  from  59  milli- 
volts ( mv )  to  338  mv  immediately  following  the  addition  of 
500  ul  of  50%  H  02  (equivalent  to  2000  mg/L  H^).   This 
addition  of  H  O      was  sufficient  to  maintain  an  increase  of  1 
mg/L  over  the  ambient  DO  in  the  reaction  flask.   The 
addition  of  sterile  aquifer  material  from  the  Lake  Alfred 
site  increased  the  DO  of  the  reaction  flask  immediately 
after  introduction.   This  indicated  the  catalytic  ability  of 
the  Lake  Alfred  aquifer  and  was  most  likely  associated  with 
the  presence  of  iron  salts  (Britton  1985)  although  iron 
concentrations  were  not  determined  for  the  Lake  Alfred 
aquifer  material.   This  small  scale  study  helped  to 
determine  the  reactivity  of  the  hydrogen  peroxide  in  the 
Lake  Alfred  aquifer  system,  and  provided  ranges  for  use  of 
hydrogen  peroxide  in  the  biological  experiments.   Gas 
chromatographic  analyses  of  aromatic  compounds  during  the 
course  of  this  experiment  showed  that  there  were  no 
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measurable  changes  in  the  concentrations  of  these  solutes 
resulting  from  the  addition  of  hydrogen  peroxide.   The 
absence  of  aromatic  hydrocarbon  removal  in  this  experiment 
should  not  imply  that  HO-  is  an  ineffective  method  of 
oxygen  augmentation.   The  dissolved  oxygen  levels  did 
increase  in  this  experiment.   However,  removal  of 
hydrocarbons,  abiotically,  via  oxidation  of  the  aromatic 
molecules  appeared  not  to  be  significant.   This  experiment 
was  not  designed  to  measure  the  microbial  removal  of 
hydrocarbons.   The  time  scale  of  this  experiment  was  too 
short  to  observe  a  microbial  effect.   Batch  biodegradation 
experiments  (sections  5.8  and  5.9)  demonstrated  that  several 
days  were  required  for  microbial  adaptation  to  hydrogen 
peroxide  . 

The  data  from  this  experiment  simply  confirm  that  there 
are  sufficient  catalysts  available  to  mediate  the  conversion 
of  ti2°2    to    °2*   Application  of  H  02  to  the  Lake  Alfred 
aquifer  is  underway.   Preliminary  data  indicate  that 
dissolved  oxygen  levels  were  increased/  and  that 
hydrocarbons  concentrations  were  reduced  following  an 
adaptation  period. 

5 . 9   Batch  Biodegradation  Experiment  #1 


Two  separate  batch  biodegradation  experiments  were 
performed.   The  first  experiment  was  designed  to  investigate 
the  affect  of  various  combinations  of  hydrogen  peroxide  and 
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ammonium  chloride  on  the  biodegradation  of  the  aromatic 
compounds  in  the  Lake  Alfred  aquifer,  and  to  assess  the 
degradation  of  these  compounds  in  the  presence  of  dissolved 
oxygen.   The  average  concentrations  of  aromatic  hydrocarbons 
in  the  microcosms  of  experiment  #1  over  the  time  course  of 
the  experiment  (31  days)  are  shown  in  Table  5-16.   A 
detailed  presentation  and  statistical  analysis  of  these  data 
are  presented  in  Appendix  E.   These  data  were  fit  to  several 
rate  equations.  Zero  order,  first  order,  an  empirically 
based  first  order  rate  equation  (Thomas-slope  method), 
second  order,  and  a  mixed  order  (zero  to  first  order)  rate 
equations  were  fit  to  the  biodegradation  data.   Only  the 
first  order  rate  equations  gave  adequate  fit  to  the  data, 
based  on  an  analysis  of  the  coefficients  of  determination 
for  the  various  models  (Global  F-test ,  0.05  significance 
level).   The  results  of  linear  regression  analysis  of  the 
data  to  the  first  order  models  are  shown  in  Table  5-17 
(first  order)  and  Table  5-18  (Thomas-slope).   In  both  tables 
the  rate  constants,  calculated  half  lives  and  the 
coefficients  of  determination  for  each  solute  under  each 
treatment  condition  are  presented. 

Regression  analyses  and  rate  data  for  both  types  of 
first  order  rate  equations  are  presented,  since  neither 
method  yields  coefficients  of  determination  which 
consistently  provide  superior  fit  to  the  data.   First  order 

and  Thomas  slope  rate  equations  each  provide  adequate  fit  to 

2 
the  data  as  evidenced  by  the  r   values.  In  the  sections 
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Table  5-17. 


Biodegradation  rate  constants,  half -lives  and 
correlation  coefficients  for  the  fit  of 
biodegradation  experiment  #1  data  to  a 
first  order  rate  equation. 


Treatment   Benzene   Toluene   m,p-Xylene    o-Xylene   3,4-ET£ 


1A 

k 

2 

r 

0.233 
2.97 

0.219 
3.17 

0.295 
2.35 

0.282 
2.46 

0.235 
2.95 

0.779 

0.844 

0.993 

0.995 

0.989 

IB 

k 

2 

r 

0.233 
2.97 

0.262 
2.65 

0.190 
3.65 

0.208 
3.33 

0.135 
4.47 

0.887 

0.901 

0.909 

0.876 

0.928 

1C 

k 

th 

2 
r 

0.242 
2.86 

0.261 
2.66 

0.271 
2.56 

0.284 
2.44 

0.227 
3.05 

0.970 

0.994 

0.913 

0.975 

0.967 

ID 

k 
th 

2 
r 

0.212 
3.27 

0.254 
2.73 

0.216 
3.21 

0.183 
2.79 

0.157 
4.42 

0.942 

0.920 

0.753 

6.868 

0.921 

IE 

k 

th 

2 
r 

0.154 
4.50 

0.200 
3.47 

0.200 
3.47 

0.076 
9.17 

0.171 
4.05 

0.992 

0.952 

0.753 

0.969 

0.936 

IF 

k 

th 

2 
r 

0.038 
18.3 

0.107 
6.48 

0.185 
3.75 

0.036 
19.47 

0.134 
5.17 

0.885 

0.945 

0.706 

0.917 

0.897 

1G 

k 

th 

2 

0.010 
69.31 

0.013 
51.7 

0.020 

35.55 

0.017 
40.29 

0.030 
23.42 

0.370 


0.401 


0.625 


0.584 


0.724 


IB    k     0.145     0.102      0.193        0.126 

5.50 


1C    k     0.097     0.146      0.214        0.133 

5.21 


k 

0.223 
3.11 

0.207 
3.35 

0.249 
2.78 

2 
r 

0.935 

0.982 

0.989 

k 

0.145 
4.78 

0.102 
6.80 

0.193 
3.59 

2 
r 

0.936 

0.937 

0.939 

k 
th 

0.097 
7.17 

0.146 
4.75 

0.214 
3.24 

2 

r 

0.927 

0.994 

0.777 

k 
th 

0.034 
20.21 

0.052 
13.30 

0.172 
4.03 

2 
r 

0.775 

0.982 

0.632 

k 

0.053 
13.05 

0.065 
10.75 

0.164 
4.23 

2 

r 

0.946 

0.946 

0.591 

k 

0.043 
16.16 

0.036 
19.42 

0.183 
3.79 

2 

r 

0.820 

0.868 

0.673 

k 
2 

0.030 
23.50 

0.026 
26.36 

0.030 
23.11 

a3,4-Ethyltoluene 
1 , 3 , 5-Trimethylbenzene 

Q 

2 -Ethyl toluene 

1,2, 4-Trimethylbenzene 

e 
1,2, 3-Trimethylbenzene 

Dissolved  oxygen 


0.230 
3.01 

0.986  0.013 


0.933  0.093 


0.990  0.111 


ID        k  0.034  0.052  0.172  0.033 

21.00 


0.875  0.0 


IE         0.053     0.065      0.164        0.052 

13.46 


0.931     0.005 

IF    k     0.043     0.036      0.183        0.032 

21.66 

0.799  0.042 

1G         0.030     0.026      0.030        0.023 

30.27 

0.725  0.683  0.728  0.655  0.429 
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Table  5-10.   Biodegradation  rate  constants,  half -lives  and 
correlation  coefficients  for  the  fit  of 
biodegradation  experiment  #1  data  to  the 
Thomas  slope  rate  equation. 


Treatment   Benzene   Toluene   m,p-Xylene    o-Xylene   3,4-ET 


1A 

t^b 

0.224 

0.239 

0.188 

0.139 

0.195 

3.09 

2.91 

3.68 

4.97 

3.56 

2 

r 

0.964 

0.958 

0.974 

0.999 

0.967 

IB 

k 

-3.406 

-0.155 

0.135 

-0.155 

0.111 

-0.204 

-4.46 

5.14 

-0.107 

6.27 

2 
r 

0.764 

0.129 

0.838 

0.129 

0.745 

1C 

k 

0.094 

0.153 

0.183 

0.087 

0.195 

7.37 

4.54 

3.79 

7.93 

3.55 

2 
r 

0.899 

0.986 

0.987 

0.947 

0.984 

ID 

k 

0.118 

0.171 

0.228 

0.096 

0.205 

t* 

5.90 

4.06 

3.04 

7.21 

3.38 

2 
r 

0.978 

0.975 

0.966 

0.712 

0.972 

IE 

k 

-0.027 

0.048 

0.144 

0.045 

0.151 

th    - 

■25.99 

14.30 

4.80 

15.33 

4.60 

2 

r 

0.040 

0.199 

0.946 

0.926 

0.982 

IF 

k 

0.101 

0.142 

0.116 

0.130 

0.160 

th 

*>* 

6.88 

4.89 

5.96 

5.31 

4.33 

2 
r 

0.568 

0.989 

0.822 

0.816 

0.983 

1G 

k 

0.130 

0.177 

0.176 

0.157 

0.173 

2 

5.33 

3.92 

3.93 

4.41 

4.00 

0.488     0.780      0.863        0.761     0.883 
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Table  5-18.   Continued, 


Treatment   1,3,5-TMB   2-ET    1,2,4-TMB    1,2,3-TMB    DO 


1A 

k 

2 

r 

0.140 
4.96 

0.165 
4.20 

0.193 
3.59 

0.140 
4.94 

0.382 
1.81 

0.988 

0.955 

0.978 

0.984 

0.990 

IB 

k 

2 
r 

0.032 
21.66 

0.012 
56.98 

0.172 
4.04 

0.078 
8.85 

0.259 
2.68 

0.139 

0.012 

0.925 

0.526 

0.952 

1C 

k 

2 

r 

0.158 
4.39 

0.129 

5.38 

0.197 
3.52 

0.122 
5.70 

0.262 
2.64 

0.975 

0.949 

0.986 

0.903 

0.955 

ID 

k 
tJj 

2 

r 

0.166 
4.17 

0.106 
6.57 

0.235 
2.96 

0.148 
4.67 

0.310 
2.24 

0.969 

0.889 

0.964 

0.933 

0.939 

IE 

k 

0.127 
5.45 

0.102 
6.81 

0.210 
3.44 

0.110 
6.30 

0.324 
2.14 

2 
r 

0.896 

0.786 

0.985 

0.785 

0.940 

IF 

k 
t% 

0.176 
3.95 

0.151 
4.60 

0.167 
4.15 

0.169 
4.10 

0.074 
9.37 

2 

r 

0.924 

0.850 

0.981 

0.923 

0.183 

1G 

k 

2 

r 

0.171 
4.05 

0.176 
3.95 

0.170 
4.07 

0.174 
4.00 

0.176 
3.94 

0.897 

0.890 

0.896 

0.885 

0.801 

day 
days 
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below,  each  treatment  is  discussed  individually/  prior  to  an 
overall  analysis  of  these  experiments. 

5.9.1  Treatment  1A. 

Data  are  plotted  in  Figures  5-14  and  5-15.   The 
dissolved  oxygen  in  these  microcosms  was  not  artificially 
increased,  other  than  by  aeration  during  the  transfer  and 
filling  of  the  solute  containing  water  into  the 
biodegradation  vials.   The  DO  of  these  vials  at  time  =  0  was 
7.3  mg/L .   An  examination  of  the  half  lives  of  these  solutes 
showed  that  benzene  (3.1  days)  and  toluene  (2.9  days)  were 
readily  removed  from  the  microcosm  compared  to  a  half  life 
of  4.96  for  1 , 3 , 5-trimethylbenzene  .   This  is  in 
contradiction  to  studies  which  note  the  recalcitrance  of 
these  compounds  to  biodegradation  (Bossert  and  Bartha, 
1984).   The  ortho  isomer  of  xylene  was  more  resistant  to 
microbial  attack  than  were  the  meta  and  para  isomers.   This 
confirmed  the  data  of  Kappeler  and  Wuhrmann  (1978a,  1978b). 
Complete  degradation  of  all  solutes  was  achieved  by  31  days 
(detection  limit  0.5  ug/L)  and  in  many  cases  degradation  was 
complete  in  15  days.   Toluene  was  degraded  particularly 
rapidly.   The  depletion  of  oxygen  in  these  microcosms 
suggested  that  this  loss  was  microbially  mediated  (Figure  5- 
16)  . 

5.9.2  Treatment  IB. 

Data  for  this  treatment  are  shown  in  Appendix  E  . 
These  data  demonstrate  the  effect  of  17  mg/L  hydrogen 
peroxide  treatment  on  the  degradation  of  the  aromatic 
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solutes.   The  main  feature  of  these  plots  is  the  eight  day 
lag  phase  in  the  removal  of  benzene,  toluene  and  o-xylene. 
However/  with  time/  degradation  of  these  compounds  was 
essentially  complete.   The  CgH,2  compounds  also  displayed  a 
lag  up  to  eight  days  in  length.   The  half  lives  were 
slightly  increased  over  treatment  1A/  indicating  the  time 
involved  in  the  adaption  of  the  microorganisms  to  the 
hydrogen  peroxide. 

5.9.3  Treatment  1C 

The  effects  of  increasing  the  concentration  of  H0_  to 
68  mg/L  was  shown  in  this  treatment.   There  was  no  apparent 
increase  in  toxicity  over  the  17  mg/L  treatment.   Benzene/ 
toluene  and  m/p-xylene  were  completely  removed.   Again, 
1/ 2 /4-trimethylbenzene  exhibited  the  most  rapid  degradation 
of  the  CgH,~  compounds.   The  initial  DO  in  these  microcosms  . 
was  10.5  mg/L.   The  increased  half  lives  of  the  study 
compounds  reflect  the  adaptation  to  hydrogen  peroxide  noted 
in  treatment  2B.   Dissolved  oxygen  concentrations  are  shown 
in  Figure  5-16. 

5.9.4  Treatment  ID 

Ammonium  chloride  (18mg/L)  was  added  to  these  vials. 
This  concentration  was  chosen  based  on  the  calculated  amount 
of  nitrogen  required  to  completely  degrade  the  aromatic 
solutes  assuming  a  C:N  ratio  of  10:1.   The  initial  DO 
concentration  was  9  mg/L.   Oxygen  consumption  appeared  rapid 
but  fell  after  seven  days.   Addition  of  ammonium  chloride 
significantly  reduced  the  rate  of  microbial  degradation  of 
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benzene,  toluene,  1 , 3 , 5-trimethylbenzene ,  2-ethyltoluene , 
and   1 , 2 , 3-trimethylbenzene .   Both  m,p-xylene  (half  life  =  3 
days)  and  1 , 2 , 4-trimethylbenzene  (half  life  =  2.96  days) 
were  less  affected  and  were  each  degraded  to  1.7  ug/L .   At 
the  end  of  31  days,  31%  of  1 , 2 , 3-trimethylbenzene  remained 
but  1 , 3, 5-trimethylbenzene ,  2-ethyltoluene  and  1,2,4- 
trimethylbenzene  although  not  completely  degraded  were  well 
removed . 

5.9.5  Treatment  IE 

Graphical  representation  of  this  treatment  is  shown  in 
Appendix  E.   These  vials  are  treated  with  18  mg/L  NH  CI  and 
17  mg/L  H^O-.   This  combination  of  treatments  caused  an 
increase  in  the  1/2  lives  for  the  Cfi-CR  compounds,  but  the 
removal  rates  of  the  CgH-,_  compounds  were  improved  relative 
to  treatment  with  ammonium  chloride  alone.   Ultimate  removal 
of  compounds  was  good  except  for  the  more  recalcitrant  o- 
xylene,  1 ,  2 , 3-trimethy lbenzene ,  3,4-ET  and  1,3,5- 
trimethylbenzene .   The  lag  in  degradation  of  the  aromatic 
solutes  was  reflected  in  the  reduced  consumption  of  DO  in 
the  microcosms  (Figure  5-17). 

5.9.6  Treatment  IF 

Treatment  with  68  mg/L  HO   and  18  mg/L  NH  CI  is  shown 
in  Appendix  E  and  average  data  are  shown  in  Table  5-16. 
There  was  an  obvious  lag  in  the  DO  profile  (Figure  5-17), 
and  the  half  life  for  dissolved  oxygen  was  increased  from 
1.8  days  in  treatment  1A  to  9.4  days  in  this  treatment. 
This  was  also  reflected  in  the  concentrations  of  benzene,  o- 
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xylene,  2-ethyl toluene ,  1 ,  3 , 5-trimethylbenzene  and  1,2,3- 

trimethylbenzene.   All  compounds  showed  a  substantial 

increase  in  half  lives,  except  for  m,p-xylene  and  1,2,4- 

trimethylbenzene .   Throughout  this  study  these  compounds 

were  well  degraded.   These  compounds  were  the  most  rapidly 

removed  compounds  in  the  work  of  Kappeler  and  Wuhrmann 

(1978a,  1987b).   The  changes  in  the  degradation  rates  for 

the  other  aromatic  compounds  may  reflect  a  change  in  the 

community  structure  of  the  well  water  bacteria.   The 

bacteria  from  the  field  site  may  be  adapted  to  degrade 

aromatic  compounds  that  are  usually  thought  to  be 

recalcitrant  or  degraded  slowly  (i.e.,  benzene).   In 

treatment  1A,  benzene  was  rapidly  and  easily  removed  by  the 

bacterial  community.   However,  by  disrupting  this  community 

by  the  addition  of  additional  nutrients  and  hydrogen 

peroxide,  only  the  compounds  which  are  more  easily 

biodegraded  are  removed.   It  may  be  possible  that  if  these 

experiments  were  carried  out  for  longer  incubation  periods, 

the  microbial  communities  may  have  adapted  to  efficiently 

degrade  the  remaining  aromatic  compounds. 

5.9.7   Treatment  1G 

This  treatment  was  a  sterile  control.   The  DO  in  this 

system  remained  at  the  initial  levels  (Figure  5-17), 

indicating  the  sterility  of  the  system.   Losses  from  the 

system  are  in  the  range  of  25-50%  for  Cr-co  solutes,  and  55- 

6   8 

65%  for  CgH12  solutes.   These  losses  were  most  likely  the 
result  of  diffusion  of  the  volatile  components  through  the 
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teflon  septa.   The  loss  of  volatile  compounds  was 
exacerbated  by  storage  at  20  C.   However,  this  does  not 
account  for  the  increased  reduction  of  C_H10  concentrations 
relative  to  the  more  volatile  C,-CQ.   C^-Cn  values  would 

DO        DO 

normally  exhibit  greater  volatile  losses.   The 
biodegradation  data  were  not  corrected  for  these  losses. 


5.10   Batch  Biodegradation  Experiment  #2 


The  second  series  of  batch  biodegradation  experiments 
were  performed  to  assess  the  efficacy  of  several  methods  of 
oxygen  addition  to  the  Lake  Alfred  aquifer,  and  to  repeat 
some  of  the  previous  treatments.   In  this  series  of 
experiments,  additional  sterile  controls  were  added  to 
better  assess  the  losses  exhibited  in  biodegradation 
experiment  #1  (treatment  IF).   In  addition,  microbial 
activity  was  measured  by  quantifying  the  microbial  reduction 
of  INT  to  INT-formazan. 

The  data  from  this  experiment  were  fit  to  several  rate 
models  as  described  in  section  5.9.   No  single  model  gave  a 
consistently  good  fit  to  the  data.   Mixed  order  (zero  order 
to  first  order)  and   zero  order  rate  equations  did  not  match 
the  data  well  as  evidenced  by  low  coefficients  of 
determination . 

Average  values  of  hydrocarbons  in  the  microcosms  are 
shown  in  Table  5-19.   Rate  constants,  and  regression 
coefficients  for  the  fit  of  biodegradation  #2  data  to  the 
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first  order  model  are  shown  in  Table  5-20  and  via  the  Thomas 
slope  method  in  Table  5-21.   Examination  of  correlation 
coefficients  from  the  fit  of  the  Thomas  slope  model  to  the 
data  in  this  experiment  suggest  that  rate  constants  from 
this  model  selectively  give  an  adequate  estimate  of  the 
biological  degradation  constant. 

5.10.1  Treatment  2A 

This  is  a  repeat  of  treatment  1A.   The  general  trends 
were  the  same  and  the  rate  constants  were  equivalent.   Again 
m,p-xylene  and  1 , 2 ,4-trimethylbenzene  showed  rapid  loss  as 
noted  in  treatment  1A.   The  extent  of  degradation  was  not  as 
complete  as  in  treatment  1A.   This  may  be  the  result  of  the 
increased  concentration  of  toluene  (7759  ug/L  vs  2171  ug/L) 
and  the  concentration  of  m,p-xylene  (7241  ug/L  vs  4823  ug/L 
in  1A).   There  was  a  rapid  loss  in  DO  over  the  first  two 
days,  indicating  microbial  activity. 

5.10.2  Treatment  2B 

This  is  a  repeat  of  treatment  ID.   The  half  lives  were 
generally  higher,  reflecting  a  decrease  in  hydrocarbon 
removal.   There  was  a  pronounced  lag  phase  of  8  days  for 
benzene  (Figure  5-18).   The  lag  in  degradation  for  toluene, 
o-xylene,  m,p-xylene  and  2-ethyltoluene  was  two  days.   The 
DO  profile  in  this  treatment  was  consistent  with  that  of  2A 
indicating  C>2  removal  during  the  first  several  days. 
However,  microbial  activity  was  increased  over  treatment  2A. 
This  microbial  activity  did  not  result  in  significant 
degradation  of  hydrocarbons. 
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Table  5-20.   Biodegradation  rate  constants  and  correlation 
coefficients  for  the  fit  of  biodegradation 
experiment  #2  data  to  a  first  order  rate 
equation. 


Treatment   Benzene   Toluene    Ethbza   m,p-Xylene  o-Xylene 


r 


2A    k^    0.042     0.071      0.070 


2         0.462  0.486  0.344 


2B    k  0.049  0.089  0.082 

r2  0.887  0.942  0.919 

2C    k  0.003  0.006  0.013 

r2  0.105  0.317  0.619 

2D    k  0.043  0.073  0.096 

rZ  0.621  0.877  0.884 

2E    k  -0.004  0.006  0.017 

r2  0.169  0.219  0.138 


2F    k  -0.001  0.001  0.024 

r2  0.003  0.012  0.177 

2G    k  0.067  0.036  0.004 

r2  0.234  0.057  0.002 

2H    k  0.059  0.072  0.135 

r2  0.640  0.683  0.469 

21    k  -0.002  0.007  0.049 

r2  0.056  0.211  0.160 


0 

.046 

0 

.052 

0 

.618 

0 

.972 

0 

.114 

0 

.070 

0 

.955 

0 

.894 

0 

.013 

0 

.007 

0 

.619 

0 

403 

0 

138 

0 

046 

0 

986 

0 

800 

0 

017 

0 

002 

0 

009 

0 

038 

0. 

007 

0 

002 

0 

179 

0. 

029 

0. 

081 

0. 

057 

0. 

266 

0 

406 

0. 

095 

0. 

057 

0. 

398 

0. 

500 

0. 

010 

0, 

004 

0. 

873 

0. 

740 
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Table  5-20..   Continued 


Treatment    3,4-ET   1,3,5-TMB   2-ET     1,2,4-TMB   1,2,3-TMB 


2A 

k 
r2 

0.046 

0.049 

0.048 

0.049 

0.048 

0.698 

0.867 

0.917 

0.574 

0.923 

2B 

k 

0.105 

0.086 

0.073 

0.130 

0.072 

r2 

0.951 

0.889 

0.886 

0.949 

0.8  71 

2C 

k 

0.017 

0.015 

0.013 

0.016 

0.012 

r2 

0.653 

0.553 

0.557 

0.572 

0.539 

2D 

k 

0.078 

0.057 

0.053 

0.091 

0.045 

r2 

0.691 

0.885 

0.833 

0.586 

0.858 

2E 

k 

0.019 

0.010 

0.008 

0.022 

0.008 

r2 

0.360 

0.139 

0.136 

0.404 

0.128 

2F 

k 

0.013 

0.011 

0.009 

0.012 

0.006 

r2 

0.333 

0.240 

0.192 

0.256 

0.110 

2G 

k 

0.103 

0.096 

0.058 

0.085 

0.085 

r2 

0.531 

0.585 

0.578 

0.255 

0.763 

2H 

k 

0.091 

0.070 

0.067 

0.087 

0.056 

r2 

0.421 

0.529 

0.649 

0.297 

0.572 

21 

k2 
r 

0.024 

0.018 

0.021 

0.016 

0.014 

0.627 

0.560 

0.627 

0.614 

0.651 

Ethylbenzene 

day 
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Table  5-21.   Biodegradation  rate  constants  and  correlation 
coefficients  for  the  fit  of  biodegradation 
experiment  #2  data  to  the  Thomas  slope  rate 
equation. 


Treatment  Benzene   Toluene    Ethbz   m,p-Xylene  o-Xylene 


2A 

r2 

0.239 
0.915 

0.248 
0.929 

0.250 
0.948 

0.222 
0.893 

0.096 
0.962 

2B 

k 
r2 

0.004 
0.020 

0.081 
0.902 

0.194 
0.921 

0.133 
0.971 

0.028 
0.373 

2C 

k 
r2 

0.379 
0.850 

0.139 
0.384 

0.232 
0.920 

0.122 
0.802 

0.139 
0.466 

2D 

k 
r2 

0.032 
0.196 

-0.016 
0.026 

0.170 
0.975 

0.124 
0.997 

0.031 
0.249 

2E 

k 
r2 

0.597 
0.979 

0.153 
0.832 

0.256 
0.931 

0.201 
0.687 

0.298 
0.388 

2F 

k 
r2 

0.228 
0.631 

0.212 
0.730. 

0.249 
0.923 

0.237 
0.740 

0.212 
0.880 

2G 

$ 

0.230 
0.948 

0.232 
0.935 

- 

0.227 
0.938 

0.215 
0.908 

2H 

k 
r2 

0.098 
0.883 

0.174 
0.974 

0.174 
0.977 

0.246 
0.937 

0.158 
0.960 

21 

k 
r2 

: 

- 

0.287 
0.957 

0.123 
0.651 

- 
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Table  5-21,   Continued. 


Treatment    3 , 4 -ET   1,3,5 -TMB   2 -ET    1,2,4 -TMB   1,2,3 -TMB 


2A 

k 

0.216 

0.152 

0.130 

0.228 

0.162 

r2 

0.909 

0.960 

0.961 

0.914 

0.915 

2B 

k 

0.156 

0.161 

0.144 

0.172 

0.138 

r2 

0.942 

0.890 

0.863 

0.961 

0.843 

2C 

k 

0.162 

0.177 

0.154 

0.173 

0.157 

r2 

0.858 

0.850 

0.813 

0.852 

0.754 

2D 

> 

0.116 

0.129 

0.079 

0.140 

0.153 

0.987 

0.986 

0.910 

0.999 

0.880 

2E 

k 

0.223 

0.255 

0.254 

0.216 

0.255 

r2 

0.832 

0.881 

0.858 

0.801 

0.856 

2F 

k 

0.222 

0.234 

0.235 

0.230 

0.254 

r2 

0.852 

0.859 

0.786 

0.846 

0.704 

2G 

k 

0.225 

0.177 

0.211 

0.227 

0.206 

r2 

0.940 

0.961 

0.943 

0.940 

0.935 

2H 

k 

0.240 

0.149 

0.165 

0.251 

0.096 

r2 

0.946 

0.968 

0.989 

0.934 

0.718 

21 

k 

0.211 

0.217 

0.193 

0.206 

0.195 

r2 

0.905 

0.898 

0.841 

0.886 

0.871 

day 
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5.10.3   Treatments  2D  and  2E 


These  data  assessed  the  effect  of  Ho   addition  on  the 
microbial  community.   These  were  essentially  repeats  of 
treatments  1C  and  IF.   The  addition  of  58  mg/L  H  O 
(treatment  2D)  produced  a  2  day  lag  phase  for  benzene, 
toluene/  m,p-xylene  and  2-ethyltoluene  relative  to  no  HO 
treatment  (Figure  5-19).   The  extent  of  treatment  was 
comparable  to  treatment  2A  (air  addition).   The  DO  profiles 
are  shown  in  Figure  5-20.   The  lag  in  bioactivity  was 
paralleled  by  a  lag  in  consumption  of  DO  for  the  same  2  day 
period.   The  addition  of  NH  CI  in  treatment  2E  produced  a 
toxic  effect/  and  the  hydrocarbon  data  were  equivalent  to 
the  sterile  control  losses.   However/  the  DO  profile  (Figure 
5-20)  showed  consumption  of  dissolved  oxygen,  and  this 
implied  some  microbial  activity.   The  INT  data  also 
demonstrated  increased  bioactivity  (Figure  5-21)  following  a 
lag  of  at  least  two  days.  Microbial  activity  decreased  after 
14  days.   This  decreased  INT  reduction  was  not  noted  for 
treatment  2D.    Again,  there  was  microbial  activity  and  O 
consumption  without  substantial  reduction  in  hydrocarbon 
concentrations. 
5.10.4   Treatments  2G,  2H,  21. 

These  treatments  showed  the  effect  of  the  addition  of 
oxygen  gas  to  the  microcosms.   One  consequence  of  the  oxygen 
sparging  was  to  reduce  the  initial  concentrations  of 
hydrocarbons  from  a  total  of  26,000  ug/L  to  5000  ug/L. 
Benzene  and  toluene  are  preferentially  removed,  owing  to 
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their  higher  vapor  pressure  (95  torr  for  benzene  and  29  torr 
for  toluene  compared  with  6  torr  for  m-xylene  at  25  C).   The 
half  lives  for  biodegradation  under  these  conditions  were  in 
the  same  range  as  those  for  the  higher  concentrations  (i.e. 
treatment  2A).   Degradation  to  the  low  ug/L  level  was  noted 
for  all  compounds.   The  most  recalcitrant  was  o-xylene.   The 
addition  of  NH .CI  produced  the  same  effect  noted  in  other 
treatments.   The  lag  phase  seemed  to  be  especially 
significant  for  benzene,  o-xylene,  1 , 3 , 5-trimethylbenzene , 
2-ethyltoluene  and  1 , 2 , 3-trimethylbenzene .   As  in  previous 
experiments,  1 , 2 , 4-trimethylbenzene  was  well  degraded  under 
all  non-sterile  conditions.   The  reduced  treatment  of 
solutes  in  2H  was  accompanied  by  consumption  of  oxygen  and 
an  apparent  increase  in  microbial  activity.   The  DO  profiles 
show:   that  the  concentration  of  dissolved  oxygen  remained 
above  4  mg/L  throughout  the  entire  study,  indicating  no 
oxygen  limitation. 
5.10.5   Treatments  2C,  2F,  21  (Sterile  Controls) 

Sterility  in  treatments  2C,  2F  and  2G  was  indicated  by 
the  lack  of  0_  consumption,  low  INT  reduction  and  the 
persistence  of  aromatic  hydrocarbons.   The  losses  (assumed 
to  result  from  diffusion  through  the  teflon  septa)  were  in 
the  range  of  0-25%  for  C  -CQ  hydrocarbons  and  20-50%  for 
CgH,„  compounds.   Sorption  losses  were  accounted  for  in 
these  data  by  adding  the  amount  lost  to  sorption  to  the 
aqueous  concentration  data.   The  drop  in  DO  in  treatment  2G 
indicated  that  oxygen  also  diffuses  through  the  teflon 
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septa.   Based  on  this  loss,  the  decrease  in  DO  in  treatments 
2G  and  2H  was  not  attributed  solely  to  microbial  activity. 
The  diffusion  of  oxygen  across  the  teflon  septa  may  also 
account  for  the  increased  degradation  of  the  solutes  in  this 
study.   In  general/  two  parts  of  oxygen  are  required  to 
remove  one  part  of  hydrocarbon.   Using  these  calculations/ 
the  microcosms  should  have  been  oxygen  limited.   The  absence 
of  oxygen  limitation/  particularly  in  treatments  with  no 
added  source  of  oxygen,  was  likely  the  result  of  oxygen 
diffusion  into  the  microcosms.   This  suggests  that  in  future 
studies,  a  more  suitable  microcosm  should  be  employed. 
5.10.6   Discussion  of  Batch  Biodegradation  Data 

The  rates  of  biodegradation  (k  values  and  half-lives) 
were  highest  in  both  biodegradation  experiments  with  air 
augmentation  (1A  and  2A)  or  with  the  addition  of  oxygen 
(2G).   Almost  complete  degradation  (to  below  detection  limit 
or  less  than  2  ug/L)  was  shown  in  each  case.   These  data 
indicated  that  the  limiting  substance  was  oxygen.   Analysis 
of  water  quality  data  supports  this  hypothesis.   Total 
phosphate  in  the  groundwater  averages  0.6  rag/L  and  nitrate 
averages  0.25  mg/L . 

Addition  of  NH  Cl  to  the  microcosms  may  have  spurred 
the  process  of  nitrification.   This  would  account  for  the 
reduction  in  dissolved  oxygen  values  without  the  concomitant 
loss  of  hydrocarbons.   This  effect  was  seen  in  treatment  2B. 
In  this  case,  the  DO  dropped  from  7.5  to  2.5  mg/L  with  no 
loss  in  hydrocarbon  concentration.   The  presence  of  nitrate 
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in  the  aquifer  at  the  Lake  Alfred  site  and  the  low  dissolved 
oxygen  levels  in  the  aquifer  indicate  that  there  may  be 
significant  levels  of  nitrifying  bacteria  which  may  be 
stimulated  by  the  addition  of  NH  CI .   Indirect  evidence  for 
this  hypothesis  is  seen  in  the  results  of  the  INT  studies 
for  treatments  2A  and  2B.   Treatment  2B  showed  a  seven-fold 
increase  in  electron  transport  activity  although  there  was 
no  substantial  decrease  in  hydrocarbon  concentration  (Figure 
5-22). 

The  microorganisms  in  this  study  are  able  to  degrade 
aromatic  hydrocarbons  rapidly  down  to  the  1  ug/L  range  given 
sufficient  oxygen.   These  data  confirm  the  work  of  Jensen  et 
al .  (1985)  who  show  the  degradation  of  aromatic  hydrocarbons 
in  petroleum  contaminated  groundwater  to  1  ug/L  or  less. 
The  rates  of  biodegradation  determined  from  these  batch 
biodegradation  studies  were  significantly  faster  than  that 
of  Kappeler  and  Wuhrmann  (1978b).   In  that  study  benzene 
required  12  days  to  degrade  completely.   Delfino  and  Miles 
(1985)  showed  an  eight  day  lag  phase  for  benzene  dosed  to 
clean  groundwater.  In  this  study,  benzene  was  completely 
removed  in  a  minimum  of  two  days,  indicating  the  adaptation 
of  bacteria  from  the  Lake  Alfred  site  to  degrade  aromatic 
hydrocarbons . 
5.10.7   Hydrogen  Peroxide 

The  introduction  of  hydrogen  peroxide  into  the 
microcosms  produced  a  lag  phase.   The  microbial  populations 
required  several  days  to  adapt  to  the  changing  environmental 
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conditions.   This  was  reflected  in  the  electron  transport 

activity  graphs  where  H202  treatment  reduced  INT-formazan 

production.   This  may  result  from  the  toxicity  of  H  0   to 

the  microbial  population. 

The  addition  of  hydrogen  peroxide  to  the  batch 
microcosms  did  increase  the  DO,  especially  in  biodegradation 
experiment  #1.   Initial  DO  levels  were  increased  by  2  to  3 
mg/L  following  addition  of  H20   in  treatments  IB  and  1C. 
However,  hydrogen  peroxide  treatment  did  not  increase  the 
rates  of  degradation,  or  improve  the  extent  of  removal  of 
aromatic  hydrocarbons.  Air  or  oxygen  addition  was  shown  to 
be  more  effective  methods  of  oxygen  supplementation,  since 
they  did  not  require  an  adaption  period. 

This  conclusion  may  not  be  significant  under  field 
conditions,  where  an  adaption  period  may  not  be  an  important 
drawback,  given  the  economy  of  HO   relative  to  air  or 
oxygen  addition.  These  studies  demonstrate  that  H  0   does 
not  limit  the  extent  of  biodegradation,  but  that  there  is  an 
adaption  period  associated  with  its  use.   No  column 
experiments  were  performed  with  hydrogen  peroxide.   These 
experiments  would  have  assessed  hydrogen  peroxide  reactivity 
under  flowing  conditions. 


5.11   Column  Biodegradation  Experiments 
Several  authors  employed  flow-through  soil  columns  to 
study  the  degradation  of  organic  contaminants  (Kuhn  et  al . , 
1985).   Kuhn  et  al .  (1985)  noted  that  if  one  expects  to 
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apply  laboratory  derived  rate  constants  to  field  conditions 
that  input  concentrations  should  be  similar  to  field 
concentrations.   The  columns  used  in  this  work  were  packed 
with  soil  from  the  Lake  Alfred  field  site,  and  were  supplied 
with  well  water  from  Lake  Alfred  so  that  field  and 
laboratory  conditions  were  closely  matched.   Two  flow  rates 
were  used  to  better  simulate  the  varied  flow  conditions 
present  in  the  Lake  Alfred  aquifer. 

The  results  of  the  columns  run  at  low  flow  velocities 
(0.01  cm/min)  are  shown  in  Figure  5-23.   These  data  were 
analyzed  to  determine  rate  constants,  and  these  are  shown  in 
Table  5-22.  The  rate  constants  were  calculated  with  equation 
[4.7],   The  rates  of  degradation  were  much  higher  in  the 
column  system  than  in  the  batch  biodegradation  system.   The 
rate  constants  of  the  aromatic  compounds  in  the  column 
system  were  one  to  two  orders  of  magnitude  greater  than 
batch  study  constants.   The  increased  rate  of  hydrocarbon 
removal  was  the  result  of  improved  transport  of  carbon 
sources,  nutrients  and  oxygen  to  the  microbial  community  in 
the  column  system.   The  microcosms  in  the  batch  studies  were 
not  continuously  mixed  so  that  portions  of  the  microcosm  may 
have  been  nutrient  or  oxygen  limited.   Removal  efficiencies 
of  85-95%  were  seen  for  all  compounds  except  for  benzene, 
which  was  degraded  more  slowly  and  showed  a  60%  percent 
removal.   The  order  of  degradation  was  m,p-xylene  > 
ethylbenzene  >  o-xylene  >  toluene  >  benzene.   This  order  was 
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Table  5-22-   First  order  biological  rate  constants  and 

half-lives  of  aromatic  hydrocarbons  for  the 
biodegradation  column  with  flow  at  0.90  mL/hr, 


Compound  C/Co      ka   .         t 


day-1       day 


Benzene  0.36  5.80  0.120 

Toluene  0.144  11.01  0.063 

Ethylbenzene  0.086  13.94  0.050 

m,p-Xylene  0.051  16.82  0.041 

0-Xylene  0.128  11.69  0.059 

acalculated  from  the  following  data: 

length  =  2.5  cm 

bulk  density  =  1.8  g/mL 

particle  density  =  2.6  g/mL 

pore  water  velocity  =  14.21  cm/day 

volumetric  water  content  =  0.31 
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consistent  with  the  removal  of  side  chains  prior  to  attack 
on  the  aromatic  nucleous. 

Degradation  rates  for  columns  run  at  the  higher  flow 
rate  are  shown  in  Table  5-23.   The  faster  flow  rate 
decreased  the  half  lives  of  the  aromatic  contaminants.   This 
resulted  from  improved  transport  of  oxygen  and  subtrate . 
The  flow  rate  used  in  this  column  was  equivalent  to 
groundwater  velocities  in  portions  of  the  Lake  Alfred 
aquifer  (Killan,  1987)  and  may  closely  reflect  field 
removals.   Breakthrough  curves  for  benzene,  toluene  and 
1 , 2 ,4-trimethylbenzene  are  shown  in  Figures  5-24  ,5-25  ,  and 
5-26.   Only  10%  of  the  benzene  and  30%  of  the  toluene  were 
removed  at  this  flow  rate  although  other  removals  are  in  the 
range  of  50-60%.   The  difference  in  the  degradation  for 
benzene  is  seen  by  comparison  of  Figures  5-24  with  Figure  5- 
26.   Benzene  almost  breaks  through  the  column  completely, 
reflecting  the  time  involved  for  the  microbes  to  degrade 
this  solute.   The  branched  aromatic  compounds  are  more 
rapidly  degraded,  which  is  consistent  with  the  results  from 
the  0.680  cm/min  column.   The  rates  of  degradation  of  the 
aromatic  compounds  in  the  0.01  cm/min  columns  were  in  the 
order  Cg  >  Cg  >  C^    >    Cg_   Benzene  was  the  most  recalcitrant 
with  a  half  life  of  0.12  days  (2.88  hours).   These  data  are 
consistent  with  the  aromatic  degradation  process  described 
by  Evans  (1977),  and  the  literature  on  the  fate  of  aromatic 
hydrocarbons  in  soils  (Bossert  and  Bartha,  1984). 
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Table  5-23. 


First  order  biological  rate  constants  and 
half-lives  of  aromatic  hydrocarbons  for  the 
biodegradation  column  with  flow  at  1  mL/min 


Compound 


day 


-1 


■1/2 
day 


Benzene 

Toluene 

Ethylbenzene 

m/p-Xylene 

o-Xylene 

Isopropylbenzene 

n-Propylbenzene 

3  or  4  Ethyltoluene 

1,3,5-Trimethylbenzene 

2-Ethyltoluene 

1 , 2,4-Trimethylbenzene 

1 ,  2,  3-Trimethyl benzene 


17.7 
71.1 
125.3 
133.1 
115.2 
166.7 
192.4 
172.4 
178.1 
158.4 
172.3 
150.6 


0.040 
0.010 
0.006 
0.005 
0.006 
0.004 
0.004 
0.004 
0.004 
0.004 
0.004 
0.005 


Calculated   with    the    following    data: 

length  =  5.0  cm 

pore  water  velocity  =  0.680  cm/min 

particle  density  =  2.6  g/mL 

bulk  density  =  1.82  g/mL 

volumetric  water  content  -    0.30 
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Based  on  these  data/  it  is  evident  that  a  well  adapted, 
standing  microbial  population  from  Lake  Alfred  is  capable  of 
degrading  aromatic  hydrocarbons  at  relatively  high  loadings 
and  short  contact  times.   Degradation   may  be  aided  by  the 
development  of  an  efficient  biofilm  (Bouwer  and  McCarty, 
1984).   Benzene  was  degraded  to  a  lesser  extent  than  the 
alkyl  aromatics  and  o-xylene  was  more  resistant  to  microbial 
degradation  than  the  meta  and  para  isomers.   Kuhn  et  al . 
(1985)  in  column  studies  with  all  three  xylene  isomers  noted 
the  same  phenomenon.   The  increased  resistance  of  o-xylene 
was  also  seen  in  the  batch  biodegradation  data  presented  in 
this  study. 

In  the  batch  studies,  benzene  was  degraded  at  higher 
rates  than  most  of  the  branched  aromatic  compounds.   This  is 
contrary  to  the  column  data.   This  difference  in  removal 
rates  may  be  the  result  of  longer  contact  time  in  the  batch 
studies.   This  phenomenon  has  important  consequences  for  the 
degradation  of  hydrocarbons  in  field  studies,  where  high 
flow  rates  may  provide  insufficient  contact  time  for 
biological  removal.   Flow  rates  were  seen  to  significantly 
affect  the  rates  of  degradation  in  the  experimental  column 
system.   The  half  life  for  benzene  decreased  from  0.120  days 
at  0.01  cm/min  to  0.039  days  at  a  flow  velocity  of  .680 
cm/min.   The  half  life  for  toluene  decreased  from  0.063 
hours  to  0.0097  hours. 

However,  these  rates  are  derived  from  non-limiting 
conditions  and  these  data  may  not  completely  represent 
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biodegradation  in  the  field,  where  oxygen  limitation  reduces 
the  biological  removal  of  hydrocarbons. 

The  rapid  biodegradation  of  solutes  in  the  batch  and 
column  biodegradation  studies  supports  the  hypothesis  that 
the  microbial  community  in  the  Lake  Alfred  aquifer  is  well 
adapted  to  remove  aromatic  hydrocarbons  from  gasoline 
sources.   In  instances  where  adaptation  has  occurred/ 
biotransformations  may  be  so  rapid  that  they  are  considered 
instantaneous  relative  to  the  rate  of  groundwater  flow. 
This  shifts  the  quantitative  prediction  of  biological 
activity  from  a  consideration  of  biological  kinetics  to  a 
consideration  of  the  extent  of  utilization  (Wilson  et  al., 
1983b).   Thus  geochemical  constraints  become  the  controlling 
factors  in  the  biodegradation  process.   This  is  the 
situation  at  the  Lake  Alfred  site,  where  oxygen  is  the 
limiting  substance. 

5.12   Field  Data 

It  is  beyond  the  scope  of  this  work  to  thoroughly 
describe  all  the  field  data  collected  during  the  past  18 
months.   Rather,  the  major  findings  of  this  dissertation 
will  be  correlated  with  selected  aspects  of  the  field  data. 
These  topics  include  analysis  of  field-scale  dispersion, 
field-scale  solute  transport  and  the  enumeration  of 
microbial  populations. 
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5.12.1  Dispersion 

The  breakthrough  curve  for  the  NH  CI  tracer  at  the  Lake 

Alfred  site  is  shown  in  Figure  5-27.   The  dispersion 

2 
coefficient  calculated  from  these  data  was  0.546  cm  /min. 

Compared  with  the  dispersion  coefficients  from  the  column 

experiments  performed  at  0.680  cm/min  (0.01-0.07  cm  /min), 

the  field  scale  dispersion  is  an  order  of  magnitude  larger 

than  the  column  data.   The  calculated  P   is  approximately 

25.   This  value  should  be  useful  in  the  modeling  of  the  Lake 

Alfred  aquifer. 

5.12.2  Solute  Transport 

An  analysis  of  solute  transport  at  the  Lake  Alfred 
research  site  is  hampered  by  geologic  and  man-made 
obstacles.   The  presence  of  a  swale  running  through  the  site 
and  of  a  dual  flow  pattern  around  the  pump  house  (Building 
12)  have  produced  preferential  flow  in  the  aquifer.   Also, 
sewer,  drainage,  steam  and  telecommunications  lines  criss- 
cross contaminated  portions  of  the  aquifer,  complicating  the 
flow  path  of  dissolved  hydrocarbons.   It  is  also  probable 
that  gasoline  storage  tanks  in  between  Buildings  10  and  12 
and  transfer  equipment  south  of  the  wash  rack  added  unknown 
quantities  of  gasoline  to  portions  of  the  study  area 
(Killan,  1987).   This  skews  the  distribution  of  hydrocarbons 
in  the  aquifer  and  makes  determination  of  solute  transport 
difficult.   Finally,  the  pumping  wells  (UF-2M,  RAP-1  and 
RAP-3)   distort  the  transport  of  contaminants.   Less 
retarded  solutes  appear  to  be  more  retained  by  reversal  of 
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the  hydraulic  gradient.   Conversely,  solutes  which  are  more 
highly  sorbed  move  more  quickly  towards  the   pumping  wells 
as  a  result  of  increased  convective  flow  of  the  mobile  phase 
(ie.,  groundwater).   With  these  caveats  in  mind, 
interpretation  of  the  field  data  becomes  very  complex. 

The  distribution  of  benzene  at  the  field  site  is 
highest  toward  the  swamp  (well  UF-3W)  as  is  shown  in  Figure 
5-28.   The  major  concentrations  of  benzene  seem  to  have 
migrated  substantially  faster  than  the  other  compounds. 
This  is  expected  from  the  relatively  low  retardation  factor 
demonstrated  in  the  laboratory  studies,  and  also  from  the 
inefficient  biodegradation  of  benzene  at  high  flow,  as 
determined  in  the  biodegradation  column  run  at  0.680  cm/min. 
In  this  column,  only  10%  of  the  benzene  was  removed,  whereas 
the  branched  aromatics  were  more  easily  degraded. 

Hydrocarbon  data  from  each  monitoring  well  are 
presented  in  Appendix  G.   Examination  of  these  data 
demonstrate  that  the  areal  distribution  of  other  aromatic 
solutes  are  less  distinct  than  benzene.   Analyses  of  these 
data  do  not  yield  a  sequential  distribution  of  compounds 
suggested  by  retardation  factors  found  in  laboratory 
experiments.   This  results  from  the  increased  susceptibility 
to  microbial  attack,  and  the  probability  of  multiple  spill 
sites  . 

Ortho-xylene  (Figure  5-29)  has  a  measured  retardation 
factor  of  1.6,  which  is  20%  larger  than  that  of  benzene. 
However,  the  areal  distribution  of  o-xylene  is  much 
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UF-3W 


Wetland 


Figure    5-2! 


Distribution    of    benzene     (ug/L)     at    the 
Lake    Alfred    field    site. 
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UF-2M 


Wetland 


Figure  5-29.   Distribution  of  o-xylene  (ug/L)  at  the 
Lake  Alfred  field  site. 
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different  than  that  of  benzene.  Based  on  a  comparison  of  R 
values,  the  distribution  of  o-xylene  should  be  skewed 
towards  the  wetlands,  as  is  benzene.   One  factor  which  may 
account  for  the  apparent  increase  in  retardation  for  o- 
xylene  is  the  higher  rate  constant  for  biological  removal 
(half  life  =  0.144  hours  at  0.680  cm/min)  compared  to 
benzene  (0.940  hours).   This  results  in  improved  removal  of 
o-xylene  as  it  moves  through  the  aquifer,  so  that  the  solute 
front  appears  retarded,  relative  to  benzene.   These  type  of 
analysis,  although  somewhat  crude,  indicates  the  importance 
of  biological  data  in  the  interpretation  of  solute  transport 
under  field  conditions. 
5.12.3   Microbial  Enumeration  and  biodegradation 

Microbial  populations  at  the  field  site  are  given  in 
Tables  5-24  to  5-27.   Generally  they  range  from  105  to  106 
organisms/gram  dry  weight  soil.   Given  sufficient  oxygen, 
this  population  is  sufficient  for  complete  biodegradation  as 
shown  in  the  batch  laboratory  studies;  no  additional 
nutrients  are  necessary.   The  extent  of  oxygen  limitation  is 
noted  in  Table  5-28,  where  the  range  in  concentrations  of  DO 
over  the  study  period  are  presented.   Wells  which  exhibit 
low  DO  also  exhibit  high  concentrations  of  hydrocarbons.   An 
example  of  the  apparent  lack  of  significant  bioactivity  is 
seen  in  an  examination  of  well  data  from  OHM-3.   This  well 
shows  consistently  low  DO  (0.4  mg/L),  with  no  significant 
changes  in  the  concentrations  of  any  of  the  individual 
solutes . 
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Table  5-24.   Microbial  populations  in  a  soil  core 

taken  south  of  the  paint  shop  (Bldg.  54), 
June,  1986. 


a     5 
CFU/gdw  x  10 

Depth,         

feet  avg.    std.  dev. 


CFU/gdw:  colony  forming  units  per  gram  dry  weight 


Table  5-2  5.   Microbial  populations  in  a  soil  core 
taken  in  the  spray  field,  June,  1986. 


Depth, 
feet 

CFU/gdwax  105 

avg . 

std.  dev. 

0.5 

82.3 

5.18 

1.0 

30.8 

0.54 

1.5 

8.77 

0.66 

2.0 

4.21 

0.36 

3.0 

4.59 

2.83 

CFU/gdw:   colony  forming  units  per  gram  dry  weight 
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Table  5-26.   Microbial  populations  in  a  soil  core 

taken  south  of  the  pump  house  (Bldg.  12) , 
July,  1986. 


Depth, 

CFU/gdw 

x  106 

feet 

avg.    s 

td.  dev. 

Comments 

0.5 

3.6 

0.73 

1.0 

4.6 

1.21 

1.5 

6.0 

2.25 

2.0 

2.6 

0.36 

2.5 

2.1 

0.91 

4.0 

4.5 

2.16 

5.0 

1.6 

0.16 

6.0 

3.3 

0.60 

Saturated  zone 
gasoline  odor 
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Table  5-2  7.   Microbial  populations  from  samples 
collected  during  installation  of 
monitoring  wells  RAP-5  and  RAP-6, 
September,  1986. 


CFU/gdw  x 

106 

Depth, 

RAP-4 

RAP 

-5 

feet 

avg . 

std.  dev. 

avg. 

std.  dev. 

Comments 

2.0 

a 
ns 

ns 

4.9 

0.34 

3.0 

8.9 

0.96 

3.2 

0.25 

4.0 

4.1 

0.22 

3.1 

0.82 

5.0 

ns 

ns 

8.4 

0.31 

6.0 

4.2 

0.11 

2.0 

0.85 

Saturated 
zone 

no  sample 
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The  opposite  appears  to  be  true  in  data  from  well  P-7. 
This  well  shows  significant  decreases  in  some  aromatic 
hydrocarbons,  and  the  dissolved  oxygen  levels  start  to 
increase  after  November/  1986.   This  may  be  the  result  of 
the  increased  recharge  of  aerated  water  and  change  in 
pumping  conditions  established  by  Killan  (1987).   Decreases 
in  hydrocarbon  concentrations  concomitant  with  increasing  DO 
are  also  noted  up  gradient  of  P-7  prior  to  November  22/ 
1986.   Wells  RAP-6  and  RAP-5  exhibit  rapid  removals  of 
hydrocarbons  following  the  start  of  increased  flushing  with 
aerated  water  (22  October/  1986).   However/  the  data  are 
insufficient  to  conclude  whether  this  change  results  from 
the  increased  supply  of  oxygenated  water  and  subsequent 
biodegradation /  or  if  it  is  caused  by  the  more  rapid  solute 
transport  owing  to  increasing  the  flow  velocity  of  the 
aquifer . 


CHAPTER  VI 
SUMMARY  AND  CONCLUSIONS 


6.1   Sumraar1 


Hydrolysis,  sorption  and  biodegradation  reactions  of 

aromatic  hydrocarbons  (all  isomers  of  C-H--C-H,-)  under 

do   9  12 

water  saturated  soil  conditions  were  investigated.   Several 
treatment  techniques  (H202,  02  gas  and  NH4C1  )  were  evaluated 
to  enhance  the  microbial  degradation  of  the  aromatic 
compounds.   These  studies  were  performed  with  a 
multicomponent  solute  system  of  aromatic  hydrocarbons, 
resulting  from  the  partial  solubilization  of  gasoline  into 
groundwater,  obtained  from  a  field  site  where  a  gasoline 
spill  had  occurred.   The  sorbent  used  in  these  studies  was 
collected  from  the  field  site  in  a  non-contaminated  portion 
of  the  aquifer.   This  site  was  typical  of  sandy  surficial 
aquifers  in  Florida,  and  was  characterized  by  the  low 
organic  carbon  content  (0.015%)  of  the  aquifer  material. 
The  aquifer  was  composed  primarily  of  medium  to  fine  grained 
sands . 

Contaminated  well  water  from  the  surficial  aquifer  was 
employed  as  the  source  of  solutes  for  the  majority  of 
experiments.   The  major  components  of  this  water  were  the 
hydrocarbons  (C6H5-C9H12)  used  in  this  study. 
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6.1.1  Hydrolysis  studies 

Hydrolysis  did  not  account  for  substantial  losses  of 
the  aromatic  hydrocarbons  in  this  study.   These  solutes  were 
resistant  to  hydrolysis  even  under  extreme  (relative  to  the 
environment)  conditions  of  pH  (2,9,12)  or  temperature  (60 
O. 

6.1.2  Sorption  studies 

Sorption  of  Cg-C   aromatic  solutes  to  the  aquifer 
material  employed  in  this  work  was  relatively  rapid.  Rate 
studies  to  determine  the  approach  to  equilibrium  revealed 
that  there  was  an  initial  period  of  rapid  sorption  and  that 
equilibrium  conditions  were  established  within  four  to  eight 
hours . 

Multicomponent  sorption  experiments  of  the  dissolved 
aromatic  hydrocarbons  in  this  study  were  performed  in  batch 
isotherms  and  in  leaching  column  experiments.   Surficial 
well  water  was  used  as  the  source  of  these  solutes.   Aquifer 
material  from  the  Lake  Alfred  site  was  used  as  the  sorbent. 
Batch  sorption  employed  a  3:1  solids  to  solution  ratio  to 
approximate  aquifer  conditions,  and  to  maximize  the  change 
in  solution  concentration  resulting  from  sorption  of  the 
solutes  to  the  aquifer  material.   Equilibrium  batch  isotherm 
data  was  evaluated  with  the  Freundlich  model,  the  linear 
model  and  the  linear  model  with  suppressed  intercept. 
Sorption  coefficients  (Kd)  for  the  two  linear  models  were 
equivalent,  and  the  Freundlich  model  gave  similar  sorption 
coefficients  (K  ).   Kd  values  from  the  linear  models  ranged 
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from  0.49  for  toluene  to  0.142  for  1 , 3 , 5-tr imethylbenzene . 
It  was  concluded  that  the  sorption  process  was  reversible  in 
these  studies  with  no  significant  hysteresis. 

Column  sorption  experiments  were  performed  under 
saturated  steady  flow  conditions.   Sorption  coefficients 
from  batch  isotherm  studies  and  column  sorption  studies  were 

closely  matched.   The  average  dispersion  in  these  columns 

2 
was  0.044  cm  /min  with  a  standard  deviation  of  0.029 

2 
cm  /minute.   Retardation  factors  were  determined  for  all  C  - 

b 

C„  aromatic  compounds  from  the  multicomponent  system  of 
dissolved  aromatic  hydrocarbons  in  well  water.   Retardation 
factors  ranged  from  1.36  for  benzene  to  2.40  for  n- 
propylbenzene .   These  data  demonstrated  the  relatively  low 
retardation  of  aromatic  solutes  by  the  aquifer  materials. 
The  flow  velocity  used  in  the  column  experiments  (0.680 
cm/min)  was  close  to  seepage  velocities  measured  at  the 
field  site  (0.02-0.38  cm/min).   The  column  breakthrough 
curves  (BTC)  exhibited  some  non-equilibria  as  a  result  of 
slow  sorption  kinetics,  and  this  process  may  also  affect  the 
transport  of  solutes  in  the  field. 

The  influence  of  competing  solutes  was  investigated  by 
comparing  retardation  values  for  benzene  in  a  single  solute 
system  with  the  breakthrough  curve  for  benzene  in  the 
multicomponent  system.   The  retardation  factors  for  benzene 
in  both  column  systems  were  similar  (1.36  vs.  1.40).   Solute 
competition  for  sorbing  sites  was  not  a  significant  factor 
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in  this  study.   This  was  likely  the  result  of  the  use  of  low 

concentrations  ( less  than   2%  of  the  water  solubility). 

Sorption  mechanisms  were  evaluated  by  comparison  of  K 

oc 

data  from  column  studies  in  this  study  with  partitioning  and 

molecular  topology  models.  Regression  analysis  of  K    data 

oc 

versus  literature  values  for  K    demonstrated  that 

ow 

partitioning  alone  did  not  adequately  describe  the  sorption 

2 

process  (r   =  0.857).   Regression  analysis  of  K    data  with 

oc 

first  order  molecular  connectivity  indices  indicated  that 

sorption  may  be  partially  described  as  a  surface  area 

2 

dependent  phenomena  (r   =  0.839).   Both  models  gave 

equivalent  fit  to  the  K    data,  and  this  suggested  that 

oc  3S 

sorption  was  the  result  of  several  processes. 
6.1.3   Biodegradation  Studies 

Rate  constants  for  the  biodegradation  of  selected 
aromatic  hydrocarbons  were  determined  from  batch  and  column 
studies.   Column  studies  yielded  higher  rate  constants  than 
the  batch  studies  indicating  more  rapid  removal  of  solutes. 
Half  lives  for  columns  run  at  0.01  cm/min  were  between  0.120 
days  for  benzene  to  0.041  days  for  m,p-xylene.   Rate 
constants  derived  from  columns  with  velocities  of  0.680 
cm/min  were  higher,  indicating  increased  removal  of 
hydrocarbons.   These  half  lives  ranged  from  0.940  hours  for 
benzene  to  0.086  hours  for  n-propylbenzene .   The  increased 
removal  at  the  higher  flow  rate  was  the  result  of  improved 
transport  of  oxygen  and  nutrients  to  the  microbes.   The  1 
mL/min  flow  rate  (0.680  cm/min)  was  equivalent  to  seepage 
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velocities  at  the  field  site,  and  suggested  that  the  contact 
time  was  suitable  for  complete  degradation  of  aromatic 
solutes,  under  non-limiting  conditions.   This  indicated  the 
need  for  oxygen  augmentation  at  the  field  site  to  increase 
the  biodegradation  rates  of  aromatic  contaminants. 

Batch  biodegradation  experiments  were  performed  to 
assess  the  efficacy  of  various  methods  to  increase  the 
biological  degradation  of  dissolved  aromatic  hydrocarbons  at 
the  Lake  Alfred  field  site.   Laboratory  experiments  with 
hydrogen  peroxide  indicated  the  ability  of  the  microbial 
community  and  the  aquifer  materials  to  catalyze  the 
reduction  of  hydrogen  peroxide  to  yield  oxygen  gas.  No 
hydrocarbon  oxidation  was  apparent  as  a  result  of  the 
hydrogen  peroxide  decomposition. 

The  half  lives  for  biological  removal  of  the  selected 
aromatic  hydrocarbons  with  the  addition  of  air  (7-8  mg/L  0„) 
in  the  batch  biodegradation  studies  (Treatment  1A,  Thomas 
slope  rate  equation)  ranged  from  2.91  days  for  toluene  to 
4.96  days  for  1 ,  3 , 5-trimethylbenzene .   Benzene,  toluene  and 
1 , 2, 4-trimethylbenzene  were  degraded  most  rapidly. 
Augmentation  of  oxygen  in  the  form  of  air  or  oxygen  gas  was 
most  effective  for  increasing  the  biodegradation  of  aromatic 
hydrocarbons.   These  data  indicated  that  the  microbes  from 
the  Lake  Alfred  site  were  well  adapted  to  aromatic  gasoline 
hydrocarbons,  and  were  limited  only  by  the  availability  of 
oxygen.   Treatment  with  hydrogen  peroxide  (17  mg/L  to  68 
mg/L)  increased  the  dissolved  oxygen  levels  in  the 
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microcosms,  but  did  not  increase  the  rates  of  degradation  or 
the  extent  of  hydrocarbon  removal  in  batch  studies. 
Ammonium  chloride  addition  produced  conditions  favorable  to 
nitrifying  bacteria,  which  resulted  in  the  depletion  of 
oxygen  without  substantial  hydrocarbon  removal.   Treatment 
with  ammonium  chloride  resulted  in  a  lag  phase,  during  which 
the  microbial  community  adapted  to  the  changed  nutrient 
conditions.   The  combination  of  NH  cl  and  H„0   exhibited  a 
toxic  effect  to  the  microflora  from  the  Lake  Alfred  aquifer. 
This  suggests  that  the  use  of  NH  Cl  as  a  field  tracer  should 
not  be  followed  by  application  of  H_05  to  increase  dissolved 
oxygen.   Treatment  with  HO   at  concentrations  of  17  mg/L 
and  68  mg/L  produced  lag  times  of  up  to  eight  days,  with 
slightly  reduced  rate  constants. 

VOA  vials  were  employed  as  microcosms  for  the  batch 
biodegradation  experiments.   These  vials  allowed  substantial 
abiotic  losses  of  solutes,  as  determined  by  control  samples, 
and  were  not  recommended  for  use  in  future  experiments. 
6.1.4   Microbial  Field  Data 

Oxygen  was  the  limiting  substance  at  the  Lake  Alfred 
Site.   Dissolved  oxygen  concentrations  in  the  main  area  of 
hydrocarbon  contamination  was  only  0.1  to  0.4  mg/L.   The 
average  value  for  total  phosphorus  was  0.6  mg/L  and  the 
average  value  for  nitrate  was  0.3  mg/L  over  the  study  area. 
These  data  indicate  that  there  were  sufficient  phosphorus 
and  nitrogen  sources  to  initiate  the  biological  removal  of 
hydrocarbons,  given  sufficient  oxygen.   This  hypothesis  was 
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confirmed  by  the  laboratory  biodegradation  experiments/ 
where  the  microbes  from  the  field  site  were  well  adapted  and 
were  able  to  rapidly  degrade  the  aromatic  solutes  to  less 

than  0.5  ug/L.   Microbial  populations  at  the  field  site  were 

5       6 
determined  to  be  in  the  range  of  10   to  10   organisms  per 

gram  dry  weight  of  soil. 

6.1.5   Field  Scale  Solute  Transport 

The  distribution  of  contaminants  at  the  Lake  Alfred 

field  site  was  complicated  by  physical  structures/  numerous 

spills  and  underground  utilities.   This  made  the 

interpretation  of  field  data  difficult.   The  major 

concentrations  of  benzene  were  found  at  the  western  boundary 

of  the  field  site,  down  gradient  from  the  spill  area.   This 

distribution  of  benzene  demonstrates  that  the  combination  of 

low  retardation  and  inefficient  biodegradation  can  result  in 

increased  migration,  relative  to  solutes  which  were  more 

retarded  and  more  easily  degraded. 
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6.2   Conclusions 


1.  Hydrolysis  was  not  a  significant  removal  mechanism 
for  aromatic  hydrocarbons. 

2.  Solute-sorbent  equilibrium  was  established  in  batch 
sorption  vials  in  four  to  eight  hours. 

3.  Equilibrium  batch  sorption  isotherms  were  linear. 

4.  Column  breakthrough  curves  exhibited  apparent 
nonequilibria  . 

5.  Solute-solute  competition  for  sorbing  sites  was  not 
observed . 

6.  Equilibrium  batch  sorption  isotherm  data  and 
breakthrough  curve  data  yielded  similar  estimates  of  solute 
retardation  and  sorption. 

7.  Partitioning  models  and  the  first  order  molecular 
connectivity  model  gave  equivalent  fit  to  the  sorption  data. 
This  supports  the  hypothesis  that  sorption  results  from 
several  processes  depending  on  the  sorbent  and  the  solute. 

8.  Bacteria  from  the  Lake  Alfred  site  are  adapted  to 
aromatic  hydrocarbons  from  gasoline  sources. 

9.  Enzymatic  and  nonbiolog ical  hydrogen  peroxide 
catalysts  are  present  at  the  Lake  Alfred  site. 

10.   Hydrogen  peroxide  does  not  oxidize  aromatic 
hydrocarbons . 
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11.  Addition  of  air  or  oxygen  gas  was  the  most 
effective  method  for  stimulating  microbial  degradation  of 
the  solutes  in  this  study. 

12.  Hydrogen  peroxide  was  effective  in  increasing  the 
dissolved  oxygen  level  in  the  microcosms. 

13.  Hydrogen  peroxide  addition  to  the  microcosms  did 
not  increase  the  rate  or  extent  of  aromatic  hydrocarbon 
removal . 

14.  Ammonium  chloride  addition  to  the  microcosms  caused 
nitrification/  resulting  in  oxygen  consumption  in  the 
microcosms  without  hydrocarbon  removal. 

15.  The  combination  of  ammonium  chloride  (18  mg/L)  and 
hydrogen  peroxide  (17  mg/L)  additions  to  biodegradation 
microcosms  produced  toxic  conditions. 

16.  Column  biodegradation  studies  yielded  higher  rate 
constants  than  the  batch  studies/  reflecting  improved 
transport  of  nutrients  and  oxygen  to  bacteria. 

17.  Rates  of  biodegradation  for  aromatic  compounds  were 
in  the  order  C0  >  CQ  >  C_  >  Cc.   Benzene  was  the  most 
recalcitrant  solute  in  the  column  studies. 

18.  Benzene,  toluene,  m,p-xylene  and  1,2,4- 
trimethylbenzene  showed  the  most  rapid  biodegradation  in  the 
batch  studies. 

19.  Microbial  communities  in  the  Lake  Alfred  aquifer 

5    6 
were  in  the  range  of  10  -10   colony  forming  units  per  gram 

dry  weight  of  aquifer  material. 
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20.  Microbial  populations  at  the  Lake  Alfred  site  were 
oxygen  limited,  but  not  phosphorus  or  nitrogen  limited. 

21.  Horizontal  dispersion  at  the  field  site  was 
calculated  to  be  0.546  cm  /rain. 

22.  Low  retardation  (1.36),  high  flow  velocity  (0.680 
cm/min)  and  a  degradation  rate  of  17.65  ug/L/day  explain 
the  distribution  of  benzene  at  the  field  site. 
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APPENDIX  A 

CHROMATOGRAPHIC  CONDITIONS  AND  QUALITY  CONTROL 

PARAMETERS  FOR  THE  ANALYSIS  OF  AROMATIC  HYDROCARBONS 


This  appendix  lists  the  chromatographic  conditions  used 
in  the  gas  chromatographic  analysis  of  aromatic  hydrocarbons 
on  the  Perkin  Elmer  4100  gas  chromatograph.   Following  the 
GC  parameters,  summary  quality  control  data  is  presented  for 
chromatographic  analyses  performed  during  the  course  of 
these  studies. 
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METHOD    6        ANGLEY     2 


DATE     LAST     WRITTEN     87/05^12 


SECTION     1 


GC     CONTROL 


1 

2 

3 

4 

OUEN  TEMP  <DEG  C) 

50 

70 

34 

200 

ISO  TIME   CMIN> 

5.0 

7.0 

9.0 

0.0 

RAMP  RATE  <DEG  C/MIN) 

3.6 

3.0 

30.0 

FID   SENS        HIGH 
DET  ZERO  ON 

BET   TEMP  300 


FLOW 

A   55  ML/MIN 

CARRIER  GAS  HE 

EQUIL 

IB  TIME      0.5  MIN 

TOTAL 

RUN  TIME    30. 1  MIN 

SECTION 

2   TIMED  EUENT 

TIME 

EUENT 

0.01 

WIDTH          5 

8.02 

SET  ZERO 

2?. 80 

INTEG        OFF 

SECTION     3        DATA     HANDLING 
DATA    ACQUISITION 


START   TIME 
END    TIME 

WIDTH 
SKIM   SENS 
BASELINE   CORP 

AREA   SENS 
BASE    SENS 


0.00   MIN 
38.19   MIN 


100 

B-B 


us 


PEAK  IDENTIFICATION 


COMPONENT  LIST 


REFORT 

CALC  TYPE 

INT  STD 

AREA'HT  CALC 

AREA 

PRINT  TOL 

0.0000 

OUTPUT 

SCREEN 

NO 

PRINTER 

YES 

QUANTITATION  kCALIE  AUG  OF  12) 


UNRETD 
AREA/HT 

»EAK  TIME 
REJECT 

0.00 

0.0000 

MIN 

SCALING  FACTOR 
RF  FOR  UNKNOWNS 

1.0000 
10.0086 

REF  PK: 

TIME 
TIME  TOL 

13.19 
0.05 

MIN 
MIN 

STD  COMPNT  NAME 
SMP  AMOUNT 

CHLOROBN 

1.0000 

COMPNT: 

TOL  AES 
TOL  ;: 

0.05 
0.50 

STD  AMOUNT 

S2. 0000 

RT 


RF 


5 

S3 

7 

4916 

9 

68 

7 

5224 

13 

5b 

10 

0000 

14 

S2 

7 

8808 

15 

48 

7 

S3S0 

17 

25 

7 

??ee 

20 

23 

{ 

"~->~ir)r> 

22 

66 

7 

5404 

23 

38 

6 

764b 

24 

00 

g 

7483 

24 

70 

6 

7756 

2  j 

86 

7 

4867 

9"* 

44 

7 

9088 

D 

AMT 

NAME 

QR 

7 

1424 

BENZENE 

1 

4 

8768 

TOLUENE 

1 

62 

8000 

CHLOROBNZ 

8 

5 

4360 

ETHYL  BENZENE 

1 

5 

5038 

M.P  XYLENE 

1 

5 

632S 

0  XYLENE 

1 

5 

6840 

ISOPROPYL  BNZ 

2 

4 

1608 

N  PROPYL  EHZ 

2 

11 

8616 

3,4  ETHYLTOL 

2 

3 

2  528 

1.3.5  T  M  BNZ 

2 

5 

1312 

2,ETHYLTGLUENE 

-2 

4 

1664 

1,2.4  TI1BNZ 

5 

2448 

1,2,2  TMENZ 

2 
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Precision  and  accuracy  data  for  the  analysis  of  aromatic 
hydrocarbons  in  groundwater  by  EPA  method  6  02/ 
(modified) . 


Compound 

Benzene 

Toluene 

Ethylbenzene 

m,p-Xylene 

o-Xylene 

Isopropylbenzene 

n-Propyl benzene 

3 , 4-Ethy ltoluene 

1,3, 5-Trimethy lbenzene 

2-Ethyltoluene 

1,2, 4-Trimethy lbenzene 

1,2, 3-Trimethy lbenzene 


Precision 

Accuracy 

%  RSD, 

sd 

%  R, 

sd 

====== 

============ 

========= 

====== 

6.6, 

8.7 

98.0, 

13.4 

6.2, 

9.3 

99.3, 

18.0 

5.1, 

7.0 

95.2, 

12.7 

4.2, 

4.2 

95.3, 

12.1 

4.8, 

4.7 

100.8, 

15.9 

5.8, 

13.2 

89.8, 

11.3 

5.7, 

8.3 

86.5, 

10.7 

4.5, 

6.1 

90.4, 

10.2 

7.3, 

16.2 

89.6, 

12.7 

4.0, 

4.7 

83.8, 

10.2 

4.4, 

3.6 

90.0, 

15.2 

6.1, 

5.5 

93.0, 

12.2 

APPENDIX  B 
FIELD  SAMPLING  PROCEDURES 

The  sampling  procedures  employed  during  this  research 
were  described  in  the  Lake  Alfred  Quality  Assurance/Quality 
Control  (QA/QC)  plan  (Killan,  1987).   Section  six  of  the  QA/QC 
plan  is  presented  in  the  following  pages. 
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6.0   SAMPLING  PROCEDURES 

6.1  Cleaning  Procedures 

6.1.1  Volatile  Organics. 


Bottle  type:  -  water:  60  mL  glass  vial  with 

teflon  lined  septum  caps. 
-  soil:  1  quart  mason  jars. 

Soap:  Alconox 

1)  Wash  caps,  liners  and  vials  in  hot  soapy 

water  . 

2)  Rinse  liberally  with  tap  and  DI  water. 

3)  Rinse  with  pesticide  grade  methanol. 

4)  Dry  caps,  septa,  and  vials  in  oven  at  105 

C  for  no  more  than  60  minutes. 

5)  Cool  in  inverted  position,  and  cap 
immediately  when  bottles  are  cool  enough 
to  handle. 

6.1.2   Labels 

1)  After  cleaning  the  appropriate  label  is 
attached  to  each  bottle,  and  the  date  cleaned 
is  entered  . 

6.2  Field  Documents  and  Records. 

6.2.1  Field  sheets. 

The  field  sheet  (see  attachments)  is  filled 
in  with  the  following  information  upon 
sampl ing : 

1)  date 

2)  time 

3)  sample  type 

4)  preservation 

5)  well  number  (for  well  samples) 

6)  well  casing  and  diameter 

7)  depth  of  water  at  time  of  sampling 

8)  depth  of  core  (if  applicable;  soil 
samples) 

9)  note  special  characteristics  of 
sample . 

10)  field  number 

All  field  measurements  are  recorded  in  the 
bound  field  notebook  or  on  the  data  sheets. 
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All  samples  are  assigned  a  consecutive  field 
number  and  this  is  recorded  in  the  field 
notebook,  the  field  data  sheet  and  on  the 
sample  bottle. 

6.3  Water  Sample  collection  -  Field  procedures. 

6.3.1  Well  preparation. 

The  volume  of  water  in  each  well  is 
determined.  The  well  is  bailed  for  three 
times  the  calculated  volume.  This  is  done 
using  a  PVC  bailer,  or  a  battery  operated 
pump.  The  bailed  volume  is  measured  in  a 
calibrated  container.   Glass  tubing  is 
attached  to  the  end  of  tygon  tubing  to 
present  a  glass  surface  to  the  water  in  the 
well.  Tygon  tubing  is  attached  directly  to 
the  pump.   Each  well  is  supplied  with  a 
dedicated  piece  of  tygon  tubing  and  a 
dedicated  glass  insert. 


6.3.2  Volatile  Organic  Samples 

An  unopened  field  blank  is  taken  into  the 
field  . 

Wells  are  bailed  as  above. 

Vials  are  filled  by  lowering  the  container 

directly  into  the  well.  The  vial  is  filled  to 

overflowing.  The  teflon  side  of  the  septum  is 

placed  on  the  meniscus, and  capped  tightly. 

The  sample  is  inverted  and  examined  for  air 

bubbles . 

If  air  is  present  the  sample  is  discarded, 

and  the  well  resampled. 

Samples  are  taken  in  duplicate  and  chilled 

immediately  and  placed  in  the  dark.  Blanks 

are  stored  with  the  samples. 

No  preservation  of  the  samples  other  than 
cooling  is  performed. 


6.4  Soil  Sample  Collection 

6.4.1   Chemical  Analysis 

Soil  samples  destined  for  chemical  analysis 
for   gasoline  will  be  collected  with  a 
stainless  steel  auger.   Samples  will  be 
placed  into  a  one  quart  mason  jar.   Samples 
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will  be  tightly  packed  to  reduce  headspace  in 
the  sample  container.   Samples  are 
transported  on  ice  and  stored  at  4  C  in  the 
dark . 

6.4.2   Biological  Analysis 

Soil  samples  for  biological  analysis  will  be 
collected  with  sterilized  (100  ppm  chlorine 
solution)  stainless  steel  auger.  Samples  will 
be  placed  in  a  one  quart  mason  jar,  cleaned 
in  the  same  manner  as  the  described  in 
section  6.1.1.  Samples  will  be  stored  at  4  C 
in  the  dark. 


6.5   Measurement  of  field  parameters. 

6.5.1  Temperature. 

Temperature  will  be  measured  in  the  wells 
using  a  thermometer,  calibrated  against  an 
NBS  standard  thermometer. 

6.5.2  pH. 

The  pH  of  the  well  water  will  be  measured 
using  a  portable  pH  meter  (Orion  Research 
model  401)  connected  to  a  Fisher  AccupHast 
microprobe  combination  electrode. 

6.5.3  Dissolved  Oxygen. 

DO  will  be  measured  using  a  portable 
dissolved  oxygen  meter  (YSI  54A  )  with  a  YSI 
5739  oxygen  probe. 


APPENDIX  C 
ISOTHERM  DATA  FOR  THE  SORPTION  OF  STUDY  COMPOUNDS 
TO  LAKE  ALFRED  AQUIFER  MATERIAL 


Batch  isotherm  data  is  presented  for  each  compound 
in  this  study.   Freundlich  isotherms  for  each  compound  are 
also  shown. 
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APPENDIX  D 
BREAKTHROUGH  CURVE  DATA  FOR  THE  SORPTION  OF  STUDY 
COMPOUNDS  TO  LAKE  ALFRED  AQUIFER  MATERIAL 


This  appendix  presents  the  breakthrough  data  for  the 
aromatic  solutes  used  in  this  study.   The  data  is  presented 
as  concentrations  (ug/L)  and  as  C/C   (effluent  concentration/ 
influent  concentration) .   Following  these  data,  plots  of  C/C 
vs.  pore  volumes  are  presented  for  those  compounds  not  shown 
graphically  in  the  body  of  the  dissertation.   The  final  page 
of  this  appendix  presents  the  data  for  the  single  solute 
breakthrough  of  benzene. 
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Column  Sorption  Data 
L  =  5  cm 

v  =  0.2 0 4  cm/min 
pv  =  6.9  mL 

Values  are  as  ug/L 

mL     BNZ         TOL     ETHBZ   m,p-XYL     o-XYL 


1.76 

3 

4 

2 

2 

2 

3.74 

2 

3 

1 

2 

2 

4.71 

48 

15 

3 

4 

7 

5.71 

480 

158 

42 

37 

83 

6.7 

1332 

514 

197 

174 

335 

7.7 

2040 

825 

413 

313 

598 

8.7 

2777 

1200 

677 

525 

922 

9.69 

2948 

1345 

831 

655 

1115 

10.69 

3341 

1511 

951 

756 

1249 

12.69 

3877 

1826 

1238 

994 

1553 

13.69 

4081 

1935 

1343 

1070 

1658 

14.69 

4037 

1954 

1380 

1116 

1696 

15.69 

4332 

2124 

1482 

1188 

1831 

16.69 

4420 

2165 

1551 

1246 

1905 

17.69 

4589 

2189 

1641 

1315 

2013 

18.69 

4930 

2462 

1793 

1431 

2171 

19.69 

4531 

2368 

1801 

1472 

2183 

20.69 

4727 

2468 

1889 

1528 

2266 

21.69 

4871 

2505 

1883 

1534 

2267 

22.69 

4523 

2415 

1904 

1537 

2301 

23.69 

4615 

2377 

1802 

1460 

2210 

24.69 

4851 

2621 

2104 

1713 

2499 

26.69 

5454 

2985 

2385 

1952 

2905 

27.69 

5177 

2794 

2192 

1802 

2696 

28.69 

5563 

2999 

2386 

1947 

2897 

30.7 

5345 

2838 

2213 

1822 

2759 

31.7 

5197 

2  77  6 

2168 

1776 

2666 
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Column  Sorption  Data 
L  =  5  cm 

v  =  0.204  cm/min 
pv  =  6.9  mL 

Values  are  as  ug/L 

ispBZ      npBZ    3,4  ET   135  TMB      2  ET    124TMB    123TMB 


2  14  0  2  2  2 

113  0  2  13 

2  2  4  15  12 

5  2  7  2  9  4  12 

34  7  38  11  47  23  72 

90  22  101  33  113  61  168 

189  53  225  72  218  127  309 

271  85  359  112  314  188  431 

337  117  462  145  380  236  511 

509  198  695  229  542  357  721 

571  226  784  257  594  395  767 

616  254  838  282  637  431  821 

661  276  926  306  695  461  893 

706  302  990  331  730  499  941 

767  337  1109  368  799  555  1030 

845  370  1164  398  854  589  1085 

876  404  1281  431  894  636  1136 

922  420  1291  442  917  646  1152 

925  423  1306  448  932  654  1179 

937  439  1404  461  960  683  1224 

883  406  1270  435  909  642  1173 

1072  511  1561  532  1079  773  1350 

1229  579  1774  617  1270  907  1618 

1127  523  1586  562  1156  822  1478 

1233  581  1743  617  1258  897  1596 

1125  534  1662  580  1196  859  1552 

1121  539  1639  579  1178  845  1513 


Column  Sorption  Data  as  C/Co 
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PV 


BNZ 


TOL 


Ethyl  Bz  m,p-XYL   O-XYL 


0.253 

0.001 

0.002 

0.001 

0.001 

0.001 

0.537   . 

0.001 

0.001 

0.001 

0  .001 

0.001 

0.677 

0.010 

0.006 

0.002 

0.002 

0.003 

0.820 

0.101 

0.062 

0.023 

0.021 

0.037 

0.963 

0.280 

0.201 

0.108 

0.101 

0.150 

1.106 

0.429 

0.322 

0.226 

0.181 

0.269 

1.250 

0.584 

0.469 

0.370 

0.305 

0.414 

1.392 

0.620 

0.525 

0.455 

0.380 

0  .500 

1.536 

0.703 

0.590 

0.520 

0.439 

0.560 

•     1.823 

0.816 

0.713 

0.677 

0.577 

0.697 

1.967 

0.859 

0.756 

0.735 

0.621 

0.744 

2.111 

0.849 

0.763 

0.755 

0.648 

0.761 

2.254 

0.912 

0.830 

0.810 

0.690 

0.822 

2.398 

0.930 

0.846 

0.848 

0.723 

0.855 

2.542 

0.966 

0.855 

0.897 

0.763 

0.903 

2.685 

1.037 

0.962 

0.980 

0.830 

0.974 

2.829 

0.953 

0.925 

0.985 

0.854 

0.980 

2.973 

0.995 

0.964 

1.033 

0.887 

1.017 

3.116 

1.025 

0.979 

1.030 

0.890 

1.018 

3.260 

0.952 

0.944 

1.041 

0.892 

1.033 

3.404 

0.971 

0.929 

0.985 

0.848 

0.992 

3.547 

1.021 

1.024 

1.150 

0.994 

1.122 

3.835 

1.148 

1.166 

1.304 

1.133 

1.304 

3.978 

1.089 

1.091 

1.198 

1.046 

1.210 

4.122 

1.171 

1.171 

1.305 

1.130 

1.300 

4.411 

1.125 

1.109 

1.210 

1.057 

1.238 

4.555 

1.094 

1.084 

1.185 

1.031 

1.197 

Column  Sorption  Data  as  C/Co 
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ispBZ 


npBZ 


3,4  ET    135  TMB   2  ET 


124  TMB   123  TMB 


0.002 

0.002 

0.002 

0.000 

0.002 

0.002 

0.001 

0.001 

0.002 

0.002 

0  .000 

0.002 

0.001 

0.003 

0.001 

0.003 

0.002 

0.002 

0.005 

0.002 

0.002 

0  .005 

0.003 

0  .004 

0.003 

0  .009 

0  .005 

0.010 

0.034 

0.013 

0.023 

0.021 

0.047 

0.029 

0.059 

0.089 

0.038 

0.062 

0.061 

0.114 

0.079 

0.137 

0.188 

0.094 

0.138 

0.135 

0.219 

0.166 

0.253 

0.269 

0.151 

0.220 

0.211 

0.316 

0.245 

0.352 

0.335 

0.207 

0.283 

0.273 

0.382 

0.308 

0.418 

0.506 

0.350 

0.426 

0.431 

0.545 

0.466 

0.590 

0.568 

0.402 

0.481 

0.483 

0.598 

0.516 

0.627 

0.612 

0.450 

0.514 

0.531 

0.641 

0.563 

0.671 

0.657 

0.490 

0.568 

0.575 

0.699 

0.602 

0.730 

0.702 

0.535 

0.607 

0.623 

0.734 

0.651 

0.769 

0.762 

0.597 

0.681 

0.691 

0.804 

0.725 

0.842 

0.840 

0.656 

0.714 

0.747 

0.859 

0.769 

0.887 

0.871 

0.717 

0.786 

0.810 

0.899 

0.830 

0.929 

0.916 

0.744 

0.792 

0.831 

0.923 

0.843 

0.942 

0.919 

0.749 

0.801 

0.841 

0.938 

0.854 

0.964 

0.932 

0.777 

0.861 

0.867 

0.966 

0.892 

1.001 

0.878 

0.720 

0.779 

0.817 

0.914 

0.838 

0.959 

1.065 

0.906 

0.958 

1.000 

1.086 

1.009 

1.104 

1.222 

1.027 

1.088 

1.160 

1.278 

1.184 

1.323 

1.120 

0.927 

0.973 

1.056 

1.163 

1.073 

1.209 

1.226 

1.030 

1.069 

1.160 

1.266 

1.171 

1.305 

1.118 

0.947 

1.020 

1.090 

1.203 

1.121 

1.269 

1.114 

0.956 

1.006 

1.088 

1.185 

1.103 

1.237 
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Single  Solute  Breakthrough  Data 
For  Benzene 

L  =  5.0  cm 

Volume  =  24.5  cu.  cm. 

Bulk  Density  =  1.82  g/cu.  cm. 

1  PV  =  6.96  mL 

v  =  0.2  04  cm/min. 

Benzene 
mL     PV     (ug/L)    C/Co 


1.43 

0.205 

9 

0.009 

3.43 

0.493 

12 

0.013 

4.43 

0.636 

25 

0.027 

5.43 

0.780 

35 

0.038 

6.43 

0.924 

84 

0.092 

7.43 

1.068 

250 

0.272 

8.43 

1.211 

409 

0.446 

9.43 

1.355 

529 

0.576 

10.43 

1.499 

656 

0.714 

11.43 

1.642 

734 

0.800 

12.43 

1.786 

779 

0.848 

13.43 

1.930 

784 

0.854 

14.43 

2.073 

825 

0.899 

15.43 

2.217 

853 

0.929 

16.43 

2.361 

820 

0.893 

17.43 

2.504 

918 

1.000 

23.43 

3.366 

909 

0.990 

24.43 

3.510 

917 

0.999 

25.43 

3.654 

887 

0.966 

26.43 

3.797 

748 

0.815 

27.43 

3.941 

578 

0.629 

28.43 

4.085 

408 

0.445 

29.43 

4.228 

298 

0.325 

30.43 

4.372 

240 

0.261 

31.43 

4.516 

181 

0.197 

32.43 

4.659 

143 

0.156 

33.43 

4.803 

113 

0.123 

35.43 

5.091 

88 

0.095 

37.43 

5.378 

81 

0.089 

39.43 

5.665 

63 

0.068 

49.43 

7.102 

33 

0.036 

APPENDIX  E 
BATCH  BIODEGRADATION  DATA 


This  appendix  presents  data  from  batch  biodegradation 
experiments  one  and  two.   All  hydrocarbon  concentration 
values  are  in  units  of  ug/L.   Zero  values  indicate  that  the 
concentration  of  hydrocarbons  was  below  0.5  ug/L.   Dissolved 
oxygem  values  are  in  units  of  mg/L. 
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APPENDIX  F 
COLUMN  BREAKTHROUGH  DATA  FOR 
BIODEGRADATION  COLUMNS 


Breakthrough  data  for  columns  with  flow  rates  of 
1  mL/min  and   0.9  mL/hr   are  presented  in  tabular  form, 
Breakthrough  curves  for  each  compound  for  the  1  mL/min 
column   are  presented  following  the   tabular  data. 
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Column  Biodegradat ion  Data 

values  as  ug/L 
CUMM   ML      BNZ        TOL        EBZ        MPX         OX 


Co 

=  1711 

7034 

1061 

3848 

1931 

0 

0 

0 

0 

0 

0 

3.897 

72.2 

59.8 

36 

0 

17 

5.553 

177.4 

163.5 

50 

0 

6.975 

316 

272 

30 

33 

8.775 

433 

412 

34 

49 

9.828 

437 

485 

40 

64 

68 

10.8 

455 

499 

43 

69 

11.7 

490 

548 

54 

87 

13.05 

467 

590 

43 

62 

100 

14.4 

537 

731 

54 

93 

33 

15.75 

491 

690 

48 

92 

141 

17.1 

513 

726 

52 

110 

169 

18.225 

490 

751 

55 

119 

177 

22.725 

385 

649 

61 

133 

206 

25.425 

479 

818 

85 

182 

213 

28.125 

523 

853 

86 

181 

205 

29.25 

658 

1092 

93 

203 

249 

44.325 

599 

1085 

95 

192 

246 

50.85 

661 

1009 

88 

195 

259 

54.45 

644 

1033 

95 

227 

276 

L  =  2.5  cm 

v  =  0.003  cm/min 
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Column  Biodegradation  Data 

Values  are  C/Co 
TOL       EBZ       MPX 


PV 


BZ 


OX 


0.000 
0.570 
0.812 
1.020 
1.283 
1.437 
1.579 
1.711 
1.908 
2.105 
2.303 
2.500 
2.664 
3.323 
3.718 
4.112 
4.277 
6.481 
7.435 
7.962 


0.000 
0.042 
0.104 
0.185 
0.253 
0.255 
0.266 
0.286 
0.273 
0.314 
0.287 
0.300 
0.286 
0.225 
0.280 
0.306 
0.385 
0.350 
0.386 
0.376 


0.000 
0.009 
0.023 
0.039 
0.059 
0.069 
0.071 
0.078 
0.084 
0.104 
0.098 
0.103 
0.107 
0.092 
0.116 
0.121 
0.155 
0.154 
0.143 
0.147 


0.000 
0.034 
0.047 


0.038 
0.041 
0.051 
0.041 
0.051 
0.045 
0.049 
0.052 
0.057 
0.080 
0.081 
0.088 
0.090 
0.083 
0.090 


0.000 
0.000 
0.000 
0.008 
0.009 
0.017 


0.016 
0.024 
0.024 
0.029 
0.031 
0.035 
0.047 
0.047 
0.053 
0.050 
0.051 
0.059 


0.000 
0.009 
0.000 
0.017 
0.025 
0.035 
0.036 
0.045 
0.052 

0.073 
0.088 
0.092 
0.107 
0.110 
0.106 
0.129 
0.127 
0.134 
0.143 


L  =  2.5  cm 
v  =  0.003  crn/min 
pore  water  velocity- 
bulk  density  = 
particle  density  = 


14.21290 
1.8  g/ml 
2.6  g/ml 


cm/ day 
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L  =  5  cm 

v  =  0.204  cm/min 

1  pore  volume  =  6.94  ML 

Values  are  as  ug/L 

CUMM  ML       BNZ        TOL 


EBZ 


MPX 


O-XYL 


Co  =  5000 

2862 

2449 

1.78 

0.7 

2.9 

0.9 

3.78 

0.7 

6.4 

0.9 

4.78 

0.7 

9.2 

0.8 

5.8 

54.4 

12.9 

2.3 

6.8 

465.6 

122.3 

26 

7.8 

1346.1 

398.7 

122 

8.8 

2231 

769 

331.6 

9.8 

3104 

1186 

608 

10.8 

3472 

1373 

775 

11.8 

3822 

1552 

922 

13.8 

4206 

1740 

1082 

15.8 

4360 

1851 

1173 

16.8 

4401 

1846 

1170 

17.8 

4631 

1948 

1240 

18.8 

4414 

1921 

1264 

19.8 

5027 

2093 

1313 

20.8 

4774 

2106 

1373 

21.8 

4302 

1897 

1245 

22.8 

4577 

1968 

1266 

25.8 

4668 

2056 

1342 

26.8 

4642 

2035 

1320 

27.8 

4385 

1889 

1205 

28.8 

4632 

2011 

1307 

29.8 

4541 

2064 

1386 

30.8 

4378 

1809 

1080 

31.8 

4730 

2149 

1437 

32.8 

4579 

2065 

1383 

33.8 

4648 

2119 

1422 

34.8 

4914 

2241 

1513 

2092 


2805 


2.8 

3.4 

2.1 

1.9 

2.9 

2.5 

2.6 

4.9 

23.4 

50.4 

106 

211 

251 

484 

473 

817 

610 

1008 

728 

1170 

874 

1338 

951 

1428 

961 

1445 

1007 

1516 

1038 

1531 

1071 

1593 

1145 

1633 

1016 

1498 

1036 

1515 

1097 

1607 

108  7 

1578 

986 

1487 

1067 

1569 

1147 

1658 

900 

1371 

1193 

1710 

1134 

1636 

1173 

1681 

1237 

1780 

Pore  water  velocity  =  0.68  0  cm/min. 
Particle  density  =  2.6  g/mL . 
Bulk  density  =  1.8  2  g/mL. 
Volumetric  water  content  =  0.30 
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Values  are  as  ug/L 
ISPBZ     NPBZ    3,4ET   135TMB      2ET   124TMB   123TMB 

Co  =  151     850     2318      836     1479     1145     1786 


1.9 

1.6 

10.5 

1.4 

2.8 

2.2 

9.1 

0 

1.1 

5.7 

0.9 

1.9 

1.3 

5.3 

0 

0.6 

2.7 

0.8 

4.2 

1.5 

1.7 

6 

0.9 

7.3 

0 

1.6 

1.4 

1.9 

3.1 

2 

8.4 

2.2 

9.2 

2.9 

9.3 

21.8 

5 

27 

7.9 

34 

13 

44 

68 

16 

81 

26 

89 

45 

130 

165 

44 

204 

65 

194 

109 

267 

260 

78 

333 

107 

290 

170 

386 

359 

120 

462 

155 

382 

231 

501 

484 

193 

727 

228 

497 

335 

645 

557 

245 

774 

273 

563 

389 

712 

564 

255 

799 

284 

578 

409 

748 

607 

276 

860 

303 

612 

428 

776 

635 

302 

920 

324 

641 

463 

805 

658 

307 

943 

330 

654 

467 

819 

711 

353 

1048 

367 

710 

521 

899 

627 

306 

915 

323 

631 

460 

798 

638 

312 

933 

326 

636 

451 

800 

673 

333 

977 

348 

672 

491 

837 

667 

334 

926 

345 

668 

481 

838 

592 

289 

889 

313 

616 

445 

793 

659 

325 

1040 

340 

663 

464 

845 

709 

360 

1058 

373 

718 

529 

896 

522 

240 

789 

274 

549 

349 

702 

735 

364 

1111 

382 

736 

539 

908 

710 

355 

1057 

365 

702 

567 

871 

730 

368 

1106 

378 

725 

532 

895 

776 

393 

1144 

402 

767 

559 

948 
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Column  Biodegradat ion  Data 


C/Co 

C/Co 

C/Co 

C/Co 

C/Co 

PV 

BNZ 

TOL 

EBZ 

MPX 

o-XYL 

0.256 

0  .000 

0.001 

0  .000 

0.001 

0.001 

0.545 

0  .000 

0.002 

0  .000 

0.001 

0.001 

0.689 

0  .000 

0.003 

0  .000 

0.001 

0.001 

0.836 

0.011 

0.005 

0.001 

0.001 

0.002 

0.980 

0.093 

0.043 

0.011 

0.011 

0.018 

1.124 

0.269 

0.139 

0.050 

0.051 

0.075 

1.268 

0.446 

0.269 

0.135 

0.120 

0.173 

1.412 

0.621 

0.414 

0.248 

0.226 

0.291 

1.556 

0.694 

0.480 

0.316 

0.292 

0.359 

1.700 

0.764 

0.542 

0.376 

0.348 

0.417 

1.988 

0.841 

0.608 

0.442 

0.418 

0.477 

2.277 

0.872 

0.647 

0.479 

0.455 

0.509 

2.421 

0.880 

0.645 

0.478 

0.459 

0.515 

2.565 

0.926 

0.681 

0.506 

0.481 

0.540 

2.709 

0.883 

0.671 

0.516 

0.496 

0.546 

2.853 

1.005 

0.731 

0.536 

0.512 

0.568 

2.997 

0.955 

0.736 

0.561 

0.547 

0.582 

3.141 

0.860 

0.663 

0.508 

0.486 

0.534 

3.285 

0.915 

0.688 

0.517 

0.495 

0.540 

3.718 

0.933 

0.718 

0.548 

0.524 

0.573 

3.862 

0.928 

0.711 

0.539 

0.520 

0.563 

4.006 

0.877 

0.660 

0.492 

0.471 

0.530 

4.150 

0.926 

0.703 

0.534 

0.510 

0.559 

4.294 

0.908 

0.721 

0.566 

0.548 

0.591 

4.438 

0.875 

0.632 

0.441 

0.430 

0.489 

4.582 

0.946 

0.751 

0.587 

0.570 

0.610 

4.726 

0.916 

0.722 

0.565 

0  .542 

0.583 

4.870 

0.929 

0.740 

0.581 

0.561 

0.599 

5.014 

0.983 

0.783 

0.618 

0.591 

0.635 
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278 


C/Co 

C/Co 

C/Co 

C/Co 

C/Co 

C/Co 

C/Co 

ISPBZ 

NPBZ 

3,4ET 

135TMB 

2ET 

124TMB 

123TMB 

0.001 

0.002 

0  .005 

0.002 

0.002 

0.002 

0  .005 

0.000 

0.001 

0.002 

0.001 

0.001 

0.001 

0.003 

0.000 

0.001 

0.001 

0.001 

0.003 

0.001 

0.001 

0.004 

0.001 

0.003 

0.000 

0.001 

0.001 

0.001 

0.002 

0.002 

0  .004 

0.003 

0  .006 

0.003 

0.005 

0.014 

0  .006 

0.012 

0.009 

0.023 

0.011 

0.025 

0.045 

0.019 

0.035 

0.031 

0  .060 

0.039 

0.073 

0.109 

0.052 

0.088 

0.078 

0.131 

0.095 

0.149 

0.172 

0.092 

0.144 

0.128 

0.196 

0.148 

0.216 

0.237 

0.141 

0.199 

0.185 

0.258 

0.202 

0.280 

0.320 

0.227 

0.314 

0.273 

0.336 

0.292 

0.361 

0.368 

0.288 

0.334 

0.327 

0.381 

0.340 

0.399 

0.373 

0.300 

0.345 

0.340 

0.391 

0.357 

0.419 

0.401 

0.325 

0.371 

0.363 

0.414 

0.374 

0.434 

0.420 

0.355 

0.397 

0.388 

0.433 

0.404 

0.451 

0.435 

0.361 

0.407 

0.395 

0.442 

0.408 

0.458 

0.470 

0.415 

0.452 

0.439 

0.480 

0.455 

0.503 

0.414 

0.360 

0.395 

0.387 

0.427 

0.402 

0.447 

0.422 

0.367 

0.402 

0.390 

0.430 

0.394 

0.448 

0.445 

0.392 

0.421 

0.416 

0.454 

0.429 

0.469 

0.441 

0.393 

0.399 

0.413 

0.452 

0.420 

0.469 

0.391 

0.340 

0.383 

0.375 

0.416 

0.389 

0.444 

0.436 

0.382 

0.449 

0.407 

0.448 

0.405 

0.473 

0.469 

0.424 

0.456 

0.446 

0.485 

0.462 

0.502 

0.345 

0.282 

0.340 

0.328 

0.371 

0.305 

0.393 

0.486 

0.428 

0.479 

0.457 

0.498 

0.471 

0.508 

0.469 

0.418 

0.456 

0.437 

0.475 

0.495 

0.488 

0.483 

0.433 

0.477 

0.452 

0.490 

0.465 

0.501 

0.513 

0.463 

0.493 

0.481 

0.519 

0.488 

0.531 
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APPENDIX  G 
HYDROCARBON  CONCENTRATIONS  IN  MONITORING  WELLS  AT 
THE  LAKE  ALFRED  CITRUS  RESEARCH  AND  EDUCATION  CENTER 


All  concentration  values  are  in  units  of  ug/L.   Blank 
spaces  within  the  body  of  each  table  indicate  that  the 
concentration  was  below  the  limit  of  detection  (o.5  ug/L) 
of  the  analytical  system.   Figure  4-1  shows  the  location 
of  each  well. 
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